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ABSTRACT 


Pulsed  weather  radars  can  be  used  to  depict  meteorological  conditions  such  as 
rainfall  rate  and  wind  velocity.  These  quantities  are  calculated  from  measurements  of 
reflectivity,  mean  radial  velocity  and  velocity  spread  using  echo  signal  samples  from 
weather  targets.  These  radar  measurements  derive  from  modified  radio  frequency  (RF) 
echo  signal  parameters,  including  amplitude,  frequency  and  phase,  returned  to  the  radar 
from  the  weather  target.  RF  scattering  and  propagation  effects  modify  echo  signal 
parameters.  Bias  and  variance  in  the  weather  signal  parameter  estimates  naturally 
influence  the  accuracy  of  all  subsequent  quantities  produced.  For  meteorological 
products  to  be  as  accurate  as  possible,  the  amount  of  uncertainty  in  each  estimated 
quantity  must  be  minimized.  If  radar  system  parameters  are  not  accurately  known,  the 
reflectivity  estimate  will  be  biased.  A  well-controlled  calibration  process  is  therefore 
critical  to  reduce  the  bias  of  the  reflectivity  estimate.  This  thesis  presents  the  design  and 
implementation  of  one  such  calibration  system,  specifically  for  use  with  the  MWR-05XP 
(a  Mobile  Phased-Array  Pulse-Doppler  X-band  Weather  Radar  first  created  at  the  Naval 
Postgraduate  School  in  2005),  although  the  general  results  are  applicable  to  all  radars. 
The  calibration  system  presented  is  an  active,  external  calibrator  intended  to  verify  end- 
to-end  radar  system  performance.  Within  this  thesis,  a  background  of  the  radar 
calibration  problem  along  with  the  research  objectives  for  this  specific  project  and  related 
works  are  introduced  and  discussed.  The  theory  behind  the  operation  of  weather  radar 
(how  the  three  principle  quantities  are  measured  and  related  to  signal  parameters)  is  also 
presented  and  explained.  The  density  function  for  precipitation,  relation  between  signal 
correlation  and  velocity  spread,  and  fundamentals  of  weather  radar  signal  parameter 
estimation  are  given.  The  specific  design  is  described  in  explicit  detail  and  measured 
results  provided  in  appendices.  Conclusions  on  the  success  of  the  implementation  are 
drawn  and  recommendations  for  future  work  are  presented. 


V 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


VI 


TABLE  OF  CONTENTS 


1.  INTRODUCTION . 1 

A.  BACKGROUND . 1 

B.  RESEARCH  OBJECTIVE . 1 

C.  RELATED  WORK . 3 

D.  ORGANIZATION  OF  THE  THESIS . 5 

IL  BACKGROUND  AND  SUPPORTING  THEORY . 7 

A.  WEATHER  RADAR  VS.  TRADITIONAL  RADAR . 7 

B.  RADAR  RESOLUTION  CELL . 8 

C.  SWERLING  TARGET  TYPES . 9 

D.  SAMPLE  CORRELATION . 10 

E.  ESTIMATION  OF  POWER  DENSITY  AND  REFLECTIVITY  (Z) . 12 

1,  Antenna  Gain  and  Beam  effective  Solid  Angle . 12 

2,  Physics  of  Precipitation . 13 

a.  Size  and  Shape . 13 

b.  Drop  Size  Distribution  (DSD) . 14 

c.  Rainfall  Rate . 14 

3,  Weather  Radar  Range  Equation  and  Estimation  of  Reflectivity,,  15 

F.  RADAR  PRINCIPLES  AND  WEATHER  RADAR  EQUATION 

CORRECTION . 16 

G.  ESTIMATION  OF  VELOCITY . 19 

H.  ESTIMATION  OF  VELOCITY  SPREAD . 20 

I.  CALIBRATION  REQUIREMENTS . 20 

J.  CALIBRATION  METHODS  AND  SELECTED  METHOD . 21 

K.  SUMMARY . 22 

III.  DESIGN,  IMPLEMENTATION  AND  OPERATION . 23 

A,  DESIGN . 23 

B,  HARDWARE  IMPLEMENTATION . 24 

1,  System  Description  Overview . 24 

2,  Hardware  List . 25 

3,  Voltage  Regulator  Circuit . 25 

4,  RE  Components . 26 

5,  Control  Hardware . 26 

6,  Hardware  Limitations . 27 

C,  ACS  FREQUENCY/DOPPLER  CALIBRATION  AND  SOFTWARE 

IMPLEMENTATION . 28 

D,  ADDING  VELOCITY  SPREAD . 31 

E,  ACS  AMPLITUDE/REFLECTIVITY  CALIBRATION . 31 

1.  Measurement  Premise . 32 

2,  Measurement  Method . 33 

a.  Initial  Calculations . 34 

b.  Anechoic  Chamber  Measurements . 37 

vii 


IV.  CONCLUSION  AND  FUTURE  WORK . 47 

A.  CONCLUSION . 47 

B.  FUTURE  WORK . 48 

APPENDIX  A,  ACS  PICTURES . 53 

APPENDIX  B,  DIGITAL  PHASE  SHIFTER  DATA . 65 

APPENDIX  C.  SCALAR  AND  VECTOR  NETWORK  ANALYZER 

MEASUREMENT  DATA . 67 

APPENDIX  D,  PROGRAMMABLE  ATTENUATOR  DATA . 85 

APPENDIX  E.  ANTENNA  DATA . 87 

APPENDIX  F.  ANTENNA  PATTERN  DATA . 95 

A.  VNA  ANTENNA  PATTERN  DATA . 95 

B.  ACS  ANTENNA  PATTERN  DATA . 101 

C.  MATLAB  PLOTTING  FILES . 106 

D.  ABSOLUTE  GAIN  AND  3  DB  BEAMWIDTH  FOR  X-BAND  AND 

PASTERNACK  HORN  ANTENNAS  (CALCULATIONS  AND 
MEASUREMENTS) . 110 

E.  EXPLANATION  OF  DATAFILE  MEASUREMENTS  AND 

CALCULATIONS . 114 

F.  ANTENNA  ISOLATION  TEST  DATA . 116 

APPENDING.  ANECHOIC  CHAMBER . 119 

APPENDIX  H.  COMPACTRIO . 121 

A.  COMPACTRIO  REAL-TIME  EMBEDDED  CONTROLLER  (NI 

CRIO-9004,  64/512  MEMORY) . 121 

B.  COMPACTRIO  RECONFIGURABLE  EMBEDDED  CHASSIS  (NI 

CRIO-9101,  IM  GATE  4-SLOT  CHASSIS) . 124 

C.  COMPACTRIO  I/O  MODULES  FOR  DIGITAL  INPUT  AND 

COUNTER/TIMER  (NI  9401,  8  CH  5  V/TTL  HIGH  SPEED 
BIDIRECTIONAL  DIGITAL  I/O  MODULE) . 126 

D.  NI  DEVELOPER  SUITE  WITH  REAL-TIME  AND  FPGA 

SOFTWARE  FOR  COMPACTRIO . 128 

APPENDIX  I.  LABVIEW  FOR  THE  DIGITAL  PHASE  SHIFTER . 129 

A.  PHASESHIFTERCONTROL . 129 

B.  SLAVECONTROL . 130 

C.  MASTERCONTROL . 131 

LIST  OF  REFERENCES . 149 

INITIAL  DISTRIBUTION  LIST . 151 


viii 


LIST  OF  FIGURES 


Figure  1.  Illustration  of  Radar  Resolution  Cell.  [From  1] . 8 

Figure  2.  Bloek  Diagram  of  the  ACS . 23 

Figure  3.  Diagram  of  Component  Layout  inside  the  ACS  Enelosure . 24 

Figure  4.  Voltage  Regulator  Cireuit  Diagram . 26 

Figure  5.  Two  Cyeles  of  the  Periodie,  Linear,  Modulating  Sawtooth  Function . 29 

Figure  6.  Graph  for  Calculating  Radar  Cross  Section  of  a  Sphere.  [After  11] . 37 

Figure  7.  Time  vs.  Amplitude  for  Empty  Chamber,  Sphere  and  ACS.  (ESC_plotTA)  ..39 

Eigure  8.  Erequency  vs.  Amplitude  for  Empty  Chamber,  Sphere  and  ACS. 

(ESC_plotE) . 40 

Eigure  9.  Time  vs.  Amplitude  for  Sphere,  Zoomed.  (S_plotTA=Z) . 42 

Eigure  10.  Time  vs.  Amplitude  for  ACS,  Zoomed.  (C_plotTA=Z) . 43 

Eigure  1 1 .  Erequency  vs.  Equivalent  RCS  of  the  ACS.  (RCS_plot) . 44 

Eigure  12.  Calibrator  Eront  View  with  RAM  (straight  on) . 53 

Eigure  13.  Calibrator  Eront  View  with  RAM  (looking  down  angle) . 53 

Eigure  14.  Calibrator  Left  Side  View  (straight  on) . 53 

Eigure  15.  Calibrator  Eeft  Side  View  (looking  down  angle) . 54 

Eigure  16.  Calibrator  Eeft  Side  View  (looking  down  angle).  Zoomed . 55 

Eigure  17.  Calibrator  Eeft  Side  View  (looking  down  angle).  Zoomed  -  Top  Cover 

OEE . 56 

Eigure  18.  Calibrator  Top  View  (straight  on).  Zoomed  -  Top  Cover  OEE . 57 

Eigure  19.  Calibrator  Inside  View  (CRIO) . 57 

Eigure  20.  Calibrator  Inside  View  (Power  Supply  Circuit) . 58 

Eigure  21 .  Calibrator  Inside  View  (Phase  Shifter) . 58 

Eigure  22.  Calibrator  Inside  View  (Power  Supply) . 59 

Eigure  23.  Calibrator  Right  Side  View  (looking  down  angle)  Zoomed,  Top  Cover 

OEE . 60 

Eigure  24.  Calibrator  Right  Side  View  (looking  from  above  at  external  in-line  amp) . 61 

Eigure  25.  Calibrator  Rear  Eeft  Side  View  (looking  down  angle) . 61 

Eigure  26.  Right  Horn  Antenna  (close-in  view,  polarization  is  horizontal) . 62 

Eigure  27.  12”  Sphere  (used  to  calibrate  the  ACS) . 63 

Eigure  28.  Specifications  Sheet  for  Digital  Phase  Shifter  from  G.T.  Microwave . 66 

Eigure  29.  RE  Amplifier  Gain  Single  VNA  Plot . 67 

Eigure  30.  Digital  Phase  Shifter  Insertion  Eoss  (Eeft)  and  Phase  Value  (Right)  Dual 

VNA  Plots.  Zero  volts  applied  to  all  pins.  Active  Eow.  All  bits  “1”. 

(11111111) . 68 

Eigure  31.  Digital  Phase  Shifter  Insertion  Eoss  (Eeft)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  MSB  pin  only,  five  volts  others. 

(10000000) . 69 

Eigure  32.  Digital  Phase  Shifter  Insertion  Eoss  (Eeft)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  second  MSB  pin  only,  five  volts  others. 
(01000000) . 70 


IX 


Figure  33.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  third  MSB  pin  only,  five  volts  others. 

(00100000) . 71 

Figure  34.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  fourth  MSB  pin  only,  five  volts  others. 

(00010000) . 72 

Figure  35.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  fifth  MSB  pin  only,  five  volts  others. 

(00001000) . 73 

Figure  36.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  sixth  MSB  pin  only,  five  volts  others. 

(00000100) . 74 

Figure  37.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  seventh  MSB  pin  only,  five  volts  others. 

(00000010) . 75 

Figure  38.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  LSB  pin  only,  five  volts  others. 

(00000001) . 76 

Figure  39.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 

VNA  Plots.  Five  volts  applied  to  all  pins.  (00000000) . 77 

Figure  40.  Digital  Phase  Shifter  Forward  (left)  and  Reverse  (right)  Insertion  Loss 
Dual  VNA  Plots.  Demonstration  of  deviee’s  bidirectional  property. 
Losses  in  the  forward  and  reverse  are  equal.  Five  volts  applied  to  all  pins. 

(00000000) . 78 

Figure  41.  SNA  Plot  of  Scalar  Network  Analyzer’s  Straight-Through  State. 

“Calibrated  to  Zero” . 79 

Figure  42.  SNA  Plot  of  Cable  used  with  Digital  Phase  Shifter  during  SNA  testing . 80 

Figure  43.  Digital  Phase  Shifter  SNA  Plot.  Fixed  Phase  State,  no  Doppler  Shift . 81 

Figure  44.  Digital  Phase  Shifter  SNA  Plot.  Periodic  Rotation  through  Phase  States, 

Doppler  Shift . 82 

Figure  45.  Digital  Phase  Shifter  Phase  Value  Single  VNA  Plot.  Zero  volts  applied  to 

MSB  pin  only,  five  volts  others.  (10000000)  180°  Phase  Shift . 83 

Figure  46.  Digital  Phase  Shifter  Phase  Value  Single  VNA  Plot.  Zero  volts  applied  to 

second  MSB  pin  only,  five  volts  others.  (01000000)  90°  Phase  Shift . 84 

Figure  47.  Specifications  Sheet  for  Programmable  Attenuator  from  American 

Microwave . 86 

Figure  48.  Pasternack  Sectoral  Horn  Antenna  Data  Sheet . 88 

Figure  49.  Pasternack  Sectoral  Horn  Antenna  Outline  Drawing . 89 

Figure  50.  Pasternack  Waveguide  to  Coaxial  Adapter  Outline  Drawing . 90 

Figure  5 1 .  Narda  Standard  Gain  Horn  Antenna  Datasheet  (page  1) . 91 

Figure  52.  Narda  Standard  Gain  Horn  Antenna  Datasheet  (page  2) . 92 

Figure  53.  Narda  Standard  Gain  Horn  Antenna  Datasheet  (page  3) . 93 

Figure  54.  Narda  Standard  Gain  Horn  Antenna  Datasheet  (page  5) . 94 

Figure  55.  X-band  Horn  Antenna  #1  Pattern  at  9370MHz  (E  plane  in  red  and  H  plane 

in  blue) . 96 


X 


Figure  56. 
Figure  57. 
Figure  58. 

Figure  59. 

Figure  60. 

Figure  61. 

Figure  62. 

Figure  63. 

Figure  64. 

Figure  65. 

Figure  66. 

Figure  67. 

Figure  68. 

Figure  69. 

Figure  70. 

Figure  71. 

Figure  72. 

Figure  73. 

Figure  74. 
Figure  75. 
Figure  76. 
Figure  77. 
Figure  78. 

Figure  79. 
Figure  80. 
Figure  81. 
Figure  82. 


X-band  Horn  Antenna  #2  Pattern  at  9370MHz  (E  plane  in  red) . 

X-band  Horn  Antenna  #2  Pattern  at  9370MHz  (H  plane  in  blue) . 

X-band  Horn  Antenna  #1  Pattern  at  9680MHz  (E  plane  in  red  and  H  plane 

in  blue) . 

X-band  Horn  Antenna  #2  Pattern  at  9370MHz  (E  plane  in  red  and  H  plane 

in  blue) . 

X-band  Horn  Antenna  #1  Pattern  at  9990MHz  (E  plane  in  red  and  H  plane 

in  blue) . 

X-band  Horn  Antenna  #2  Pattern  at  9990MHz  (E  plane  in  red  and  H  plane 

in  blue) . 

Standard  Gain  Horn  Antenna  Pattern  at  9370MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Standard  Gain  Horn  Antenna  Pattern  at  9680MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Standard  Gain  Horn  Antenna  Pattern  at  9370MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Pasternaek  Horn  Antenna  #1  Pattern  at  9370MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Pasternack  Horn  Antenna  #2  Pattern  at  9370MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Pasternack  Horn  Antenna  #1  Pattern  at  9680MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Pasternack  Horn  Antenna  #2  Pattern  at  9680MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Pasternack  Horn  Antenna  #1  Pattern  at  9990MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Pasternack  Horn  Antenna  #2  Pattern  at  9990MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Standard  Gain  Horn  Antenna  Pattern  at  9370MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Standard  Gain  Horn  Antenna  Pattern  at  9680MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Standard  Gain  Horn  Antenna  Pattern  at  9990MHz  (E  plane  in  red  and  H 

plane  in  blue) . 

Plot  of  Absolute  Gain  for  X-band  Horn  Antennas . 

Plot  of  3dB  Beamwidth  for  X-band  Horn  Antennas . 

Plot  of  Absolute  Gain  for  Pasternack  Horn  Antennas . 

Plot  of  3dB  Beamwidth  for  Pasternack  Horn  Antennas . 

Narda  Standard  Gain  Horn  Absolute  Gain  Plot  (zoomed  to  area  of  interest 

and  annotated  with  values  of  interest) . 

Anechoic  Chamber  Diagram . 

Eigure  of  CompactRIO  Software  Architecture  from  ni.com . 

Picture  of  cRlO-9004  from  ni.com . 

Datasheet  for  CompactRIO  Real-Time  Embedded  Controllers  from 
ni.com . 


..96 

..97 

..97 

..98 

..98 

..99 

..99 

100 

100 

101 

102 

102 

103 

103 

104 

104 

105 

105 

111 

112 

113 

114 

115 
120 
121 
122 

123 


XI 


Figure  83.  Picture  of  cRIO-9103  from  ni.com . 124 

Figure  84.  Datasheet  for  CompactRIO  Reconfigurable  Embedded  Chassis  from 

ni.com . 125 

Figure  85.  Picture  of  cRlO-9401  from  ni.com . 126 

Figure  86.  Datasheet  for  C  Series  Digital  Output  and  Relay  Modules  from  ni.com . 127 

Figure  87.  datalO  CODE . 134 

Eigure  88.  datalO  UI . 135 

Eigure  89.  fpgaRX  test  CODE . 135 

Eigure  90.  fpgaRX  test  Ul . 136 

Eigure  91.  fpgaTX  test  CODE . 136 

Eigure  92.  fpgaXX  test  UI . 136 

Eigure  93.  KNORR-CAE_project . 137 

Eigure  94.  ManualControl  CODE . 138 

Eigure  95.  ManualControl  Ul . 139 

Eigure  96.  MasterControl  CODE . 139 

Eigure  97.  MasterControl_Ul . 140 

Eigure  98.  NetworkRX  test  CODE . 140 

Eigure  99.  NetworkRX_test_Ul . 141 

Eigure  100.  NetworkTX_test_CODE . 141 

Eigure  101.  PhaseShifterControl6-bit_CODEl . 142 

Eigure  102.  PhaseShifterControl6-bit_CODE2 . 142 

Eigure  103.  PhaseShifterControl6-bit_UI . 143 

Eigure  104.  PhaseShifterControlCODEl . 143 

Eigure  105.  PhaseShifterControl_CODE2 . 144 

Eigure  106.  PhaseShifterControlUl . 144 

Eigure  107.  SlaveControlCODEO . 145 

Eigure  108.  SlaveControlCODEl . 145 

Eigure  109.  SlaveControl_CODE2 . 146 

Eigure  110.  SlaveControlUI . 146 

Eigure  111.  tDcalc  CODE . 147 

Eigure  112.  tDcalc_Ul . 148 


xii 


LIST  OF  TABLES 


Table  1 .  Swerling  Target  Types  and  Probability  Density  Funetions . 10 

Table  2.  Summary  of  Constants  Used  in  Estimated  Equivalent  RCS  Caleulations . 35 

Table  3.  Summary  of  Estimated  Equivalent  RCS  Caleulations . 35 

Table  4.  Summary  for  Caleulations  of  Radar  Cross  Section  of  a  Sphere . 37 

Table  5.  Naming  Convention  Summary . 38 

Table  6.  Power  and  Eogic  Connections  for  Digitally  Controlled  Phase  Shifter . 65 

Table  7.  Power  and  Eogic  Connections  for  Programmable  Attenuator . 85 

Table  8.  Datafile  for  X-band  Horn  Antennas . 110 

Table  9.  Datafile  for  Pasternack  Horn  Antennas . 112 

Table  10.  Table  Summary  of  Eab VIEW  Vis . 133 


xiii 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


XIV 


EXECUTIVE  SUMMARY 


Pulsed  weather  radars  can  be  used  to  quantify  meteorology  into  amounts  such  as 
rainfall  rate  maps  and  wind  velocity  fields.  These  quantities  are  calculated  from 
measurements  of  reflectivity,  mean  radial  velocity  and  velocity  spread  using  echo  signal 
samples  from  weather  targets.  These  radar  measurements  derive  from  modified  radio 
frequency  (RF)  echo  signal  parameters,  including  amplitude,  frequency  and  phase, 
returned  to  the  radar  from  the  radar’s  target.  RF  scattering  and  propagation  effects  are 
the  mechanisms  which  modify  echo  signal  parameters.  Bias  and  variance  in  the  weather 
signal  parameter  estimates  naturally  influence  the  accuracy  of  all  subsequent  quantities 
produced.  For  meteorological  products  to  be  as  accurate  as  possible,  the  amount  of 
uncertainty  in  each  estimated  quantity  must  be  minimized.  If  radar  system  parameters  are 
not  accurately  known,  the  reflectivity  estimate  will  be  biased.  A  well-controlled 
calibration  process  is  therefore  critical  to  reduce  the  bias  of  the  reflectivity  estimate.  This 
thesis  presents  the  design  and  implementation  of  one  such  calibration  system,  specifically 
for  use  with  the  MWR-05XP  (a  Mobile  Phased- Array  Pulse-Doppler  X-band  Weather 
Radar  first  created  at  the  Naval  Postgraduate  School  in  2005),  although  the  general 
results  are  applicable  to  all  radars.  The  calibration  system  presented  is  an  active,  external 
calibrator  intended  to  verify  end-to-end  radar  system  performance. 

The  approach  is  a  repeater  that  will  produce  a  known  Radar  Cross  Section  (RCS) 
via  loop  gain  control  and  create  Doppler  velocities  and  Doppler  spread  via  phase 
modulation  with  a  sawtooth.  Proper  control  of  the  sawtooth  frequency  will  also  allow 
introduction  of  Doppler  spread.  Although  a  constant  amplitude  return  would  suffice,  use 
of  a  programmable,  constant  phase  attenuator  would  also  allow  control  of  the  amplitude 
statistics  of  the  return  to  the  radar.  Thus,  with  full  control  of  both  amplitude  and  phase, 
one  could  replicate  the  return  for  each  of  the  Swerling  target  types,  although  it  is  only 
precipitation  return  that  is  really  of  interest  in  our  application.  A  block  diagram  of  the 
implementation  of  this  proposed  Active  Calibrator  System  (ACS)  design  is  presented 
below. 
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The  reflectivity  for  a  weather  radar  contains  radar  system  parameters  which  may 
not  be  precisely  known.  The  objective  of  calibration  is  to  determine  the  composite  of 
these  parameters.  The  composite  of  these  factors  is  also  known  as  the  radar  system 
calibration  constant.  In  Dr.  Jeffrey  Knorr’s  paper,  “Weather  Radar  Equation  Correction 
for  Frequency  Agile  and  Phased  Array  Radars,”  it  has  been  shown  that,  for  the  case 
where  the  reference  target  is  on  boresight,  6  =  0,  the  measurement  is  made  at  the 
reference  frequency,  X  =  Xq,  and  the  range  is  short  {La  =  1)  the  radar  system  calibration 
constant  can  be  described  as: 


4  J 

1 

b 

_ 1 

From  this  equation,  it  is  clear  that  a  measurement  of  receiver  output  power  and 
range  for  a  reference  target  of  known  cross-section  {6  =  0,  X  =  Xo,  La  =  1)  permits  the 
radar  system  constant  on  the  left  side  to  be  determined  and  the  radar  calibrated  for 
reflectivity  estimation.  Phase  modulation  is  required  to  produce  Doppler  shift,  according 
to: 


v{t)- Acos'\joj  +  O{t)  =y4Re 


For  0(t)  periodic  where  the  phase  modulated  signal  has  a  spectrum  given  by  a  sawtooth 
waveform  of  amplitude  adjusted  to  produce  a  phase  change  that  is  a  multiple  of  In  (see 
waveform  pictured  below)  then, 

v(t)  =  where  (A^  =  m2;r) . 
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The  signal  is  translated  (Doppler  shifted)  from  frequency  to  frequency  . 


By  cycling  through  360°  (In)  of  phase  at  this  periodic  rate  a  known 
frequency/Doppler  shift  is  produced  against  which  the  weather  radar  can  be  calibrated. 
This  periodic  phase  cycling  (and  thus  frequency  translation)  is  implemented  in  the  ACS 
using  a  digital  phase  shifter.  Velocity  spread  can  be  added  to  the  Doppler  shifted  signal 
by  varying  the  frequency  of  the  sawtooth  waveform  used  in  the  phase  modulator.  The 
frequency  should  be  varied  in  accordance  with  a  Gaussian  random  variable  that  has  been 
low  pass  filtered  to  produce  a  signal  with  the  desired  correlation  time. 


On  the  basis  of  these  principles,  the  ACS  was  constructed  (RF  box  internal 
configuration  digital  photograph  follows): 


All  control  of  the  ACS  was  achieved  through  National  Instruments  Lab  VIEW 
(user  interface  shown  below). 
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After  construction,  measurements  were  made  to  determine  the  equivalent 
reflectivity  (RCS)  of  the  ACS.  The  result  of  those  measurements  and  that  entire  testing 


process  are  presented  below. 


Time  vs.  Amplitude  for  Empty  Chamber,  Sphere  and  ACS.  (ESC_plotTA) 
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Empty  Chamber,  Sphere,  &  Calibrator  -  Frequency  Domain 


Not  Gated  -  S 
Gated  -  S 
Not  Gated  -  C 
Gated  -  C 
Not  Gated  -  E 
Gated  for  Calibrator  -  E 
Gated  for  Sphere  -  E 


Frequency  vs.  Amplitude  for  Empty  Chamber,  Sphere  and  ACS.  (ESC_plotF) 

The  first  plot  is  a  composite  of  all  data  sets  in  the  time  domain,  not  zoomed  (all). 
The  second  plot  is  a  composite  of  all  data  sets  in  the  frequency  domain.  All  data  sets 
include  the  measured  data  for  the  three  targets:  an  Empty  Anechoic  Chamber  (chamber), 
the  12”  Aluminum  Sphere  (sphere)  and  the  ACS  (calibrator).  Timing  gates  were  placed 
around  the  sphere’s  location  (42  ns)  and  calibrator’s  location  (57  ns).  These  timing  gates 
eliminated  undesirable  echoes.  One  such  undesirable  echo  was  at  <2  ns.  This  echo 
corresponds  to  the  coupling  between  transmit  and  receive  antennas  which  were  placed 
about  18”  apart.  There  are  seven  graphed  data  sets  on  each  plot:  the  sphere  without  a 
time  gate  and  with  a  time  gate  centered  on  its  42  ns  return,  the  calibrator  without  a  time 
gate  and  with  a  time  gate  centered  on  its  57  ns  return,  the  empty  chamber  without  a  time 
gate  and  with  a  time  gate  for  the  calibrator,  then  with  a  time  gate  for  the  sphere. 


For  the  sphere,  the  time  gate  (G2)  was  centered  at  41.8  ns  and  spanned  39.8  to 
43.8  ns.  For  the  calibrator,  the  time  gate  (Gl)  was  centered  at  56.8  ns  and  spanned  54.8 
to  58.8  ns.  The  time  gates  correlate  directly  to  target  return  distances.  The  distance 
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(range)  from  transmit  and  receive  antennas  in  the  anechoic  chamber  to  the  target  pedestal 
inside  the  chamber  is  19  feet  (or  38.6  ns  in  time).  The  distance  is  traversed  twice  (once  to 
the  target  and  then  from  the  target)  at  the  speed  of  light.  The  reason  the  time  gate  for  the 
sphere  was  observed  at  and  consequently  set  at  42  ns  is  the  sphere  was  suspended  from 
the  ceiling  of  the  chamber  above  and  slightly  behind  the  target  pedestal  inside  the 
anechoic  chamber.  This  added  approximately  an  additional  1.6  feet  of  distance  and  3  ns 
of  time  because  the  hypotenuse  rather  than  the  longest  leg  of  the  triangle  was  traversed. 
The  reason  the  time  gate  for  the  calibrator  was  observed  at  and  consequently  set  at  57  ns 
is  different  than  for  the  sphere.  This  time  gate  was  approximately  8  ns  greater  because  of 
the  time  delay  due  to  transmission  line  lengths  rather  than  target  distance.  The  8  ns  is 
almost  8  feet  (one-way  path),  slightly  greater  than  the  straight  line  distance  between 
transmit  and  receive  antennas  of  the  ACS  (as  expected).  See  the  example  calculations 
below  for  the  mathematical  relation  between  time  and  distance  (two-way  path). 

c  =  3.00x10*  m/s 

i?  =  5.7912  m  (19  feet) 

9  /? 

—  =  38.6  ns 
c 

t  =  41.8  ns 

^  =  6.2700  m  (20.6  feet) 

Initially,  calculations  were  made  to  predict  the  received  power  from  both  the  ACS 
and  12”  Sphere.  The  expected  (calculated)  value  for  received  power  from  the  ACS  was 
-36.5dBm.  The  measured  value  for  received  power  from  the  ACS  was  -34.0dBm.  The 
expected  (calculated)  value  for  received  power  from  the  12”  Sphere  was  -47.9dBm.  The 
measured  value  for  received  power  from  the  12”  Sphere  was  -46.7dBm.  The  expected 
(calculated)  values  very  nearly  equaled  the  measured  values. 

The  next  plot  is  from  the  data  in  the  frequency  domain  of  both  the  sphere  and 
calibrator  previously  presented,  now  used  to  calculate  the  equivalent  RCS  of  the 
calibrator.  Measurements  were  made  from  9-10  GHz.  The  Vector  Network  Analyzer 
used  in  the  test  setup  as  the  signal  source  and  receiver  in  the  anechoic  chamber  recorded 
the  magnitude  of  the  returned  amplitude  power  in  logarithmic  form  (log  MAG)  in  units  of 
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milliwatts  (dBm).  All  measurements  are  in  the  frequeney  domain  gated  for  time  (FG). 
This  ineludes  both  the  sphere  (SFG)  and  ealibrator  (CFG).  The  measured  values  were 
very  close  to  the  expected  (calculated)  values.  The  received  power  for  the  12” 
Aluminum  Sphere  was  calculated  to  be  -47.8679  dBm  and  observed  to  be  near  -47  dBm. 
The  received  power  for  the  calibrator  was  calculated  to  be  -36.5032  dBm  and  observed  to 
be  near  -35  dBm.  The  RCS  of  the  sphere  was  previously  calculated  to  be  -1 1.3688  dBsm 
(S-RCS).  To  calculate  the  RCS  of  the  calibrator  (C-RCS)  based  on  the 
observed/measured  data,  the  CFG  -  SFG=A  (dB)  was  added  to  the  sphere’s  RCS.  Over 
the  entire  frequency  band  the  C-RCS  varied  from  0.6  to  3.2  dBsm.  Note  all  these 
measurements  in  the  chamber  included  the  use  of  the  RF  amplifier  within  the  calibrator 
system.  When  used  with  the  weather  radar,  signal  levels  are  great  enough  and  ranges 
short  enough  such  that  an  RF  amplifier  is  not  needed.  Furthermore,  calculations  are 
made  to  make  sure  the  signal  return  to  the  radar  is  in  the  proper  range  to  prevent  receiver 
saturation. 


Equivalent  RCS  of  Calibrator 
(with  in-line  amplifier) 


Frequency  (GHz) 


With  these  measured  values  the  equivalent  RCS  of  the  ACS  was  observed  to 
range  from  0.6  to  3.2  dBsm  across  the  9-10  GFlz  RF  band.  On  average  (mean),  the 
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equivalent  RCS  of  the  ACS  was  1.4934  dBsm.  The  median  was  1.5101  dBsm.  The  high 
value  was  3.1742  dBsm  at  9  GHz  and  the  low  value  was  0.5862  dBsm  at  9.38  GHz. 

An  aetive  ealibrator  system  (ACS)  for  weather  radar  was  sueeessfully  designed 
and  implemented.  The  ACS  was  verified  to  be  a  funetional  weather  radar  aetive 
ealibration  system  under  laboratory  testing.  It  sueeessfully  provided  the  means  (met  the 
requirements)  to  verify  all  three  of  the  signal  parameters  measured  (estimated)  by  a 
weather  radar.  These  inelude  the  zeroth  Doppler  moment  (refieetivity),  first  Doppler 
moment  (average  radial  veloeity),  and  seeond  Doppler  moment  (veloeity  spread). 
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I.  INTRODUCTION 


A,  BACKGROUND 

Pulsed  weather  radars  produee  estimates  of  refleetivity,  mean  radial  veloeity  and 
veloeity  spread  using  eeho  signal  samples  from  a  wide  variety  of  weather  targets.  From 
these  radar  return  signals,  other  meteorologieal  quantities  sueh  as  rainfall  rate  and  wind 
field  eharaeteristies  are  also  derived.  Estimates  are  derived  from  the  parameters  of  the 
modified  radio  frequeney  (RF)  eeho  signal  (amplitude  and  frequeney/phase)  seattered 
from  the  target  baek  in  the  direetion  of  the  radar.  RF  seattering  and  propagation  effeets 
are  the  meehanisms  whieh  modify  radar  eeho  signal  parameters  and  it  is  those 
modifieations  that  eontain  the  desired  weather-related  information  as  a  modulation  of  the 
earner  wave.  Bias  and  varianee  in  the  weather  signal  parameter  estimates  naturally 
influence  the  accuracy  of  all  subsequent  quantities  produced,  and  condition  assessments 
made.  For  meteorological  products  to  reflect  weather  conditions  as  accurately  as 
possible,  the  uncertainty  of  each  estimate  must  be  minimized,  or  accurately  characterized. 
If  radar  system  parameter  accuracies  are  not  precisely  known,  the  reflectivity  estimate 
will  be  biased  and  thus,  assessments  of  true  conditions  could  be  either  misleading  or 
masked.  A  well-controlled  calibration  process  is  therefore  critical  to  reduce  reflectivity 
bias  due  to  imprecise  knowledge  of  the  radar  system  parameters.  There  may  also  be  bias 
associated  with  the  statistical  estimation  process  and  its  resulting  estimator,  but  poor 
calibration  is  the  most  serious  source  of  error.  The  variance  of  the  estimate  of  interest 
will  depend  on  the  number  of  samples  and  their  correlation. 

B,  RESEARCH  OBJECTIVE 

The  research  objective  of  this  thesis  is  to  design  and  implement  a  radar  calibration 
system  (calibrator).  The  calibrator  is  for  specific  use  with  the  Naval  Postgraduate 
School’s  MWR-05XPi  X-band  weather  radar,  but  the  general  concept  is  applicable  to  all 

1  World’s  first  and  currently  (2007)  only  Mobile  Weather  Radar  created  in  2005  operating  at  X-band 
in  Pulse-Doppler  with  a  Phased- Array.  In  acronym,  this  is  the  MWR-05XP. 
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radars.  The  calibrator  presented  is  an  active,  external  calibrator  intended  to  verify  end-to- 
end  radar  system  performance.  An  ideal  active  calibration  system  produces  a  radar  return 
with  known  RCS  and  Doppler  shift  permitting  the  determination  of  the  radar  calibration 
constant.  As  of  2005,  the  MWR-05XP  operates  with  a  calibration  constant  determined 
through  the  use  of  radar  system  parameters  that  have  previously  been  only  approximately 
known.  A  known  RCS  return  condition  can  be  obtained  when  a  metal  sphere  is  used  to 
calibrate  a  radar.  This  calibration  method  involves  the  suspension  and  isolation  of  the 
sphere  from  great  height  at  long  distance^  in  ideal  weather  conditions.  These  pristine 
conditions  are  not  always  available  when  and  where  the  MWR-05XP  is  in  use.  In  fact, 
the  MWR-05XP  is  primarily  located  at  the  Center  for  Interdisciplinary  Remotely  Piloted 
Aircraft  Studies  (CIRPAS)  in  Marina,  CA  on  the  grounds  of  the  municipal  airport  (a 
relatively  small  workspace)  where  the  heights  and  distances  required  for  calibration  are 
never  available.  The  calibration  of  the  MWR-05XP  is  done  periodically  on  a  large  test 
range  at  Camp  Roberts,  CA  over  100  miles  away.  The  calibrator  research  involved  in 
this  thesis  effort  will  eliminate  the  prohibitive  burdens  (physical  conditions  and 
workload)  required  with  flying  the  sphere  while  simultaneously  allowing  more  specific 
calibration  parameter  determinations  along  with  the  inherent  advantage  of  more  frequent 
calibration  constant  verification.  The  calibrator  also  allows  the  added  calibration  of 
Doppler-related  quantities  (which  is  not  currently  possible).  The  final  advantage  is  the 
end-to-end  method.  It  uses  the  entire  radar  system  in  normal  operating  mode. 
Calibration  systems  sometimes  rely  on  piecemeal  measurements  that  are  combined 
together  for  performance  estimates,  thus  leading  to  larger  uncertainties.  This  system  does 
not. 

Research  in  this  effort  will  implement  an  existing  design  of  an  active  calibrator 
for  the  MWR-05XP  and  will  not  involve  the  creation  of  an  innovative  design 
configuration.  It  will  neither  involve  design  optimization,  nor  investigation  of 
alternatives.  This  calibrator,  once  implemented,  will  be  used  with  the  MWR-05XP 
exclusively  (i.e.  no  other  radars  involved)  to  determine  the  radar  calibration  constant. 

2  Sphere  heights  of  0- 1000m  with  ground  distance  equal  to  700m  have  been  used  with  the  MWR-05XP 
to  cover  beam  elevation  angles  from  0-55°. 
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The  type  of  signal  used  by  the  ealibrator  in  this  proeess  will  be  of  the  minimum  required 
eomplexity,  although  the  performanee  parameters  involved  are  those  that  most  aecurately 
reflect  operational  weather  characteristics. 

The  primary  hardware  involved  with  this  thesis  effort  will  be  a  digital  phase 
shifter,  programmable  attenuator,  and  real-time  controller  with  reconfigurable  chassis. 
The  only  software  involved  will  be  Lab  VIEW  (in  the  manipulation  of  the  controller 
hardware)  and  MATLAB  (in  the  analysis  of  data  and  display  of  results). 

C.  RELATED  WORK 

The  thesis  by  Schroder  presents  a  weather  radar  signal  simulator.  According  to 
the  abstract  of  the  Schroder  thesis, 

the  simulator  was  developed  in  MATLAB  and  implemented  several 
different  functionalities  allowing  for  stepped  frequency,  multiple  pulse 
repetition  frequencies  (PRFs),  pulse  compression  using  a  chirp  waveform, 
and  variation  of  both  weather  and  radar  input  parameters.  Post-detection 
processing  capabilities  include  autocorrelation  and  Fast  Fourier  Transform 
(FFT  -  for  single  PRF  only);  estimation  of  weather  parameters  such  as 
reflectivity  (Z);  average  Doppler,  radial  velocity,  and  velocity  spread; 
pedagogical  plots  including  a  phasor  plot  of  phase  change  over  time  and  a 
velocity  histogram,  instantaneous  observed  reflectivity  and  power  for  each 
pulse  over  time.[l] 

The  output  of  Schroder’s  simulator  was  modified  to  provide  the  calibrator  controller 
software  with  Doppler  frequencies  and  amplitude  levels  indicative  of  real-world  weather 
effects  for  use  in  the  calibration  process. 

The  article  by  Marshall  and  Hitschfeld  investigates  the  random  fluctuations  of  the 
radar  echo  produced  from  a  random  array  of  scatterers  (such  as  weather  scatterers). 
According  to  the  introduction  to  this  article,  “the  intent  is  to  establish  (theoretically  and 
by  analogue  computation)  the  extent  of  the  limitation  these  fluctuations  have  on  the 
accuracy  of  the  measurement  information  obtainable.  Additionally,  what  improvements 
may  be  effected  by  averaging  or  by  modifying  the  method  of  pulsing  or  scanning  are 
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explored.”[2]  The  Marshall  and  Hitschfeld  theory  is  used  in  this  researeh  to  produee  the 
amplitude  levels  indieative  of  real-world  weather  effeets  for  use  in  the  ealibration 
proeess. 

The  report  by  Pratte  and  lee  explores  external  radar  ealibration  methods. 
Aeeording  to  the  exeeutive  summary  of  the  Pratte  and  lee  report,  “the  targeted  radar  for 
this  report  is  the  WSR-88D,  the  weather  radar  in  use  (as  of  2007)  by  the  U.S.  National 
Weather  Serviee  in  its  national  weather  radar  network.  It  reviews  availability  and 
feasibility  of  “ealibration  aids”,  equipment  and  software  that  will  assist  in  verifying 
WSR-88D  refleetivity  measurement  and  in  redueing  its  uneertainty.”[3]  The 
reeommendations  of  the  Pratte  and  lee  report  were  used  to  enhanee  understanding  of  the 
general  weather  radar  ealibration  problem  and  its  alternative  solutions. 

The  text  by  Doviak  and  Zrine  emphasizes  the  applieation  of  Doppler  radar  for  the 
observations  of  stormy  and  elear  air.  Aeeording  to  the  introduetion  (1.2)  to  the  text,  “it 
presents  a  eomprehensive  treatment  of  the  teehniques  used  in  extraeting  meteorologieal 
information  from  weather  eehoes  and  relates  radar  and  signal  eharaeteristies  to 
meteorologieal  patterns. ”[4]  The  Doviak  and  Zrine  text  served  as  a  elassie  referenee 
book  on  radar  theory  and  aeeepted  teehniques  applied  to  meteorology. 

The  paper  by  Knorr  presents  a  derivation  of  two  eorreetion  faetors  to  the  elassie 
Probert- Jones  weather  radar  equation  for  use  with  frequeney  agile,  phased  array  weather 
radars.  Aeeording  to  the  abstraet  of  this  paper,  “the  result  is  an  extended  weather  radar 
equation  for  use  with  these  advaneed  weather  radars.  The  two  derived  faetors  aeeount  for 
the  effeet  of  a  ehanging  RF  earner  and  eleetronie  beaming  steering  not  present  in  elassie 
weather  radars  when  eomputing  a  refleetivity  estimate  from  weather  signal  data 
samples. ”[5]  The  Knorr  paper,  as  well  as  personal  diseussions  with  its  author  (an  advisor 
to  this  thesis),  provided  the  neeessary  understanding  of  the  unique  ealibration 
requirements  of  these  advaneed  weather  radars  from  theory  to  praetieal  applieation. 

Chapter  Four  of  the  Federal  Meteorologieal  Handbook  diseusses  potential 
operational  uses  of  the  WSR-88D  produets  as  well  as  the  strueture  of  eertain 
meteorologieal  phenomena.  Aeeording  to  the  introduetion  to  Chapter  4  (4.1),  “it 
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provides  information  on  observing  and  foreeasting  meteorologieal  phenomena  using  this 
radar’s  products. ”[6]  Of  particular  interest  in  the  development  of  this  thesis’s  research 
was  the  discussion  in  early  sections  of  operationally  troublesome  radar  returns  due  to 
such  things  as  ground  clutter  and  sidelobes.  The  pictorial  representations  in  the 
Meteorological  Handbook  were  invaluable. 

D,  ORGANIZATION  OF  THE  THESIS 

There  are  four  main  chapters  and  several  appendices  which  comprise  this  thesis. 
Chapter  I  gives  a  background  of  the  radar  calibration  problem  along  with  the  research 
objectives  for  this  specific  project.  Related  works  are  introduced  and  discussed.  An 
executive  summary  is  also  included.  Chapter  II  explains  the  theory  behind  the  operation 
of  weather  radar  (how  the  three  principle  quantities  of  reflectivity,  velocity  and  velocity 
spread  are  measured  and  related  to  signal  parameters).  The  density  function  for 
precipitation,  relation  between  signal  correlation  and  velocity  spread,  and  fundamentals 
of  weather  radar  signal  parameter  estimation  are  given.  Chapter  III  describes  the  specific 
calibrator  design  in  explicit  detail.  Measured  results  are  presented  and  explained  in 
appendices.  Chapter  IV  draws  conclusions  on  the  success  of  the  implementation  and 
provides  recommendations  for  future  work. 
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II.  BACKGROUND  AND  SUPPORTING  THEORY 


This  chapter  presents  and  explains  the  theory  behind  the  operation  of  weather 
radar.  The  discussion  includes  how  the  three  principle  quantities  of  reflectivity 3,  mean 
radial  veloeity^  and  veloeity  spread^  are  measured  and  are  related  to  radar  signal 
parameters.  The  density  function  for  precipitation,  relation  between  signal  eorrelation 
and  velocity  spread,  and  fundamentals  of  weather  radar  signal  parameter  estimation  are 
given. 

Significant  portions  have  been  excerpted  from  Chapter  2  “Weather  Radar  and 
Weather/Rain  Theory”  of  Ulf  Schroder’s  thesis  “Development  of  a  Weather  Radar  Signal 
Simulator  to  Examine  Sampling  Rates  and  Scanning  Schemes.”  Both  Schroder’s  thesis 
and  this  thesis  are  in  support  of  the  larger  researeh  project/interest  of  Dr.  Jeffrey  Knorr 
(advisor  to  both  theses).  Much  of  the  background  and  supporting  theory  used  in 
Schroder’s  thesis  (eompleted  2005)  applies  to  this  thesis  (completed  2008)  and  the  larger 
researeh  project/thrust. 

A,  WEATHER  RADAR  VS.  TRADITIONAL  RADAR 

Weather  radar  samples  the  scattered  RF  energy  from  distributed  targets  such  as 
precipitation,  hail,  snow,  melting  iee,  and  clear  air.  This  is  unlike  traditional  radar  which 
views  these  returns  as  clutter,  unwanted  echoes  that  interfere  with  the  observation  of  the 
return  from  a  point  target  of  interest.  Traditional  radar  usually  implements  Doppler 
filtering  to  diseriminate  between  slow  moving  weather  elutter  and  fast  moving  point 
targets.  In  this  application,  the  scattered  energy  contains  descriptive  characteristics  of  the 
weather  condition. 


3  Function  of  the  scatterer  intensity. 

Function  of  the  mean  component  of  scatterer  motion  in  the  radial  direction  from  the  radar. 

3  Function  of  the  variability  of  Doppler  velocity  values  due  to  turbulence  within  and  velocity  shear 
across  the  sampling  volume. 
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B, 


RADAR  RESOLUTION  CELL 


Weather  radar  differs  from  traditional  radar  in  another  important  way.  Traditional 
radar  seeks  to  eneapsulate  its  target  of  interest  (approximately  a  point  target)  entirely 
within  a  single  radar  resolution  cell.  For  weather  radar  this  is  neither  practical  nor 
informative  at  either  the  micro  or  macro  scales.  It  is  impractical  to  image  a  single  rain 
drop  or  an  entire  storm  and  even  if  it  were  possible  the  information  provided  would  not 
be  especially  useful.  Instead,  weather  radars  image  a  volumetric  segment  of  the 
meteorological  phenomena  of  interest.  This  volume  is  defined  by  a  segment  of 
propagated  energy  (see  Figure  1)  covering  pulse  width,  x.  The  depth  of  this  segment 
depends  on  the  pulse  width  according  to: 


Figure  1.  Illustration  of  Radar  Resolution  Cell.  [From  1] 

Schroder’s  figure  above  illustrates  and  later  expressions  (specifically,  reflectivity 
factor  Z)  use  the  volume  of  the  resolution  cell.  It  is  uniquely  defined  here  for 
convenience,  as: 
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(2.2) 


See  Chapter  II,  Section  E,  Subsection  I,  “Antenna  Gain  and  Beam  effective  Solid 

Angle,”  for  an  explanation  of  the  last  factor,  — — . 

2  In  (2) 

C.  SWERLING  TARGET  TYPES 

Meteorological  conditions  produce  a  radar  target  return  that  varies  with  time. 
Besides  the  target  itself,  fluctuations  in  the  target’s  return  may  also  be  caused  by 
variations  in  the  medium  between  the  target  and  the  radar  (atmospheric  and 
meteorological  conditions  not  the  intended  target),  radar  system  instabilities  (platform 
motion  and  equipment  instabilities),  target  aspect  changes,  and  many  other  variables.  In 
order  to  analyze  all  of  this  variation,  the  target  is  modeled  as  an  overall  system.  For 
systems  analysis  purposes,  only  the  ‘gross’  behavior  of  a  target  needs  to  be  known,  not 
every  detail  of  all  the  complicated  physics  involved  in  the  complex  scattering  situation. 
Therefore,  fluctuations  caused  by  all  of  these  many  sources  (exclusive  of  those  caused 
simply  by  a  change  in  range)  are  usually  grouped  together  and  represented  by  a  single 
quantity  for  a  specific  target  type. 

The  following  excerpt  was  relayed  from  the  radar  text  by  Skolnik,  although  nearly 
every  radar  text  contains  some  explanation  of  Swerling  target  types.  Let  cr  be  a  random 
variable  with  a  probability  density  function  (PDF)  that  depends  on  the  factors  associated 
with  target  fluctuations  other  than  a  change  in  range  and  the  overbar  represents  the 
average  value  of  the  random  variable.  Two  PDFs  are  commonly  used; 

p{cr)  =  Le-^l^  ,VTX#\  (2.3) 

cr 

and 

p{cT)  =  ,  PDF  #2  (2.4) 

cr 

These  two  PDFs  approximate  the  nature  of  the  scattering  of  RF  energy  from  many 
targets.  There  are  targets  which  consist  of  many  independent  scattering  elements  of 
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which  no  single  one  (or  few)  dominates  (PDF  #1)  and  targets  whieh  have  one  main 
seattering  element  that  dominates  in  the  presenee  of  smaller  independent  seattering 
sourees  (PDF  #2).  These  two  PDFs  along  with  the  rate  at  whieh  the  fluetuations  take 
plaee  allow  the  definition  of  four  types  of  radar  targets,  also  known  as  Swerling  target 
types.  Swerling  target  types  I  and  II  are  assoeiated  with  PDF  #1  and  types  III  and  IV 
with  PDF  #2.  The  rate  of  the  fluetuations  in  this  seattering  is  the  other  differentiating 
faetor  whieh  further  defines  the  Swerling  target  type.  There  are  targets  with  slow 
fluetuations  (sean-to-sean)  and  targets  with  rapid  fluetuations  (pulse-to-pulse).  Swerling 
target  types  I  and  III  are  of  the  former  and  types  II  and  IV  of  the  latter. 


Swerling  Target  Types  and  Probability  Density  Functions 

PDF#1 

PDF  #2 

Slow  Fluetuations 
(sean-to-sean) 

Type  I 

Type  III 

Rapid  Fluetuations 
(pulse-to-pulse) 

Type  II 

Type  IV 

Table  1 .  Swerling  Target  Types  and  Probability  Density  Funetions. 

As  a  weather  radar,  the  MWR-05XP  targets  Swerling  types  I  and  II.  The 
ealibrator  deseribed  in  this  thesis  was  eonsequently  developed  for  those  same  Swerling 
types,  but  the  general  proeess  of  the  design  and  implementation  of  a  radar  ealibrator  ean 
be  applied  to  any  radar  and  Swerling  target  not  just  weather  radars  and  Swerling  types  I 
and  II. 

D,  SAMPLE  CORRELATION 

To  paraphrase  Sehroder,  the  power  level  of  seattered  RF  energy  from  distributed 
targets  varies  over  time.  An  estimate  of  the  average  power  level  is  aehieved  by  averaging 
power  over  a  number  samples.  The  varianee  of  a  refleetivity  estimate  ean  be  improved 
simply  by  averaging  more  samples.  Flowever,  these  samples  may  not  be  independent  of 
one  another.  Independenee  signifies  a  suffieient  level  of  ehange  in  the  distribution  of  the 
target  within  the  volume  being  imaged.  For  this  ‘suffieient’  level  of  ehange  to  have 
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occurred  there  exists  a  minimum  required  period  of  time.  This  minimum  time  period 
guarantees  samples  are  not  eorrelated  from  one  observation  period  to  the  next. 


Schroder  explains,  “The  power  spectral  density  of  a  meteorological  signal  is 
approximately  Gaussian  and  can  be  written  as 


jf-ff 

2al 


(2.5) 


S{f)  =  S,e 

where  f  is  the  Doppler  frequency  and  cr^  is  the  Doppler  spectrum  variance.  Taking  the 
Fourier  transform  of  Equation  (2.5)  and  normalizing  results  in  the  correlation  coefficient 


p{t)i 


where  is  the  time  spectrum  variance  and  is  related  to  cr^  by 

1 


cr,  = 


2n(7 , 


and 


^/  = 


2cr„ 


(2.6) 


(2.7) 


(2.8) 


By  using  the  correlation  coefficient  a  measure  of  minimum  difference  in  time  can 
be  estimated  ensuring  dependent  or  independent  samples.  Nathanson  elaims  that  for 
independent  sampling  /7(r)<0.02  is  required  whereas  /7(r)>0.15  will  ensure 

dependenee.  Doviak  and  Zrnic  set  a  correlation  threshold  for  coherence  at  /7^(r)>e  ' 
which  corresponds  to  a  value  of  p(t)>  0.6  .  Using  this  threshold  renders 

Z 


T;  <0.25- 


2ra„ 


(2.9) 


For  highly  eorrelated  samples." 

A  level  of  ‘sufficient’  ehange  has  oeeurred  within  an  observation  period  for  a 
sample  so  that  any  subsequent  samples  may  be  deemed  independent  if: 

Z 


T  >0.70- 


n(j„ 


(2.10) 
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where  is  the  sample  time  period  and  <j^  is  the  RMS  veloeity  spread. 

E,  ESTIMATION  OF  POWER  DENSITY  AND  REFLECTIVITY  (Z) 


1.  Antenna  Gain  and  Beam  effective  Solid  Angle 


Antenna  gain  varies  depending  on  azimuth  and  elevation  angles  for  direetive 
antennas.  For  a  point  target,  antenna  gain  is  eonstant  only  for  a  fixed  azimuth  and 
elevation  angle  defined  between  the  radar  and  a  point  in  spaee  at  a  speeific  time.  For  a 
volumetrie  target,  antenna  gain  varies  aeross  a  range  of  azimuth  and  elevation  angles 
defined  between  the  radar  and  a  volume  of  spaee  at  a  speeifie  time.  The  antenna  pattern 
funetion,  f{6,<j)),  integrated  over  the  solid  angle  eorresponding  to  the  radar  resolution 
eell  is  used  to  estimate  a  value  of  antenna  gain  for  volumetrie  targets. 


From  Sehroder  (and  Probert-Jones  as  the  original  souree),  “The  magnitude  of  the 
antenna  pattern  funetion  is  first  squared  to  yield  a  power  relation  and  squared  again  to 
agree  with  the  two  way  propagation  of  radar. 

\\f{e,<l>)\dCL  (2.11) 

Att 


The  Gaussian  Pattern  Model  is  the  antenna  pattern  funetion  that  will  be  used.  It 
provides  a  good  approximation  of  the  main  beam  between  the  half  power  points. 
Extending  it  in  both  azimuth  and  elevation  provides 


=  \e\<U,  \<i>\<U 


(2.12) 


where  =  — —  and 


41n2 


41n2 


and  6*j,  (zij  are  the  half  power  beam  widths  in  the 


orthogonal  planes. 


Performing  integration  over  the  solid  angle  gives 

+00  +00 

\\f{04)\dn=\\e 


le^ 

1 

Ak 


ded(t>  =  ^^{\-A) 

81n2^  ’ 


(2.13) 
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where  ^  <  0.034  and  represents  the  eontributions  from  the  sidelobes.  In  a  further 
approximation,  at  its  maximum  value  the  (l-^)  term  eontributes  a  negligible  0.15  dB 
and  therefore  ean  be  ignored.” 


The  final  form  of  the  expression  used  to  approximate  beam  solid  angle  for 
volumetrie  targets  is 


4;7 


dQ.  = 


n6^(j)^ 

81n2 


(2.14) 


See  also  the  Probert- Jones  paper[7].  The  equation  above  gives  beam  solid  angle 
eontributing  to  the  volume  return.  Note  this  result  is  based  on  a  Gaussian  approximation 
to  the  beam  pattern.  As  Dr.  Knorr  shows  in  his  paper,  integration  of  an  aetual  antenna 
pattern  will  generally  result  in  a  eorreetion  whieh  may  not  be  negligible. 


2.  Physics  of  Precipitation 

The  following  presents  some  basies  in  preeipitation  physies  relevant  to 
developing  a  physieal  understanding  of  preeipitation,  speoifieally  rain  the  primary  target 
of  weather  radar.  The  following  three  seetions  have  been  seleetively  ehosen  from  a 
larger  diseussion  of  the  physies  of  preeipitation  by  Sehroder.  Eaeh  of  the  three  seetions’ 
text  has  been  paraphrased  and  equations  reprodueed  verbatim. 


a.  Size  and  Shape 

Raindrop  size  determines  the  terminal  veloeity  at  whieh  the  drop  falls  and 
the  shape  of  the  drop.  Small  drops,  with  a  diameter  D  <  0.35  mm,  are  spherieal.  Medium 
size  drops,  with  diameters  0.35  mm  <  D  <  4  mm,  have  progressively  flattened  bases  and 
are  approximated  as  spheroids.  Large  drops,  D  >  4  mm,  do  not  last  very  long  and  tend  to 
break  apart  upon  eollision.  The  largest  drops,  D  >  10  mm,  are  unstable  and  break  apart 
without  eollision. 
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b.  Drop  Size  Distribution  (DSD) 


Radar  cross  section  (RCS)  density,  rj ,  is  the  average  RCS  density  per  unit 
volume,  as  previously  mentioned,  and  for  independent  seattering  is  the  sum  of  the  RCS  of 
all  seatters  in  the  return  volume  divided  by  that  volume.  The  DSD  represents  the 
diameter  density  (mm/m  )  of  drops  for  all  diameters.  The  Marshall-Palmer  DSD 
provides  a  general  model: 

N{D)  =  N,e-^ 

A  =  4.1R  “"W  '  (2.15) 

=  8  X 1 0^  m^^mm  ' 

Applying  truneation  with  a  maximum  drop  diameter  of  yields 


/V(D) 


\N,e-^,D<D^^ 

lo, 


(2.16) 


c.  Rainfall  Rate 


Rainfall  rate  determines  the  depth  of  aeeumulated  water  per  unit  time.  It 

ean  be  derived  from  water  eontent  and  fall  speeds.  A  eloud’s  water  density  is 

00 

M  =  '^p„\D^N{D)dD  (2.17) 

where  is  the  water  density  of  a  rain  drop. 

Turning  this  into  rainfall  rate 

00 

R  =  -\D^N{D)v[D)dD .  (2.18) 

Rainfall  rate  is  estimated  from  refleetivity  using  any  one  of  a  number  of 
slightly  different  models,  Z  =  Ar^.  Z  =  300r'  for  example,  is  the  eurrent  WSR-88D 
default  model.  The  best  rainfall  model  depends  upon  the  type  of  storm  or  equivalently, 
drop  size  distribution.  Dual  polarization  radars  provide  more  information  about  drop  size 
and  thus  rainfall  rate.  Latest  results  from  JPOL  experiment  (dual  polarization  WSR-88D 
at  NSSL)  are  given  in  Ryzhkov[8]. 
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3. 


Weather  Radar  Range  Equation  and  Estimation  of  Reflectivity 


Again,  following  Schroder,  the  basie  form  of  the  radar  range  equation  (negleeting 
minor  losses): 


R  = 


{47rf  R" 


(2.19) 


“In  the  basie  form  of  the  weather  radar  equation  radar  cross  section  (RCS),  a,  is 
replaced  by  RCS  density,  rj ,  which  is  the  expeeted  RCS  density  per  unit  volume.  Sinee 
the  target  is  a  volume  target,  the  expeeted  return  needs  to  be  integrated  over  the  whole 
volume,  whieh  ean  be  divided  into  the  solid  angle  and  depth  of  the  return.  Letting  the 
power  level  at  the  target  be 

and  integrating  the  return  power  over  the  return  volume  renders 


(4^ 


ffujl/lejpa 

R  ^  Q 


1024;z-"  R^\n2 


(2.21) 


recognizing  that  R  ? 


CT 

~2' 


Radar  meteorologists  use  the  refleetivity  faetor,  Z,  instead  of  RCS  density  where 

7  =  ^K,rz  (2.22) 

and 

(2-23) 

^  ieV 

Substituting  Equation  (2.22)  into  Equation  (2.21)  and  solving  for  Z  yields  the  reflectivity 
estimator 


1024;z-^RMn2  Z" 


—  1024RMn2Z" 

^  RG,G^CT0,^y\K^f  ■ 


(2.24) 
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For  an  assumed  continuous  distribution  of  drops  sizes  with  diameter  density 
N(D)  Equation  (2.23)  becomes 

^max 

Z  =  I  D^N{D)dD.  (2.25) 

0 

Solving  the  integral  form  for  an  assumed  Marshall-Palmer  Drop  Size  Distribution 
(DSD)  Equation  (2.15),  yields 

00  ^  I  JLT 

Z={D^N(D)dD  = - (2.26) 

i  ^  ^  (4.ii?--')' 

Drops  do  not  exist  in  sizes  ranging  up  to  infinite  diameter,  but  are  likely  to  be 
limited  by  a  maximum  diameter.  Doviak  and  Zrnic  state  that  Equation  (2.26) 
overestimates  Z  and  suggest  a  truncated  DSD  with  D  =  .  In  this  case. 


Z  = 


(4.1i? 


-r{7,a) 


(2.27) 


where  a  =  AD^^^^and  y{l,a)  is  the  incomplete  Gamma  function.  Equations  (2.26)  and 


(2.27)  provide  an  analytical  means  for  computing  Z.” 


F.  RADAR  PRINCIPLES  AND  WEATHER  RADAR  EQUATION 
CORRECTION 


The  basic  results  for  phased  array  weather  radar  are  presented  from  Knorr’s 
paper.  This  paper  addresses  issues  unique  to  frequency  hopping,  phased  array  weather 
radars  and  quantifies  the  need  for  a  target  of  known  radar  cross  section  (RCS). 

“The  first  term  in  the  equation  below  indicates  those  radar  system  parameters 
which  may  not  be  precisely  known. 

~PP^Gl 

4 

where 

Xo  =  reference^  wavelength  (m) 


cos^ 


<J 


=  P.. 


(2.28) 


^  The  choice  of  reference  wavelength  within  the  operating  band  of  the  radar  is  arbitrary. 
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Ls  =  radar  system  losses  (ratio) 

La  =  atmospherie  propagation  loss  (ratio) 

Pout  ^  output  power  (W) 

R  =  range  to  target  (m) 

Pt  =  transmit  power 
Grx  =  receiver  gain 

Go  =  gain  at  the  reference  frequency, /o 
L  =  operating  wavelength  (m) 

6  =  angle  of  the  antenna  beam  with  respect  to  boresight  or  the  normal  to  the  array  face 
a  =  target  radar  cross  section  (RCS)  (m^) 

The  objective  of  calibration  is  to  determine  the  composite  of  these  factors  with  a 
single  measurement.  To  simplify  this  process,  several  conditions  can  be  selected.  These 
include:  calibration  on  boresight  with  the  antennas  of  the  active  calibrator  system  (ACS) 
and  the  radar  (0  =  0),  calibration  at  the  reference  RF  center  frequency  (^  =  ^o),  and 
calibration  at  a  range  that  is  short  enough  that  atmospheric  losses  can  be  ignored  (La  =  1). 
Then,  the  equation  simplifies  and  rearranges  to: 

~P,G^Gl 

4 

With  this  relationship,  a  measure  of  the  radar  receiver’s  output  power  and  range 
for  a  reference  target  of  known  cross-section  (the  ACS’s  equivalent  reflectivity)  will 
determine  the  radar  system  constant  (left  hand  side  of  the  equation).” 


(4;r) 


'  R^P.. 


(2.29) 


For  frequency  agile,  phased  array  weather  radars,  the  equation  for  reflectivity 
takes  the  form 


10241n2Zo 

1 

4 

L  P 

out 

7r^{cT)\K^A^ 

_P,G^Gl 

_(///o)"cos^_ 

where 

Z  =  reflectivity  (mm  /m  ) 

Ao  =  reference  wavelength  (m) 

Ls  =  radar  system  losses  (ratio) 

La  =  atmospheric  propagation  loss  (ratio) 

P  =  average  output  power  (W) 

(JUl 

R  =  range  to  target  (m) 
c  =  speed  of  light  in  a  vacuum  =  3x10^  m/s 
X  =  pulsewidth  (s) 


(2.30) 
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r 


r 


.93  (H^O) 

Sr  =  hydrometeor  relative  permittivity 

^„„and  (3„„are  the  broadside  3  dB  prineipal  plane  beamwidths  at  the  referenee 

£>\J  £>[} 

frequeney 

Pt  =  transmit  power 

Grx  =  ree  elver  gain 

Go  =  gain  at  the  referenee  frequeney, /o 
/=  operating  frequeney  (Hz) 
fo  =  referenee’^  frequeney  (Hz) 

6  =  angle  of  the  antenna  beam  with  respect  to  boresight  or  the  normal  to  the  array  face 

Note  that  the  units  of  Z  are  given  correctly  (mm  /m  )  by  the  first  two  terms  of 
Equation  (2.30)  if  range  is  in  millimeters,  Zg  in  the  numerator  of  the  first  term  is  in 
millimeters  and  the  quantities  in  the  denominator  of  the  first  term  are  in  meters. 

This  is  the  corrected  form  of  the  classical  Probert-Jones  classical  weather  radar 
equation  as  derived  by  Dr.  Jeffrey  Knorr  in  his  paper,  “Weather  Radar  Equation 
Correction  for  Erequency  Agile  and  Phased  Array  Radars.” 

The  third  term  contains  radar  system  parameters  which  may  not  be  precisely 
known.  The  objective  of  calibration  is  to  determine  the  composite  effect  of  these  factors. 
The  composite  of  these  factors  is  also  known  as  the  radar  system  calibration  constant.  In 
Dr.  Knorr’s  paper  it  has  been  shown  that,  for  the  case  where  the  reference  target  is  on 
boresight,  6  =  0,  the  measurement  is  made  at  the  reference  frequency,  2  =  Xq,  and  the 
range  is  short  {La=  1) 

4 

Prom  Equation  (2.31),  it  is  clear  that  a  measurement  of  receiver  output  power  and 
range  for  a  reference  target  of  known  cross-section  {6  =  0,  X  =  Xq,  La  =  1)  permits  the 
radar  system  constant  on  the  left  side  to  be  determined  and  the  radar  calibrated  for 
reflectivity  estimation. 


(4^)' ^'4. 

Xlcj 


(2.31) 


^  The  choice  of  reference  frequency  within  the  operating  band  of  the  radar  is  arbitrary. 


G. 


ESTIMATION  OF  VELOCITY 


According  to  Schroder,  “Pulsed  radar  ean,  by  measuring  the  phase  change 
between  eonseeutive  pulses,  estimate  radial  veloeity.  The  effeet,  ealled  Doppler,  is  the 
frequeney  ehange  that  oeeurs  due  to  the  relative  movement  of  a  target  with  respeet  to  the 
radar.  Considering  a  two-way  propagation  path  renders  a  total  phase  ehange  of 

2  R 

(f>  =  2n —  (2.32) 

X 


In  terms  of  angular  frequeney  the  result  is 

_d(l)  _  An  dR  _  Anv^ 
dt  X  dt  X 


=  ^^fd 


(2.33) 


where  is  the  radial  veloeity  and  is  the  Doppler  frequeney. 


From  a  pulse  pair,  using  a  eoherent  deteetor,  the  phase  ehange  ean  be  extraeted 


using 

V  =  I  +  jO  =  A 

n  n  J  n 

V  V*=A  ^ 

^  n+V  n 

leading  to  the  single  pulse  pair  veloeity  estimate 


4;r 


(2.34) 


(2.35) 


where  T  is  time  between  pulses  and  /  is  pulse  repetition  frequeney.’ 


Also  from  Sehroder,  eoherenee  between  pulse  pairs  is  required  to  measure  phase 
ehange  over  time  for  a  volumetrie  target  to  minimize  phase  change  from  effeets  other 
than  radial  movement.  The  varianee  of  a  veloeity  estimate  ean  be  improved  simply  by 
averaging  more  samples.  Building  on  the  single  pulse  pair  veloeity  estimate  Equation 
(2.35),  the  multiple  pulse  pair  estimator  beeomes: 


$ 

V  =  - 


1  N-\ 

— 

iv-is 


1 


N-\ 


N-\  An  tt 


(2.36) 
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H.  ESTIMATION  OF  VELOCITY  SPREAD 


Turbulence,  wind  shear,  and  differences  in  fall  speed  due  to  drop  size  are  but  a 
few  of  the  many  directly  influencing  factors  on  the  spread  in  radial  velocity  of 
precipitation.  There  are  even  indirect  factors  (motion  of  the  radar  antenna)  which  impart 
an  apparent  spread  in  velocity.  The  net  effect  of  these  factors  is  the  width  of  the  spread 
in  the  velocity  spectrum  and  the  time  for  de-correlation  between  samples. 


Schroder  explains,  “An  estimate  of  the  velocity  spread  can  be  computed  using  the 
velocity  samples  obtained  from  a  pulse  pair  comparison  between  single  pulse  pair 
velocity  estimates  Equation  (2.35).  The  estimator  for  sample  variance  is  given  by 


m 


i=\ 


S' 


N 


and  the  estimator  for  standard  deviation  is” 

. 

I,  CALIBRATION  REQUIREMENTS 


(2.37) 


(2.38) 


Factors  exist  both  internally  (hardware  and  software)  and  externally 
(environmental  and  remote  sensing  effects)  that  influence  a  radar’s  measurements  apart 
from  the  target  return  that  is  trying  to  be  quantified.  The  calibration  process  is  used  to 
reduce  the  influence  of  these  factors  on  the  composite  measurement.  Internal  calibrators 
minimize  system  variations  by  separating  the  instrument  from  variable  external 
influences  (such  as  through  the  use  of  known  reference  signals  generated  internally  to  test 
internalized  system  components).  External  calibrators,  however,  operate  completely 
outside  (and  independent  of)  the  instrument  under  test.  By  operating  completely  outside, 
external  calibrators  incorporate  those  elements  in  the  calibration  process  that  describes 
external  influences  that  would  effect  the  instrument’s  measurements,  but  are  not 
considered  with  an  internal  calibration  method.  By  operating  independently,  external 
calibrators  are  able  to  calibrate  for  the  same  factors  targeted  by  internal  calibrators 
without  the  calibrator  being  influenced  in  its  calibration  by  the  instrument  under  test. 
With  these  two  additional  traits,  external  calibrators  can  verify  whole-system  calibration. 
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The  result  is  a  true  “end-to-end”  ealibration  of  the  speeifie  radar  eonfiguration  for  the 
speeifle  external  measurements  that  are  about  to  be  aequired. 

Beeause  these  internally  and  externally  influeneing  faetors  are  so  many  and  so 
varied  they  eannot  be  individually  quantified  in  theoretieally  expressions.  Instead,  these 
eombined  effeet  faetors  are  often  denoted  by,  in  terms  of  Pratte,  “the  so-ealled  “radar 
eonstant”  C.  In  any  given  implementation,  the  radar  eonstant  C  ean  eontain  slowly 
varying  quantities,  sueh  as  reeeiver  eonversion  gain,  system  and  site  losses  (to  inelude 
propagation  path  loss  eaused  by  the  atmosphere),  and  “hidden”  assumptions.  In  WSR- 
88D  software  jargon  the  radar  eonstant  derived  from  the  engineering  ealibration  is  ealled 
“SYSCAL”  (Sirmans  1992),  and  temporal  variations  obtained  by  the  internal,  eleetronie 
ealibration-maintenanee  funetion  software  is  ealled  “delta  SYSCAL”  or  “CAL#”.” 

J.  CALIBRATION  METHODS  AND  SELECTED  METHOD 

A  summary  of  alternative  ealibration  methods  was  obtained  from  Pratte.  A  list  of 
these  methods  follows: 

1 .  Traeeability 

2.  Solar  and  Celestial  Radiative  Flux  Transfer  Methods 

3.  Radar  Calibration  Targets 

4.  Lunar  Refieetivity 

5.  Transponder  Methods 

6.  Elevated  Ground-range  Methods 

7.  Near-Field  Measurements 

8.  Refieetivity  Comparison  Methods 

9.  Rain  gauge  Comparison  Methods 

10.  Dual  Polarization  Consisteney  Methods 

Compared  to  other  ealibration  methods,  use  of  an  external  aetive  ealibration 
system  is  the  most  reliable,  independent  and  thorough  means  of  ealibration.  For  this 
reason,  this  type  of  system  was  ehosen  for  implementation. 

The  approaeh  is  a  repeater  that  will  produee  a  known  RCS  via  loop  gain  eontrol 
and  add  Doppler  return  information  via  phase  modulation  with  a  eontrolled  sawtooth. 
Proper  eontrol  of  the  sawtooth  frequeney  introduees  Doppler  spread,  an  important 
measurement  quantity  in  the  weather  radar.  This  approaeh  assumes  a  smooth  profile  in 
Doppler  return  as  a  funetion  of  the  varying  sawtooth.  A  ealibration  sphere  is  seleeted  to 
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enable  that  smooth  response  profile  (see  Figure  6).  Although  a  constant  amplitude  return 
would  seemingly  suffice,  use  of  a  programmable,  constant  phase  attenuator  would  also 
allow  software  control  of  the  amplitude  statistics  of  the  return  to  the  radar.  Thus,  with 
full  control  of  both  amplitude  and  phase,  one  could  replicate  the  return  for  each  of  the 
Swerling  target  types,  although  it  is  only  the  precipitation  return  that  is  really  of  primary 
interest  in  our  application. 

K,  SUMMARY 

This  chapter  has  lead  up  to  the  formulations  of  the  radar  equation  for  a  point 
target  and  for  a  volume  target  which  show  that  the  radar  calibration  constant  can  be 
determined  if,  having  a  target  of  known  RCS,  the  echo  power  from  that  target  is 
measured.  The  objective  of  this  thesis’s  research  is  to  determine  if  an  active  calibrator 
can  be  built  which  will  produce  the  desired  known  RCS  over  a  response  range  of  interest 
for  weather  returns.  Introducing  a  known  Doppler  shift  and  known  RMS  frequency 
spread  allows  one  to  compare  the  radar  estimates  of  those  parameters  with  the  values  set 
for  the  calibrator.  The  Doppler  shift  also  addresses  problems  associated  with  use  of  the 
calibrator  in  an  environment  where  there  is  fixed  clutter  return. 
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III.  DESIGN,  IMPLEMENTATION  AND  OPERATION 


A,  DESIGN 

This  study  implemented  an  existing  design  (see  Figure  2)  of  an  active  calibrator 
system  (ACS)  for  the  MWR-05XP.  It  did  not  involve  the  creation  of  such  a  design  nor 
did  it  involve  the  investigation  of  alternatives  to  or  optimization  of  this  design. 


Figure  2.  Block  Diagram  of  the  ACS. 

See  Figure  3  for  a  layout  of  the  actual  components  used.  These  components 
comprise  the  internals  of  the  Active  Calibrator  System  or  ACS  that  was  implemented  in 
this  study. 
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The  Active  Calibrator  (internal) 
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Figure  3.  Diagram  of  Component  Layout  inside  the  ACS  Enclosure. 


B.  HARDWARE  IMPLEMENTATION 

1,  System  Description  Overview 

See  ACS  Pictures  Appendix. 

The  active  calibrator  system  (ACS)  is  comprised  of  two  Pasternack  sectoral  horn 
antennas  mounted  on  either  end  of  a  6  ft  long  l”x4”  piece  of  lumber.  To  this  piece  of 
lumber,  another  piece  of  lumber  of  the  same  dimensions  is  attached  to  create  an  “L” 
bracket.  Midway  between  the  two  antennas  is  centered  the  ACS  enclosure,  an  aluminum 
box  (17”xl3”x5”),  to  house  the  power  supplies,  control  computer,  and  RF  components 
(programmable  attenuator,  digital  phase  shifter,  etc.).  On  the  face  of  this  enclosure, 
parallel  to  the  antennas’  mainbeams,  is  attached  a  second  piece  of  lumber  (cut  to  the 
length  of  the  enclosure  face)  to  present  a  flush,  all  wood  box  exterior.  Along  this  exterior 
and  running  the  entire  length  of  the  “L”  bracket  between  the  two  antennas,  X-band  Radar 
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Absorbent  Material  (RAM)  is  attached.  Although  this  RAM  has  its  greatest  RF 
absorbency  at  9.4  GHz,  it  provides  sufficient  RF  absorbency  across  the  entire  RF  range  of 
interest,  9-10  GHz,  to  reduce  undesired  RF  returns  greater  in  power  than  the  transmitted 
power  of  the  ACS’s  transmit  antenna  output.  This  RAM  is  sectioned  into  three  pieces 
(two  24”xl0”x2”  and  one  19”x9”x2”).  Each  piece  is  attached  via  removable  Velcro  to 
aid  in  transport  and  storage.  The  antennas  are  mounted  and  oriented  to  transmit  and 
receive  horizontally  polarized  signals.  They  are  connected  to  the  enclosure  via  3  ft  RF 
cables.  All  of  the  ACS’s  components  apart  from  the  antennas  are  contained  inside  the 
enclosure. 

2,  Hardware  List 

External  Hardware  (Quantity/Description/Dimensions,  if  applicable) 

(1)  Aluminum  box(enclosure)  (17”xl3”x5”) 

(3)  Pieces  Emerson  &  Cuming  AN77  RAM  (two,  left  and  right  side,  24”xl0”x2”  and 
one,  center,  19”x9”x2”) 

(3)  Pieces  lumber  (l”x4”,  two  6  ft  long  and  one  19”  long) 

(2)  Pasternack  Sectoral  Horn  Antennas 

(2)  Pasternack  semi-flexible  RE  Cables  (3  ft) 

(1)  Ethernet  (RJ-45)  Cable 

(1)  Dell  Precision  M65  with  2.33  GHz  Core  2  Duo  Processor,  2GB  667  MHz  RAM,  80 
GB  HDD  Laptop  and  Dell  90W-AC  Power  Supply 
(1)  AC  Power  Cable 

Miscellaneous  Mounting  Hardware  (nuts,  screws,  washers,  bolts,  etc.) 

Internal  Hardware  (Quantity/Description) 

(1)  Laptop  Power  Supply  (24V  CUI  EA1050B-240) 

(1)  Condor  DC  Power  Supply  (GLC-50A) 

(2)  25-pin-to-15-pin  Control/Power  Cables 
(1)  Voltage  Regulator  Circuit 

(1)  National  Instruments  Control  Computer  (CompactRIO  Model  Line  to  include;  (1) 
real-time  controller,  (1)  reconfigurable  embedded  chassis,  and  (2)  I/O  modules) 

3,  Voltage  Regulator  Circuit 

A  voltage  regulator  circuit  was  required  to  convert  the  DC  power  supplied  by  the 
Condor  power  supply  (-15VDC,  +15VDC,  and  +5VDC)  to  the  additional  voltages  (- 
5VDC  and  -I-12VDC)  required  by  the  RE  components. 
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Voltage  Regulator  Circuit  Hardware  List  (Quantity/Description) 
(2)  0.01  uF  Capacitors 
(1)33  ohm  Resistor 

(1)  7905C  (provides  -5VDC  from  -15VDC  input) 

(1)  7812CT  (provides  +12VDC  from  +15VDC  input) 


Figure  4.  Voltage  Regulator  Circuit  Diagram. 

4,  RF  Components 

See  Appendix  on  Digital  Phase  Shifter. 

See  Appendix  on  Programmable  Attenuator. 

See  Appendix  on  Antennas.  Appendix  describes  antenna  specifications. 

See  Appendix  on  Antenna  Testing.  Appendix  describes  antenna  patterns,  gain 
and  beamwidth  calculations,  and  isolation. 

5,  Control  Hardware 

Control  of  the  ACS  is  achieved  through  the  use  of  the  National  Instruments  (NI) 
Compact  Reconfigurable  Input/Output  (CompactRIO)  product  line.  From  the 
ni.com/compactrio  website,  “The  National  Instruments  CompactRIO  programmable 
automation  controller  (PAC)  is  a  low-cost  reconfigurable  control  and  acquisition  system 
designed  for  applications  that  require  high  performance  and  reliability.  The  system 
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combines  an  open  embedded  architecture  with  small  size,  extreme  ruggedness,  and  hot- 
swappable  industrial  I/O  modules.  NI  CompaetRIO  is  powered  by  reconfigurable  I/O 
(RIO)  FPGA  technology.”  Again  from  the  ni.eom/pac  website,  “PACs  combine 
programmable  logie  controller  (PLC)  ruggedness  with  PC  functionality  under  an  open, 
flexible  software  arehitecture.” 

The  CompaetRIO  products  used  by  the  ACS  include:  a  real-time  controller,  a 
reconfigurable  embedded  chassis,  and  two  I/O  modules  for  digital  input.  Software 
development  was  done  on  a  Windows  PC  laptop  with  NTs  Developer  Suite  with  Real¬ 
time  deployment  and  FPGA  deployment  options  installed.  Through  this  developer  suite 
and  these  options,  all  software  development  was  done  in  the  NI  Lab  VIEW™  visual 
programming  language.  For  items  not  run  on  either  the  laptop’s  Windows  0/S  or  the 
real-time  controller’s  real-time  0/S,  LabVIEW  was  auto-complied  into  the  VHSIC 
Hardware  Description  Eanguage  (VHDE)  and  run  on  the  reconfigurable  embedded 
chassis.  This  process  was  seamless  to  the  developer  and  no  specific  knowledge  of  VHDE 
was  required.  These  items  included  all  of  the  software  control  required  to  manipulate  the 
phase  shifter  and  programmable  attenuator  on  a  time  scale  at  or  near  its  fastest  possible 
switehing  speed.  A  normal  0/S  or  real-time  0/S  is  incapable  of  this  speed  of  control 
requiring  the  use  of  Eield  Programmable  Gate  Arrays  (EPGAs). 

See  Appendix  on  CompaetRIO. 

6,  Hardware  Limitations 

There  are  a  few  hardware  limitations  associated  with  the  ACS.  These  limitations 
affect  the  maximum  possible  Doppler  shift  and  therefore  radial  velocity  and  velocity 
spread  that  ean  be  practically  generated.  Eimited  factors  are  dictated  by  the  switehing 
speed  of  the  programmable  attenuator  (300  ns)  and  phase  shifter  (500  ns),  the  CRIO 
EPGA  cloek  (25ns),  and  the  response  time  of  the  digital  I/O  module  (100  ns).  Obviously, 
the  switching  speeds  of  the  programmable  attenuator  and  phase  shifter  are  the  true 
limiting  factors  as  the  CRIO  components  are  more  than  eapable  of  switching  states  at  or 
faster  than  the  fastest  switching  speeds  of  both  RE  components. 
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See  the  example  calculation  of  the  effect  of  the  phase  shifter’s  fastest  possible 
switching  speed. 

1 


2'x(500xl0 ‘') 


=  7812.5  Hz  Doppler  Frequency  Shift 


299792458 

7812.5  X - =  117.2  m/s  Radial  Velocity 

2x9990000000 

1 17.2  X  2.23693629  =  262.2  mph  Radial  Velocity 
This  limitation  actually  falls  below  what  is  physically  possible  (or  what 
has  been  observed)  and,  therefore,  could  potentially  limit  ACS  applicability  in  a  few 
applications.  The  fastest  wind  speed  ever  recorded  on  Earth  was  318  mph  in  one  of  the 
tornadoes  that  hit  the  suburbs  of  Oklahoma  City  on  3  May  1999.^  However,  the  fastest 
recorded  wind  speed  not  associated  with  a  tornado  was  measured  at  230  mph  in  1934.9 


ACS  FREQUENCY/DOPPLER  CALIBRATION  AND  SOFTWARE 
IMPLEMENTATION 


In  order  to  calibrate  the  weather  radar’s  ability  to  measure  radial  velocity  and 
velocity  spread,  the  ACS  must  be  able  to  produce  a  known  Doppler  shift  and  Doppler 
spread.  By  comparing  the  known  quantity  to  the  measured  quantity  and  adjusting  for  any 
discrepancy,  the  calibration  can  be  completed  and  the  weather  radar’s  accuracy  assured. 
This  section  describes  how  the  ACS  achieves  this. 


The  active  calibrator’s  loop  gain  determines  equivalent  reflectivity,  Zeq.^*’  Phase 
modulation  is  required  to  produce  Doppler  shift,  according  to: 


v(t)  =  Hcos[(y^t  +  6’(t)  =HRe 


(3.1) 


Following  the  derivations  in  Tan’s  thesis  [9]  and  the  paper  by  Cummings  [10],  for 
0(t)  periodic. 


^  National  Weather  Serviee  Foreeast  Offiee,  "Frequently  Asked  Questions  about  the  May  3,  1999 
Bridge  Creek/OKC  Area  Tornado,"  retrieved  2008-09-21  from 
http://www.srh.noaa.gOv/oun/storms/19990503/may3faqs.php 

9  University  Corporation  for  Atmospherie  Researeh  (UCAR),  "Windows  to  the  Universe/Wind," 
retrieved  2008-09-21  from  http://www.windows.uear.edu/tour/linkAearth/Atmosphere/wind.html 

See  seetion  on  “ACS  Amplitude/Refleetivity  Calibration”  (next  seetion)  for  a  diseussion  on  the 
determination  of  equivalent  refleetivity. 
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(3.2) 


(3.3) 


Figure  5.  Two  Cyeles  of  the  Periodie,  Linear,  Modulating  Sawtooth  Funetion. 


=  for 


- <t<+- 

2  2 


c=T 


sin(A6’/2-n;r) 


(3.4) 


(A6*/2-n;r) 

In  general,  the  phase  modulated  signal  has  a  general  signal  given  by 


(A6’/2-n;r) 


(3.5) 


If  the  amplitude  of  the  sawtooth  is  adjusted  to  produee  a  phase  change  that  is  a  multiple 
of  271  then 

A0  =  mlTT  for  (m  e  /) 

rifor77  =  m  .  (3.6) 

T  [O  for  77  m 

Therefore, 

v(z)  =  where  (A^  =  7?72;r) .  (3.7) 
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The  real  part  of  the  eomplex  exponential  expression  is  cosine  with  argument  indicating 
the  translation  (Doppler  shift)  from  frequency  to  frequency  +  mf^  .  In  other  words, 
each  spectral  line  is  separated  by  the  modulating  frequency  and  lines  exist  symmetrically 
above  and  below  the  carrier. 


By  cycling  through  360°  (2n)  of  phase  at  this  periodic  rate  a  known 
frequency/Doppler  shift  is  produced  against  which  the  weather  radar  can  be  calibrated. 
This  periodic  phase  cycling  (and  thus  frequency  translation)  is  implemented  in  the  ACS 
using  a  digital  phase  shifter.  However,  there  are  imperfections  associated  with  this 
method.  Imperfect  implementation  results  in  spurious  sidebands  due  to  finite  sawtooth 
flyback  time,  incorrect  sawtooth  amplitude,  the  digital  staircase  approximation  to  linear 
phase  shift,  phase  error,  and  finite  phase  shifter  switching  time.  [10] 


Again,  following  the  derivations  in  the  thesis  by  Tan  [8],  for  a  digital  phase  shifter 


using  N=2^  steps. 


/  \  (“!)”*  sin(;r/A)  .  . 


(3.8) 


+  l  I 

where  K  =  mN  + 1  and  m  is  an  integer.  The  frequencies  of  the  non-zero  spectral 
components  are  given  by  /  =  /^  or  /  =  X  +{tnN  +  \)  .  The  amplitude  of  the  K* 


sideband  relative  to  the  input  carrier  is  given  by 


The  first  sideband  in  the  spectrum  of  the  translated  signal  is 


(') 


sin(;r/ a) 
(;r/  a) 


Aco?,[o)^+Ko)^)t . 


The  translation  loss  is  defined  as 


(3.9) 


(3.10) 


=-201og 


sin(;z-/  A) 
(;r/ a) 


=  -20  log 


sin(^;r/2^  j 
(;r/2") 


(3.11) 


This  loss  rapidly  approaches  0  dB  as  the  number  of  bits  of  phase  shift  increase.  If  B=4 
bits,  for  example,  TLdB=0.06  dB. 
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The  amplitudes  of  the  undesired  sidebands  are  also  of  interest.  The  undesired 
sideband  with  the  highest  amplitude  is  K  =  -N  + 1 .  Relative  to  the  desired  sideband, 
K=l,  the  suppression  ratio  is 

=201og(iV-l)  =  201og(2''-l).  (3.12) 

The  suppression  ratio  also  improves  rapidly  as  the  number  of  bits  is  inereased.  If  N=  4, 
for  example,  SRdB=  23.5  dB. 

See  Appendix  I  on  Lab  VIEW  for  the  Digital  Phase  Shifter  which  gives  further 
explanation  of  how  the  phase  shifter  is  controlled  by  the  CompactRIO  to  produce  the  In 
periodic  phase  cycling. 

In  conclusion,  one  of  the  primary  advantages  of  the  ACS  is  its  ability  to  introduce 
Doppler  shift.  This  permits  the  calibrator  signal  to  be  extracted  from  mainbeam  ground 
clutter  at  elevation  angles  near  zero  and  from  sidelobe  clutter  at  higher  elevation  angles. 
In  most  cases,  convenience  will  dictate  that  the  ACS  be  operated  near  the  ground. 

D,  ADDING  VELOCITY  SPREAD 

Velocity  spread  can  be  added  to  the  Doppler  shifted  signal  by  varying  the 
frequency  of  the  sawtooth  waveform  used  in  the  phase  modulator.  The  frequency  should 
be  varied  in  accordance  with  a  Gaussian  random  variable  that  has  been  low  pass  filtered 
to  produce  a  modulating  signal  with  the  desired  correlation  time.  This  functionality  was 
not  fully  implemented  in  the  ACS  at  the  conclusion  of  this  study.  See  Chapter  IV  on 
Future  Work. 

E,  ACS  AMPLITUDE/REFLECTIVITY  CALIBRATION 

In  order  to  calibrate  the  weather  radar’s  ability  to  measure  reflectivity,  the  ACS 
must  be  able  to  reliably  simulate  a  known  reflectivity.  This  is  accomplished  by  receiving 
the  incident  radar  signal  and  retransmitting  a  signal  of  known  power  relative  to  the 
incident  signal.  From  a  measurement  of  the  echo  signal  power  received  by  the  radar  and 
knowledge  of  the  ACS  loop  gain,  the  calibration  can  be  completed  and  the  weather 
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radar’s  absolute  accuracy  in  amplitude  assured.  This  section  describes  how  the  ACS  is 
calibrated  so  it  can  in  turn  be  used  to  calibrate  a  weather  radar. 

1.  Measurement  Premise 


The  ACS  needs  to  produce  a  power  density  equivalent  to  that  scattered  by  a  target 
of  known  RCS.  This  power  will  be  the  effective  radiated  power  (or  ERP)  of  the  ACS. 
The  ERP  of  the  ACS  will  be  a  function  of  (depend  on)  the  ERP  of  the  radar.  This  is 
because  the  ACS  does  not  produce  an  ERP  from  an  independent/internal  power  source, 
but  rather  as  acts  a  repeater  (and  modifier)  of  the  power  received  from  the  radar. 


The  power  received  from  the  radar  will  depend  on  the  ERP  of  the  radar  and  the 
losses  involved  between  the  ACS  and  the  radar.  At  this  radar  frequency,  the  losses 
between  the  ACS  and  the  radar  will  primarily  depend  on  the  distance  between  the  ACS 
and  the  radar  (spreading  loss  or  free  space  loss  only,  unless  a  ground  reflection  is  present 
as  will  be  the  case  for  elevation  angles  near  zero)  and  atmospheric  attenuation 
(propagation  loss),  but  also  losses  internal  to  the  radar  (system  losses  which  cannot  be 
precisely  known).  Those  losses  which  depend  on  distance  will  affect  the  reflectivity 
observed  by  the  radar  due  to  the  calibrator  twice  -  first,  as  the  ERP  of  the  radar  is 
affected  on  its  way  to  becoming  the  received  power  of  the  calibrator  and,  second,  as  the 
ERP  of  the  calibrator  is  affected  on  its  return  to  becoming  the  received  power  of  the 


radar. 


While  the  spreading  loss  is  well  known  {L  = 


2 

)  and  can  be  accurately 


calculated,  the  propagation  loss  is  less  well  known  (as  it  depends  on  many  time  variant 
factors)  and  is  more  difficult  to  accurately  calculate.  In  order  to  isolate  this  uncertainty  at 
any  given  time,  the  equivalent  reflectivity  of  the  ACS  can  be  compared  to  a  known 
reflectivity  at  equal  distances  and  with  the  propagation  losses  minimized  as  much  as 
possible  (as  in  an  anechoic  chamber).  By  comparing  these  two  values,  the  equivalent 
reflectivity  of  the  ACS  can  be  measured.  Under  these  settings  the  losses  internal  to  the 
radar  can  be  isolated  (the  system  losses  which  cannot  be  precisely  known  will  be 
calibrated  out). 
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An  equation  for  the  ERP  of  the  ACS  and  thereby  the  equivalent  refleetivity  of  the 
ACS  ean  be  defined  as  a  funetion  of  the  reeeived  power  of  the  ACS  (from  the  radar)  and 
the  loop  gain  (system  gain  and  losses)  of  the  ACS.  These  inelude  the  gain  of  both 
antennas  (transmit  and  reeeive),  eable  losses,  phase  shifter  losses,  programmable 
attenuator  losses,  isolator  losses,  amplifier  gain,  and  other  RF  eomponents’  losses  (only 
the  amplifier  provides  gain).  On  the  whole,  the  ACS’s  ERP  and  therefore  its  equivalent 
reflectivity  can  be  controlled  (varied)  by  the  programmable  attenuator  (variable 
gain/loss).  However,  this  is  only  effective  as  long  as  the  net  gain/loss  of  all  of  the  ACS’s 
RE  components  (fixed  gain/loss)  in  total  exceed  the  reflectivity  of  the  ACS  due  to  its 
exposed  surface  area. 

The  equivalent  reflectivity  due  to  the  ACS’s  transmit  antenna  is  in  direct 
competition  with  the  actual  reflectivity  due  to  the  ACS’s  surface  area.  Both  depend  on 
the  radar  ERP,  but  the  equivalent  reflectivity  further  depends  on  the  ACS’s  loop  gain 
acting  to  re-radiate  the  radar’s  power  while  the  actual  reflectivity  depends  on  the  exposed 
surface  area  of  the  ACS  (to  include  all  sides  facing  the  radar)  acting  to  reflect  the  radar’s 
power.  If  the  controlled  power  does  not  exceed  this  value,  variable  amplitude  calibration 
(as  will  be  needed  for  amplitude  statistics)  will  not  be  possible.  Instead,  only  a  single 
reflectivity  (the  actual  reflectivity)  of  the  ACS  will  ever  be  observed  by  the  radar.  This  is 
the  same  principle  used  in  passive  calibration  systems  (such  as  the  balloon  and  sphere) 
where  a  controlled  Doppler  calibration  is  not  possible.  Radar  Absorbent  Material  (RAM) 
was  used  to  reduce  the  actual  reflectivity  to  levels  below  the  equivalent  reflectivity.  See 
Chapter  III,  Section  B,  Subsection  I  for  further  details  on  RAM  use. 

2.  Measurement  Method 

What  follows  is  the  method  used  to  calibrate  out  these  unknown  system  losses  for 
accurate  reflectivity  measurements  (estimates)  by  the  weather  radar  based  on  the 
equivalent  RCS  of  the  ACS.  See  Chapter  II,  Section  F,  on  “Radar  Principles  and 
Weather  Radar  Equation  Correction’’  for  background  and  supporting  theory  directly 
applicable  to  this  process  and  these  calculations.  The  following  development  shows  how 
the  equivalent  RCS  of  the  ACS  was  experimentally  determined  through  anechoic 
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chamber  measurements  and  subsequent  ealculations  to  satisfy  a  most  fundamental  radar 
calibration  requirement  -  namely,  that  the  ACS’s  equivalent  RCS  must  be  known  [5]. 

As  mentioned  before,  the  value  of  the  equivalent  reflectivity  (or  more  accurately 
equivalent  radar  cross  section  or  RCS)  of  the  ACS  can  be  measured  and  calculated  in 
eomparison  to  a  known  reflectivity  (or  RCS)  of  another  target.  In  this  study,  an 
aluminum  sphere  was  used  and  measurements  made  inside  an  aneehoie  ehamber  at  a 
fixed,  known  distance  to  serve  as  the  known  refieetivity.  A  12”  diameter  aluminum 
sphere  was  chosen  out  of  three  diameters  available  based  on  the  initial  caleulations  and 
estimates  performed  below. 

a.  Initial  Calculations 

Prior  to  making  measurements  in  the  aneehoie  chamber  theoretical 
calculations  were  made  using  approximately  known  or  individually  measured  values  of 
the  required  quantities.  The  signal  source  in  the  aneehoie  chamber  was  a  vector  network 
analyzer  (VNA).  The  VNA  replaces  the  radar  as  the  signal  source  when  calibrating  the 
ACS  in  the  aneehoie  chamber. 

Below  is  a  list  of  constants,  either  measured  or  known,  to  include; 

R,  the  measured  one-way  distanee  from  the  transmit  and  receive  antennas  to  the  ACS 
inside  the  aneehoie  chamber,  (measured) 

f,  the  radio  eenter  frequeney.  (known,  seleeted  on  vector  network  analyzer  or  VNA) 
c,  the  speed  of  light,  (known) 

X,  the  wavelength,  (calculated) 

Pt,  the  VNA  transmit  power  a.k.a.  the  power  input  to  the  aneehoie  transmit  antenna  from 
the  VNA  to  include  all  losses  from  the  output  of  the  VNA  to  the  antenna  (cabling,  etc.), 
(known,  selected  on  VNA) 

Gt=Gr,  the  gain  of  the  VNA  antenna  -  one  in  the  same  -  a.k.a.  the  aneehoie  transmit  and 
receive  antennas  -  both  the  same,  (measured.  See  Appendiees  on  Antennas.) 

EIRP,  the  radar’s  Effective  Isotropic  Radiated  Power.  The  product  of  the  transmit  power 
and  antenna  gain,  (ealeulated) 

Gc,  the  calibrator’s  antenna  gain.  Both  antennas  on  the  calibrator  are  the  same, 
(measured.  See  Appendices  E  and  E  on  Antennas.) 

Gamp,  the  calibrator’s  loop  gain  to  include  all  system  losses  (cabling,  phase  shifter  while 
not  phase  shifting,  conneetors,  etc.)  -  i.e.  net  gain.  The  aetual  gain  of  the  amplifier  used 
in  the  calibrator  is  +35  dB  with  -20  dB  of  system  losses  for  +15  dB  net.  (ealeulated  via 
measured  values;  amplifier  gain  and  system  losses) 

FSL,  the  Tree  Space  Eoss  or  Spreading  Eoss  associated  with  the  distance  R.  (calculated) 
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Cl,  the  calibrator’s  input  signal  power  as  measured  into  the  system  at  the  receive 
antenna’s  output  port,  (calculated) 

RI,  the  VNA  input  signal  power  as  measured  into  the  system  at  the  receive  antenna’s 
output  port,  (calculated) 

/  =  9.50xl0®  Hz 
c  =  3.00x10*  m/s 

;L  =  — =  3.16x10^"  m 

/ 

=25  dBm  =316.2278  mW 
G,  =G^  =16  dBm  =39.8107  ratio 
EIRP  =  P,G,  =  41  dBm  =  12589.2541  mW 
=  13  dBm  =  19.9526  ratio 
^ amp  =15  dBm  =  31 .6228  ratio 
7?  =  5.7912  m  (19  feet) 


FSL  = 


X 


-\2 


AnR 


=  1.8830x10^’ ratio  =-67.2516  dB 


Cl  =  EIRP X  FSL  X  G  =  4.7298  x  1 0^"  mW  =  -1 3 .25 1 6  dBm 


RI  =  Cl  X  G„„.„  X  G  X  FSL  x  G  =  2.2371  x  10  "  mW  =  -36.2032  dBm 


amp 


Constants 

Quantity: 

Range 

Frequency 

Speed  of  Light 

Wavelength 

Radar  Transmit 

Radar  Antenna 

Radar  EIRP 

Calibrator 

Calibrator 

Power 

Gain 

Antenna  Gain 

Amplifier  Gain 

Symbol: 

R 

f 

c 

A=c/f 

p, 

G,=G, 

EIRP=P,*Gt 

Gc 

^amp 

Units: 

ft 

Hz 

m/s 

m 

dBm 

dB 

dBm 

dB 

dB 

Value: 

19 

9.50E+09 

3.00E+08 

3.16E-02 

25 

16 

41 

13 

15 

Alt.  Units: 

m 

mW 

ratio 

mW 

ratio 

ratio 

Alt.  Value: 

5.79 

316.23 

39.81 

12589.25 

19.95 

31.62 

Table  2.  Summary  of  Constants  Used  in  Estimated  Equivalent  RCS  Calculations. 


Calculations 

Quantity: 

Symbol/Equation: 

Units: 

Vaiue: 

Free  Space  Loss 
FSL=[A/(4*tt*R)]^ 
ratio  dB 

1.88E-07  -67.25 

Caiibrator  Input 
CI=EIRP*FSL*Gc 
mW  dBm 

4.73E-02  -13.25 

Radar  Input 
RI=CI*G3,^p*Ge*FSL*G, 
mW  dBm 

2.24E-04  -36.50 

Table  3.  Summary  of  Estimated  Equivalent  RCS  Calculations. 


-36dBm  is  approximately  the  value  of  VNA  input  signal  power  expected 
when  measuring  the  ACS  inside  the  anechoic  chamber.  This  will  be  the  total  returned 
power  from  the  ACS  due  to  the  relayed  power  and  not  the  reflected  power.  The  reflected 
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power  is  much  smaller  (no  return  was  detected  with  the  ACS  in  the  anechoic  chamber 
and  RAM  attached)  and  has  been  minimized  by  the  use  of  X-band  Radar  Absorbent 
Material  (RAM).  Again,  see  Chapter  III,  Section  B,  Subsection  1  for  further  details  on 
RAM  use. 

The  power  received  from  the  ACS  must  be  compared  to  the  received 
power  from  a  target  of  a  known  RCS  to  calculate  an  equivalent  RCS  for  the  ACS.  The 
known  RCS  that  was  used  was  that  of  a  polished  aluminum  sphere.  Three  spheres  were 
available,  a  3”,  a  6”  and  a  12”  diameter  sphere.  The  expected  received  power  was 
calculated  for  all  three.  The  relation  for  the  monostatic  RCS  of  a  sphere  in  the  optical 
region  (frequency  region)  was  used.  The  power  had  to  be  great  enough  to  be  detected  by 
the  VNA,  but  not  so  great  as  to  cause  saturation.  The  best  suitable  candidate  from  the 
three  available  was  found  to  be  the  12”  diameter  sphere  because  of  its  larger  (but  not  too 
large,  i.e.  saturating)  return  amplitude.  See  example  calculation  below, 
a  =  6  in  =  0.1524  m 

Pa  =  — ^  =  30.3  ^  Optical  Region  (See  Figure  6.) 

A. 

=  1  for  Optical  Region  (See  Figure  6.) 
na 

G  =  na^  =7.2966x10^  sm  = -11.3688  dBsm 
Received  Power  =  1.6338x10^  mW  =  -47.8679  dBm 

Following  the  example,  for  a  3”  diameter  sphere  -59.909  dBm  and  for  a  6”  diameter 
sphere  -53.889  dBm. 
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Sphere  Calculations 

Quantity: 

Sphere 

Radius 

Radius 

OPTICAL 

OPTICAL 

RCS 

RCS 

Received 

Received 

(Diameter) 

Power 

Power 

Symbol: 

- 

a 

a 

|3a=2TTa/A 

a/rra^ 

a  =  Tia^ 

a  =  Tia^ 

Units: 

- 

in 

m 

- 

- 

m^ 

dBsm 

mW 

dBm 

Value: 

3" 

1.5 

3.81  E-02 

7.58 

1 

4.56E-03 

-23.41 

1.02E-06 

-59.91 

Value: 

6" 

3 

7.62E-02 

15.16 

1 

1.82  E-02 

-17.39 

4.08E-06 

-53.89 

Value: 

12" 

6 

1.52E-01 

30.32 

1 

7. 30  E-02 

-11.37 

1.63E-05 

-47.87 

Table  4.  Summary  for  Calculations  of  Radar  Cross  Section  of  a  Sphere. 


b.  Anechoic  Chamber  Measurements 

Measurements  were  made  in  an  anechoic  chamber  with  the  Vector 
Network  Analyzer  system.  Three  different  series  of  measurements  were  made.  One  with 
the  chamber  empty  (E),  one  with  the  sphere  in  the  chamber  (S),  and  one  with  the 
calibrator  in  the  chamber  (C).  Measurements  were  made  and  represented  both  in  the  time 
domain  (T)  and  frequency  domain  (F).ii  To  eliminate  any  possible  interference  from 
reflections  (echoes)  in  the  chamber  due  to  less  than  ideal  performance  from  the  anechoic 
chamber  (there  were  design  compromises  made  in  the  chamber  at  the  time  of 


1 1  The  basic  measurement  is  a  swept  frequency  (frequency  domain)  measurement.  The  time  domain 
response  is  then  obtained  via  TFT.  Time  domain  resolution  is  determined  by  frequency  domain  bandwidth. 
A  1  GFlz  BW  should  give  about  1  ft.  resolution.  See  reference  [10].  It  provides  further  explanation  of  this 
information. 
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construction  due  to  space  limitations)  time  based  gating  was  used  to  eliminate  undesired 
signals.  Measurements  were  made  both  with  this  time  based  gating  (G)  and  without  (N). 
In  the  time  domain,  data  sets  and  plots  were  often  zoomed  to  the  time/distanee  of  interest 
(when/where  the  calibrator’s  return  was  expected).  These  are  annotated  zoomed  (Z)  or 
not  (A,  for  all).  When  not  zoomed,  they  inelude  all  data  starting  from  time  zero.  Table  5 
summarizes  the  abbreviations,  just  discussed,  that  are  used  when  referring  to  the  various 
anechoic  chamber  measurement  data  sets.  In  all  cases,  a  total  90  ns  time  base  was  used 
(RCS90). 


S=sphere 

T=time 

domain 

G=gated 

A=all 

(0  see  to  zoomed) 

Gl=54.8-58.8  ns 

E=empty 

chamber 

E=frequency 

domain 

N=not  gated 

Z=zoomed 
(around  time  gate) 

G2=39.8-43.8  ns 

C=calibrator 

NOTE:  does  not  apply 
to  E,  always  9-10  GHz 

Table  5.  Naming  Convention  Summary. 


Figure  7  is  a  composite  of  all  data  sets  in  the  time  domain,  not  zoomed 
(all).  Figure  8  is  a  composite  of  all  data  sets  in  the  frequency  domain.  All  data  sets 
inelude  the  measured  data  for  the  three  targets:  an  Empty  Aneehoic  Chamber  (chamber), 
the  12”  Aluminum  Sphere  (sphere)  and  the  ACS  (calibrator).  Timing  gates  were  placed 
around  the  sphere’s  location  (42  ns)  and  calibrator’s  location  (57  ns).  These  timing  gates 
eliminated  undesirable  echoes.  One  such  undesirable  echo  was  at  <2  ns.  This  echo 
corresponds  to  the  coupling  between  transmit  and  receive  antennas  which  were  placed 
about  18”  apart.  There  are  seven  graphed  data  sets  on  each  plot:  (1  and  2)  the  sphere 
without  a  time  gate  and  with  a  time  gate  eentered  on  its  42  ns  return,  (3  and  4)  the 
calibrator  without  a  time  gate  and  with  a  time  gate  centered  on  its  57  ns  return,  (5  and  6) 
the  empty  chamber  without  a  time  gate  and  with  a  time  gate  for  the  calibrator,  then  (7) 
with  a  time  gate  for  the  sphere. 
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Figure  7.  Time  vs.  Amplitude  for  Empty  Chamber,  Sphere  and  ACS. 

(ESC_plotTA) 
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Empty  Chamber,  Sphere,  &  Calibrator  -  Frequency  Domain 


Not  Gated  -  S 
Gated  -  S 
Not  Gated  -  C 
Gated  -  C 
Not  Gated  -  E 
Gated  for  Calibrator  -  E 
Gated  for  Sphere  -  E 


Figure  8.  Frequency  vs.  Amplitude  for  Empty  Chamber,  Sphere  and  ACS. 

(ESC_plotF) 

Eigure  9  and  Eigure  10  are  data  sets  in  the  time  domain  of  the  sphere  and 
the  calibrator  only,  respectively.  Each  plot  is  of  the  respective  target  only  and  zoomed  to 
very  near  the  time  gate  size  of  each  target.  Eor  the  sphere,  the  time  gate  (G2)  was 
centered  at  41.8  ns  and  spanned  39.8  to  43.8  ns.  For  the  calibrator,  the  time  gate  (Gl) 
was  centered  at  56.8  ns  and  spanned  54.8  to  58.8  ns.  The  time  gates  correlate  directly  to 
target  return  distances.  As  mentioned  before,  the  distance  (range)  from  the  transmit  and 
receive  antennas  in  the  anechoic  chamber  to  the  target  pedestal  inside  the  chamber  is  19 
feet.  19  feet  is  38.6  ns  in  time.  The  distance  is  traversed  twice  (once  to  the  target  and 
then  from  the  target)  at  the  speed  of  light.  The  reason  the  time  gate  for  the  sphere  was 
observed  at  and  consequently  set  at  42  ns  is  the  sphere  was  suspended  from  the  ceiling  of 
the  chamber  above  and  slightly  behind  the  target  pedestal  inside  the  anechoic  chamber. 
This  added  approximately  an  additional  1 .6  feet  of  distance  and  3  ns  of  time  because  the 
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range  hypotenuse  of  the  distance  (height  geometry  is  involved)  rather  than  the  longest  leg 
alone.  The  reason  the  time  gate  for  the  calibrator  was  observed  at  and  consequently  set  at 
57  ns  is  different  than  for  the  sphere.  This  time  gate  was  approximately  8  ns  greater 
because  of  the  time  delay  due  to  transmission  line  lengths  rather  than  target  distance.  The 
8  ns  is  equivalent  to  8  feet  (one-way  path),  slightly  greater  than  the  straight  line  distance 
between  the  transmit  and  receive  antennas  of  the  ACS  (as  expected).  See  the  example 
calculations  below  for  the  mathematical  relation  between  time  and  distance  (two-way 
path). 

c  =  3.00x10®  m/s 

i?  =  5.7912  m  (19  feet) 

9  /? 

—  =  38.6  ns 
c 

t  =  41.8  ns 

^  =  6.2700  m  (20.6  feet) 
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Sphere  -  Time  Domain 
(Aii=Zoomed) 


Not  Gated 
Gated 


Figure  9.  Time  vs.  Amplitude  for  Sphere,  Zoomed.  (S_plotTA=Z) 
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Calibrator  -  Time  Domain 
(AII=Zoomed) 


Figure  10.  Time  vs.  Amplitude  for  ACS,  Zoomed.  (C_plotTA=Z) 

The  expected  (calculated)  value  for  received  power  from  the  ACS  was  - 
36.5  dBm.  The  measured  value  for  received  power  from  the  ACS  was  -34.8  dBm.  The 
expected  (calculated)  value  for  received  power  from  the  12”  Sphere  was  -47.9  dBm.  The 
measured  value  for  received  power  from  the  12”  Sphere  was  -46.7  dBm.  All  expected 
(calculated)  values  very  nearly  equaled  the  measured  values. 

Figure  1 1  is  from  the  data  in  the  frequency  domain  of  both  the  sphere  and 
calibrator  previously  presented,  now  used  to  calculate  the  equivalent  RCS  of  the 
calibrator.  Measurements  were  made  from  9-10  GHz.  The  VNA  recorded  the  magnitude 
of  the  returned  amplitude  power  in  logarithmic  form  (log  MAG)  in  units  of  milliwatts 
(dBm).  All  measurements  are  in  the  frequency  domain  gated  for  time  (FG).  This 
includes  both  the  sphere  (SFG)  and  calibrator  (CFG).  As  can  be  seen,  the  measured 
values  were  very  close  to  the  expected  (calculated)  values.  The  received  power  for  the 
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12”  Aluminum  Sphere  was  calculated  to  be  -47.87  dBm  and  observed  to  be  near  -47 
dBm.  The  received  power  for  the  calibrator  was  calculated  to  be  -36.50  dBm  and 
observed  to  be  near  -35  dBm.  The  RCS  of  the  sphere  was  previously  calculated  to  be 
-11.37  dBsm  (S-RCS).  To  calculate  the  RCS  of  the  calibrator  (C-RCS)  based  on  the 
observed/measured  data,  the  CFG  -  SFG=A  (dB)  was  added  to  the  sphere’s  RCS.  Over 
the  entire  frequency  band  the  C-RCS  varied  from  0.6  to  3.2  dBsm.  Note  all  these 
measurements  in  the  chamber  included  the  use  of  the  RF  amplifier  within  the  calibrator 
system.  When  used  with  the  weather  radar,  signal  levels  are  great  enough  and  ranges 
short  enough  such  that  an  RF  amplifier  is  not  needed.  Furthermore,  calculations  are 
made  to  make  sure  the  signal  return  to  the  radar  is  in  the  proper  range  to  prevent  receiver 
saturation. 


Equivalent  RCS  of  Calibrator 
(with  in-line  amplifier) 


Frequency  (GHz) 


Figure  1 1 .  Frequency  vs.  Equivalent  RCS  of  the  ACS.  (RCS_plot) 

With  these  measured  values  the  equivalent  RCS  of  the  ACS  was  observed 
to  range  from  0.6  to  3.2  dBsm  across  the  9-10  GHz  RF  band.  On  average  (mean),  the 
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equivalent  RCS  of  the  ACS  was  1.4934  dBsm.  The  median  was  1.5101  dBsm.  The  high 
value  was  3.1742  dBsm  at  9  GHz  and  the  low  value  was  0.5862  dBsm  at  9.38  GHz. 

In  summary,  the  basic  radar  measurement  is  echo  signal  power  (for  each 
signal  sample).  From  the  samples,  average  power  is  estimated.  Reflectivity  follows 
directly  from  substitution  of  average  power  in  the  weather  radar  equation.  The  ACS 
provides  a  method  of  collectively  determining  the  imprecisely  known  radar  system 
parameters  appearing  in  the  weather  radar  equation.  Amplitude  statistics  from  the 
Weather  Signal  Simulator,  like  Doppler  statistics,  are  a  planned  addition  once  the 
programmable  attenuator  arrives  and  is  installed.  See  Chapter  IV  on  Future  Work  for 
other  planned  additions. 
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IV.  CONCLUSION  AND  FUTURE  WORK 

A,  CONCLUSION 

An  active  calibration  system  (ACS)  for  weather  radar  was  sueeessfully  designed 
and  implemented.  The  ACS  was  verified  to  be  a  funetional  weather  radar  aetive 
ealibration  system  under  laboratory  testing.  It  successfully  provided  the  means  (met  the 
requirements)  to  verify  all  three  of  the  signal  parameters  measured  (estimated)  by  a 
weather  radar.  These  include  the  zeroth  Doppler  moment  (reflectivity),  first  Doppler 
moment  (average  radial  velocity),  and  second  Doppler  moment  (velocity  spread). 

The  cross-seetion  (a,  or  equivalent  RCS)  of  the  ACS  and  therefore  the  means  to 
measure  and  ealeulate  the  radar  ealibration  eonstant  (for  refleetivity  aeeording  to 
Equation  (2.31))  was  sueeessfully  determined.  The  underlying  data  to  ealeulate  the 
cross-section  was  measured  and  based  on  subsequent  calculations  determined  to  vary 
from  0.6  to  3.2  dBsm  aeross  the  frequeney  band  of  interest. 

A  eonstant  Doppler  shift  (radial  veloeity)  ereated  by  the  ACS  was  observed  and 
verified  by  the  visual  inspeetion  of  several  independent  human  operators  using  a 
speetrum  analyzer.  The  ACS  sueeessfully  produeed  a  known,  operator  seleeted  radial 
velocity.  While  it  was  not  possible  to  observe  a  variable  Doppler  shift  (velocity  spread) 
in  exaetly  the  same  manner  because  of  time  scale  limitations,  the  method  to  produee  a 
variable  Doppler  shift  was  manually  verified  in  a  similar  manner.  The  ACS  produced  an 
observable  eonstant  Doppler  shift  when  provided  a  variable  Doppler  data  file  containing 
several  entries  of  the  same  Doppler  shift  rather  than  several  entries  eaeh  of  a  different 
Doppler  shift.  If  the  eontrol  software  read  and  exeeuted  eaeh  Doppler  shift  entry 
(regardless  of  whether  they  were  all  the  same)  and  produced  the  eorresponding  Doppler 
shift  (whieh  eould  be  observed),  there  is  no  reason  to  believe  the  eontrol  software  would 
not  read  and  exeeute  each  Doppler  shift  entry  (when  they  are  different)  to  potentially 
produee  the  eorresponding  varying  Doppler  shifts  that  are  desired  in  the  same  manner  as 
was  demonstrated  for  the  test  fixture.  Based  on  this,  the  ACS  was  deemed  to  have 


successfully  produeed  a  known  veloeity  spread. 
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When  tested  with  the  MWR-05XP  weather  radar  this  should  be  suffieient  to 
remove  system  bias  from  the  refleetivity  estimate  and  to  eheek  the  aeeuraey  of  mean 
radial  velocity  and  velocity  spread  estimates. 

B,  FUTURE  WORK 

The  first  item  for  future  work  is  to  field  test  the  active  calibrator  system  (ACS) 
with  the  MWR-05XP.  It  should  be  tested  to  see  if  the  constant  Doppler  generated  by  the 
ACS  is  observed  by  the  MWR-05XP.  Then,  a  variable  Doppler  should  be  generated  by 
the  ACS  and  to  see  if  the  MWR-05XP  observes  a  Doppler  spread.  It  should  also  be 
tested  to  see  if  the  equivalent  RCS/Reflectivity  measured  is  observed  as  expected  and 
then  compared  to  the  radar  calibration  constants  obtained  from  the  balloon  and  sphere 
calibration  system. 

The  active  calibrator  system  was  not  entirely  fabricated  when  this  thesis  research 
was  concluded.  The  programmable  attenuator  had  not  yet  been  received.  Once  received, 
the  programmable  attenuator  will  need  to  be  installed.  New  RF  cables  will  have  to  be 
fabricated  to  route  the  RF  path  sequentially  through  the  phase  shifter  and  programmable 
attenuator.  A  control/power  cable  will  have  to  be  fabricated  to  link  the  CRIO’s  DIO 
module  and  the  DC  power  supply  voltages  from  the  voltage  regulator  circuit  with  the 
programmable  attenuator.  This  cable  will  be  the  same  type  (25-pin  to  15-pin)  as  the  one 
used  for  the  same  purposes  by  the  phase  shifter  differing  only  by  pin  assignments.  The 
parts  to  fabricate  this  cable  have  already  been  procured  and  received. 

The  control  software  for  the  programmable  attenuator  has,  for  the  most  part, 
already  been  written.  The  control  software  to  independently  control  the  programmable 
attenuator  is  the  same  as  for  the  phase  shifter.  It  need  only  be  duplicated  with  minor 
modifications  (e.g.  filename  for  Reflectivity  values  vs.  filename  for  Doppler  values, 
MATLAB  script  to  convert  Reflectivity  to  amplitude  shifts  vs.  Doppler  values  to  phase 
shifts,  etc.).  All  of  the  Lab  VIEW  code  could  be  duplicated  with  minor  modifications  and 
run  separately  to  control  the  programmable  attenuator  in  the  same  manner  as  the  phase 
shifter  (both  8-bit).  In  its  current  state,  the  ACS  can  calibrate  the  MWR-05XP  for  a 
single  constant  amplitude  return  (Reflectivity)  and  multiple  variable  Doppler  values. 
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With  the  addition  of  duplicate  control  code  and  the  addition  of  the  programmable 
attenuator,  the  ACS  can  be  configured  to  calibrate  the  MWR-05XP  for  multiple  variable 
amplitude  returns  to  simulate  the  changes  in  pulse  amplitude  that  would  be  observed. 
However,  controlling  these  values  independently  (one  and  then  the  other)  does  not  allow 
a  calibrator  to  simulate  “real  world”  weather  conditions.  In  order  to  calibrate  the  radar 
with  “real  world”  weather  conditions  amplitude  and  phase  values  must  be  controlled 
simultaneously  and  dependently.  To  do  this,  the  control  code  for  both  the  programmable 
attenuator  and  phase  shifter  must  be  timed  together.  This  can  be  achieved  by  controlling 
the  programmable  attenuator  and  phase  shifter  from  within  the  same  control  code.  This 
is  slightly  less  straight  forward  than  the  minor  modifications  for  independent  control  of 
the  two  devices.  However,  it  is  possible  and  desirable.  The  existing  phase  shifter  control 
code  would  have  to  be  modified  to  contain  duplicate  control  subroutines  for  the 
programmable  attenuator  from  within  its  own  control  code.  The  difference  would 
primarily  be  the  additional  use  of  nested  loops  and  sequence  frames  from  with  the  same 
control  code  (LabVIEW  VI).  With  implementation  of  this  simultaneous  dependent 
control  “real  world”  weather  signal  statistics  could  be  simulated  by  the  phase  and 
amplitude  variation  of  the  ACS. 

Testing  of  these  “real  world”  weather  patterns  with  the  MWR-05XP  is  the  final 
objective.  However,  the  ACS  could  and  should  be  tested  incrementally  with  the  MWR- 
05XP  at  various  times  as  new  features  are  implemented.  It  is  suggested  that  the  MWR- 
05XP  and  ACS  be  tested  together  first,  “as-is”  (with  the  phase  shifter  installed  only),  then 
after  the  addition  of  the  programmable  attenuator  and  added  ability  to  control  both  the 
phase  shifter  and  programmable  attenuator  simultaneously.  Final  testing  would  occur 
after  all  features  are  fully  implemented  and  integrated  into  the  ACS  (to  include  output  of 
amplitude  and  Doppler  statistics  of  “real  world”  weather  patterns  from  the  Weather 
Signal  Simulator). 

In  summary  outline.  Future  Work  should  include; 

1 .  A  test  of  the  ACS  with  the  MWR-05XP  in  its  current  state  (at  the  close  of 
this  thesis  research). 
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a.  To  include  observations  of  the  first  and  seeond  Doppler  moments, 
wind  speed  and  wind  shear  a.k.a.  radial  veloeity  and  veloeity 
spread  a.k.a.  a  constant  Doppler  and  a  variable  Doppler. 

b.  To  include  an  observation  of  the  zeroth  Doppler  moment,  a.k.a. 
reflectivity  (to  include  use  of  the  radar  calibration  constant)  or  a 
single,  eonstant  amplitude  return  (corresponding  to  the  equivalent 
RCS  of  the  ACS). 

c.  To  include  a  comparison  of  the  observation  made  in  (b.)  to  the 
observation  made  with  a  balloon  and  sphere  calibration  system. 

2.  Ineorporation  of  the  phase  invariant  programmable  attenuator,  both  the 
hardware  and  independent  software  control  (i.e.  not  affecting  the  digital 
phase  shifter’s  software  control). 

3.  Creation  of  dependent  software  eontrol  of  both  the  programmable 
attenuator  and  phase  shifter  (i.e.  control  of  phase  and  amplitude 
simultaneously). 

4.  A  test  of  the  ACS  with  the  MWR-05XP  in  its  advanced  state  (at  the  close 
of  Future  Work  (2.)  and  (3.)). 

a.  To  include  observations  of  the  zeroth  Doppler  moment, 
reflectivities  (to  include  use  of  the  radar  calibration  constants)  or 
multiple,  variable  amplitude  returns  (corresponding  to  several 
equivalent  RCSs  of  the  ACS). 

b.  To  inelude  observations  of  all  three  Doppler  moments 
simultaneously  (simultaneous  phase  and  amplitude). 

5.  Modifieation  of  the  Weather  Signal  Simulator  to  interact  directly  with  the 
ACS  control  software  to  output  “real  world”  weather  patterns  for  the 
MWR-05XP  (within  hardware  and  software  limitations,  both  for  the 
MWR-05XP  and  ACS  -  e.g.  ambiguous  velocities  due  to  PRF  and 
switehing  speed  of  ACS  hardware). 
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6.  A  test  of  the  ACS  with  the  MWR-05XP  in  its  final  state  (at  the  elose  of 
Future  Work  (5.)).  Are  “real  world”  weather  signal  statisties  observed? 

In  elosing,  it  should  also  be  noted  that  if  the  ACS  is  used  at  ground  level  (the  most 
likely  ease),  a  ground  reflection  will  occur.  The  signal  received  by  the  radar  will  then  be 
determined  by  the  amplitudes  and  relative  phase  between  the  direct  (line-of-sight)  and 
ground  reflected  signal  components  as  influenced  by  the  ground  characteristics.  The 
two-way  propagation  gain  is  12  dB  if  the  ground  reflected  signal  is  equal  in  amplitude 
and  in  phase  with  the  direct  signal.  This  effect  has  been  observed  in  measurements  made 
with  another  calibration  system  subsequent  to  the  completion  of  the  research  work 
described  here.  Compensating  for  this  effect  complicates  the  calibration  procedure  and 
will  require  that  the  ACS  antenna  be  moved  through  a  sequence  of  heights  above  ground 
to  determine  locations  where  the  direct  and  ground  reflected  signals  are  in  phase  and  180 
degrees  out  of  phase.  Signal  measurements  at  these  two  locations  will  permit 
determination  of  the  propagation  gain  for  the  case  where  the  direct  and  ground  reflected 
signals  are  in  phase. 
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APPENDIX  A.  ACS  PICTURES 


Figure  12.  Calibrator  Front  View  with  RAM  (straight  on). 


Figure  13.  Calibrator  Front  View  with  RAM  (looking  down  angle). 


Figure  14.  Calibrator  Left  Side  View  (straight  on). 
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Figure  15.  Calibrator  Left  Side  View  (looking  down  angle). 
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Figure  17.  Calibrator  Left  Side  View  (looking  down  angle),  Zoomed  -  Top  Cover 

OFF. 
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Figure  18.  Calibrator  Top  View  (straight  on),  Zoomed  -  Top  Cover  OFF. 


Figure  20.  Calibrator  Inside  View  (Power  Supply  Circuit). 


Figure  21 .  Calibrator  Inside  View  (Phase  Shifter). 


58 


Figure  22.  Calibrator  Inside  View  (Power  Supply). 
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Figure  23.  Calibrator  Right  Side  View  (looking  down  angle)  Zoomed,  Top  Cover 

OFF 
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Figure  24.  Calibrator  Right  Side  View  (looking  from  above  at  external  in-line 

amp). 


Figure  25.  Calibrator  Rear  Left  Side  View  (looking  down  angle). 
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Figure  26.  Right  Horn  Antenna  (close-in  view,  polarization  is  horizontal). 
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Figure  27.  12”  Sphere  (used  to  calibrate  the  ACS). 
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APPENDIX  B.  DIGITAL  PHASE  SHIFTER  DATA 


Digitally  Controlled  Phase  Shifter 
Manufacturer;  G.T.  Microwave,  Inc. 

Model  Number;  P8P-68N-5 

Electrical  Specifications  for  Digitally  Controlled  360°  Phase  Shifters; 
Frequency  Range  (GHz);  6.0-18.0 
Phase  Error  Vs  Freq.  Max;  +/-  10.0° 

Amplitude  Balance  Max;  +/-  1.50  dB 
Insertion  Foss  Max;  12.0  dB 
V.S.W.R.  Max;  1. 90; I 
Switching  Speed  nSec  Max;  500 
RF  Input  Power  dBm  CW;  +15 
RF  Input  Power  dBm  Max;  +20 


Power/Fogic  Connections; 


No.  of 

FOGIC  PIN 

+I5V 

-I5V 

GND 

BITS 

ASSIGNMENTS 

PIN 

PIN 

PIN 

8 

F.S.B.  @  I  to 
M.S.B.  @  8 

13 

14 

15 

Table  6.  Power  and  Fogic  Connections  for  Digitally  Controlled  Phase  Shifter. 
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Figure  28.  Specifications  Sheet  for  Digital  Phase  Shifter  from  G.T.  Microwave. 


RF  Measurements  on  Digital  Phase  Shifter;  See  Appendix  on  Scalar  and  Vector  Network 
Analyzer  Measurements. 
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APPENDIX  C.  SCALAR  AND  VECTOR  NETWORK  ANALYZER 

MEASUREMENT  DATA 


1 og  MOG 


MARKER  1 
^3.(201  GHz 
'  35.705  dS 


,9.5015  GHz: 
■  35.809  dE 


i MARKER  3 
10.002  GHz 
!  35.059  dE 


►021 

REF  35.0  dB 
2  0.5  dB/ 

V  35 . 809  dB 


START  8 . 000000000  GHz: 
STOP  12.400000000  GHz: 


.„.J  1 - 


!  31  JUL  071 
L11:24:27 


Figure  29.  RF  Amplifier  Gain  Single  VNA  Plot. 
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Figure  30.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  all  pins.  Active  Low.  All  bits  “1”. 

(11111111). 
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Figure  3 1 .  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  MSB  pin  only,  five  volts  others. 

(10000000). 
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S^i  lag  MOG  ►Sp^-^MS  phase 

REF  -10.0  dB  REF  B8.5  - 

0.2  dB/  300.0  m°y 


STRRT  B. 000000000  GHz 
STOP  10.000000000  GHz 


27'jLjL.  0- 
17:21:46 


Figure  32.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  second  MSB  pin  only,  five  volts  others. 

(01000000). 
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Figure  33.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  third  MSB  pin  only,  five  volts  others. 

(00100000). 
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S2;l_  log  MPiG  ►Sg^-^MS  phase 

REF  -10.0  dB  REF  SI . 5  ° 

0.2  dB/  200.0  m°/ 


hp 

¥; 

c 


STftRT  B. 000000000  GHz: 
STOP  10 . 000000000  GHz 


27  JUL  0- 
17: 15:34 


Figure  34.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  fourth  MSB  pin  only,  five  volts  others. 

(00010000). 
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SgjL  log  mog 

REF  -10.0  dB 
0.S  dB/ 


phase 

REF  10.5  « 

200 .0  m  °/ 


4^ 

c 


STPlRT  B. 000000000  GHz 
STOP  10 . 000000000  GHz 


27  JUL  0' 
17:12:06 


Figure  35.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  fifth  MSB  pin  only,  five  volts  others. 

(00001000). 
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Figure  36.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  sixth  MSB  pin  only,  five  volts  others. 

(00000100). 
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Sgj_  lag  MOG  ►Sg^-^MG  phase 

REF  -10.0  dB  REF  2.5  “ 

0.2  dB/  100.0  m°/ 

C 

U 

L 

r 

I 

1  p-  n 

! 

r 

1 

^  , 

r  1 

L  /I 

1 

v  n 

An 

/I  j  LLji 

— 1 

—j 

. . 

mt 

fi 

L  nL 

'  '“1 

1 

-rp| 

mr 

u 

fW  1 

li  ll 

UiH 

,u 

n 

n; 

f 

y 

H 

Ml 

i 

1 

- 1 

r„r  „ 

1 

L  

STRRT  B. 000000000  GHs: 

STOP  10 . 000000000  GHz: 

r'27' JUL  0^ 
[_17 :  03 :  SSj 

Figure  37.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  seventh  MSB  pin  only,  five  volts  others. 

(00000010). 


75 


c 

D 

S^i  log  MPG  ►Sg^/MS  phase 

REF  -10.0  dB  REF  1.5 

0.2  dB^  100.0 

/ip 

j 

i 

— 

— P 

1 

j 

— 1 — 

— 1 — 

1 

(Ul 

il  J 

i  ill 

TT7 

ww 

rf 

STRRT  9.000000000  GHz 

STOP  10 . 000000000  GHz 

!  27“ JUL  07 
i  16:55:37 

Figure  38.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Zero  volts  applied  to  LSB  pin  only,  five  volts  others. 

(00000001). 
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Sgi  log  MRG  ►Spi/MG  phase 

REF  -10.0  dB  REF  0.0  “ 

0 . 2  dB/  100 .0  m  °/ 

hp 


•El?  JUL  07 
1S:5B:49 


STRRT  B . 000000000  GHz 
STOP  10 . 000000000  GHz 


Figure  39.  Digital  Phase  Shifter  Insertion  Loss  (Left)  and  Phase  Value  (Right)  Dual 
VNA  Plots.  Five  volts  applied  to  all  pins.  (00000000). 
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c 

log  MHU  ►Sig  lag  MPG 

REF  0.0  dB  REF  0.0  dB 

10.0  dB.-"  10.0  dB^ 

/hp 

Qr-pi  p 

10.0  dB/div 

OtJOCOOCO  oT/^E 

D 

n 

r 

STRRT  9.000000000  GHz 

STOP  10 . 000000000  GHz 

a7  JUL  07] 
14:45:58  ! 

Figure  40.  Digital  Phase  Shifter  Forward  (left)  and  Reverse  (right)  Insertion  Loss 
Dual  VNA  Plots.  Demonstration  of  deviee’s  bidireetional  property.  Losses 
in  the  forward  and  reverse  are  equal.  Five  volts  applied  to  all  pins. 

(00000000). 
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Figure  43.  Digital  Phase  Shifter  SNA  Plot.  Fixed  Phase  State,  no  Doppler  Shift. 
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Figure  44.  Digital  Phase  Shifter  SNA  Plot.  Periodic  Rotation  through  Phase  States, 

Doppler  Shift. 
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Figure  46.  Digital  Phase  Shifter  Phase  Value  Single  VNA  Plot.  Zero  volts  applied 
to  second  MSB  pin  only,  five  volts  others.  (01000000)  90°  Phase  Shift. 
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APPENDIX  D.  PROGRAMMABLE  ATTENUATOR  DATA 


Phase  Invariant  Programmable  Attenuator 
Manufaeturer;  American  Microwave  Corporation 
Model  Number;  PI-DVAN-8012-8 

Electrical  Specifications: 

Frequency;  8  GHz  to  12  GHz 
Range/Resolution;  60  dB/0.25  dB 
Insertion  Loss:  4.0  dB  Typical 
Switching  speed;  300  ns  Maximum 

Maximum  Power  (No  Damage);  +27  dBm  (operating  +10  dBm) 
Phase  Variation/Shift;  0  to  20  dB  ;  +/-5° 

20  to  40dB  ;+/-10° 

40  to  60  dB  ;  +/-20° 


Power/Logic  Connections: 


No.  of 

LOGIC  PIN 

+15V 

-15V 

GND 

BITS 

ASSIGNMENTS 

PIN 

PIN 

PIN 

8 

L. S.B.  @  5  to 

M. S.B.  @  12 

13 

14 

4 

Table  7.  Power  and  Logic  Connections  for  Programmable  Attenuator. 
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ccaaows 


DESCaiPTlOK 

*MC  UOOEL  ft-OVAN-8032-8  S  A  HIGH  SPEED.  lOW  WSEPTION  LOSS. 
phase  invariant.  0  to  60  dG,  S  9iT  OlOiTAL  VARIAei.C.  SOLID  STATE. 

OCTAVr  ewD  ATrCNUATOR  that  OPERATCS  over  the  8  GM2  to  12  GH*  RANGE. 


ORIGINAL  RELEASE 


•  TREOUENCV; . -  . 8  (Wz  TO  17  CHl 

•  HANGE/RESOLUnON; . 60  <IB  /  0.74  <fS 

0  INSERTION  LOSS;-  •  •  . .  A.O  06  TYPOI..  3.80  <JB  COAt,  4.4  UB  MAXIMUM 

«  DIGITAL  CONTROL:' . 8  ST  TTl. 

«  SWfCHlNG  SPEED; . 300  ns  UAXHAJM 

c  maximum  POWER: . . 4-77  dSm  .  {{DERATING  4  10  oeml 

(NO  OAUACE) 

'  0T0I7  ■  . 36  dBm  MINMUM 

•  BIAS  VOlTACSE:- . •■»12V  O  700  mA  MAXIMUM..  '50  mA  TYPICAL 

-SV  O  TOO  mA  UAAMUU.,  74  mA  TYPICAL 

•  PHASE  VMPIATION/SHfT;  qj 

O  TO  20  <18 . 4¥  (CCMt  *4-)  4 

30  TO  40  «B . ilC  (CCML  ±8') 

40  TO  €0  dS . iTO’  (GOAL  tIS’]  -j— 

«  Tew>ewruRe . .  .q-c  to  tot 

•  RE  CONflECTORS: . SWA  fEMAUE 

•  BIAS  VOLTAGE 

AND  CONTROL .  . 13  PK  MULTlPIH  CONNECTOR 

•  SI2C:  . . -2.00  (L)  *  2.60  (W)  x  O.SO  (H)  2 

«  WEIGHT . . -<  4  OUNCES 


13  PIN  MICRO  D 
MULTI  PIN  CONtCCTOR 


«008?  THRU 
4  PLACES- 


SOL®  STATE 
VARIABLE  ATTENUATOR 
JI  unon  MO.  fl-PWM-eoig-e 


(M5  "l 

0,095  —  U- 


'  SNA  CF>  REMOVABLE 
-RF  CONNECTOR 
2  PLACES 


CONFIDENTIAL  AND  PROPRIETARY 


ENVIRONMENTAL  RATINGS: 


,  ... _ _  <rc  TO  1-70  C  (OPCRATiNC) 

-65-C  TO  -frUSC  <STORA&) 

•  HCMiOnY<  -  .  ,  .  ...  Mit-STO-20aF,  MOHOO  1 030  CONO.  8 

•  SHOCKS  .  .  ,  .  . M(L-SI0~2O2F,  MCTmOO  2136  CONO.  0 

»  VlGflATlCNi  ....  - . MIL-ST0-202E.  METHOD  7040  CONO.  8 

»  ALTITUDE: . MiL-ST0-702P.  METHOD  I05C  COND.  8 

•  TEMPERATUAE  CYCLE'  -  -  •  U:L-STO-2<12F.  HCTHOO  107D  CONO-  A 


ALL  DIMENSIONS  >IPE  IN  INCHES 
TOLERANCES; 

X.XX  rtO.QTO 
X.XXX  ±0010 


f  I  AMERICAN  MICROWAVE  CORPORATION 
LcRI  PRQDCRICK,  MARYLAND 


“J  PRODUCT  fEATuRE 

PI-0VAN-S012-8 

SOLID  state:  phase  INVA^ANT  ATTENUATOR 


Figure  47.  Specifications  Sheet  for  Programmable  Attenuator  from  American 

Microwave. 
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APPENDIX  E.  ANTENNA  DATA 


Antenna  List: 

(2)  Pasternack  Sectoral  Horn  Antennas  (used  as  part  of  ACS) 

(2)  Unknown  X-band  Quad-Ridge  Dual-Pol  Horn  Antennas  (used  as  part  of  anechoic 
chamber  ACS  testing) 

(1)  Narda  Standard  Gain  Horn  Antenna  (used  as  part  of  anechoic  chamber  antenna 
testing) 
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pASTtnMACK  eMTEmPHoa.  me. 


PO  Box  16759.  Irvine.  Ca  92623 


Toll  Free:  (866)  726-8375 
Direct:  +1  (949)  261-1920 
FAX:  +1  (949)261-7451 
Email:  techsupport@pastemack.com 


PE9850  thru  PE9864 
Standard  Gain  Horns 


Standard  Gain  Horns 

P/N 

Size 

Frequency 

(GHz) 

Gain 

(dB) 

3  dB  Width 
E-Plane  H-Plane 

PE9850-10 

WR28 

26.5-40.0 

10 

54.2” 

54.4” 

PE9850-15 

WR28 

26.5-40.0 

15 

32.1° 

31.3” 

PE9850-20 

WR28 

26.5-40.0 

20 

16.7° 

18.3” 

PE9851-10 

WR34 

22.0  -  33.0 

10 

54.1” 

53,2” 

PE9851-15 

WR34 

22,0-33.0 

15 

23.1° 

40.8” 

PE9851-20 

WR34 

22.0-33.0 

20 

17.0” 

17.4” 

PE9852-10 

WR42 

18.0-26.5 

10 

58.0” 

57.0” 

PE9852-15 

WR42 

18.0-26.5 

15 

31.3” 

31.5” 

PE9852-20 

WR42 

18.0-26.5 

20 

17.5” 

17.8” 

PE9853-10 

WR51 

15.0-22.0 

10 

55.1” 

54.2” 

PE9853-15 

WR51 

15.0-22.0 

32.0” 

31.8” 

PE9853-20 

WR51 

15.0-22.0 

20 

16.9” 

18.0” 

PE9854-10 

WR62 

12.4-18.0 

10 

55.3” 

50.9” 

PE9854-15 

WR62 

12.4-18.0 

15 

30.1” 

31.2” 

PE9854-20 

WR62 

12,4-18.0 

20 

18.8” 

18.9” 

PE9855-10 

WR75 

10.0-15.0 

10 

50.2” 

49.2” 

PE9855-15 

WR75 

10.0-15,0 

15 

35.4” 

28.5” 

PE9855-20 

WR75 

10.0-15.0 

20 

16,3” 

17.2” 

PE9856-10 

WR90 

8.20-12.4 

10 

48.5” 

47.4” 

PE9856-15 

WR90 

8.20-12.4 

15 

29,3” 

29.0” 

PE9856-20 

WR90 

8.20-12.4 

20 

16.1” 

16.5” 

PE9857-10 

WR102 

7.00-11.0 

10 

55.5” 

54,1” 

PE9857-15 

WR102 

7.00-11.0 

15 

29.6” 

29.3” 

PE9857-20 

WR102 

7.00-11.0 

20 

17.0” 

16.7” 

PE9858-10 

WR112 

7.05-10.0 

10 

56.8” 

55.2” 

PE9858-15 

WR112 

7.05-10.0 

15 

32.4” 

32.0” 

PE9858-20 

WR112 

7.05-10.0 

20 

19.3” 

19.3” 

PE9859-10 

WR137 

5.85  -  8.20 

10 

55.1” 

54.2” 

PE9859-15 

WR137 

5.85  -  8.20 

15 

33.7” 

33,2” 

PE9859-20 

WR137 

5.85  -  8.20 

20 

18.7” 

18.8” 

PE9860-10 

WR159 

4.09-7.05 

10 

59,8” 

48.3” 

PE9860-15 

WR159 

4.09-7.05 

15 

31.3” 

30.8” 

PE9860-20 

WR159 

4.09  -  7.05 

20 

14.3” 

16.9” 

PE9861-10 

WR187 

3.95-5.85 

10 

55.0” 

54.1” 

PE9861-15 

WR187 

3.95-5.85 

15 

33.8” 

33.3” 

PE9861-20 

WR187 

3.95-5.85 

20 

18.9” 

19.2” 

PE9862-10 

WR229 

3.30  -  4.90 

10 

58.6” 

51.9” 

PE9862-15 

WR229 

3.30  -  4.90 

.  15 

32.7” 

32.4” 

PE9862-20 

WR229 

3.30-4.90 

20 

17.1” 

16.7” 

PE9863-10 

WR284 

2.60-3.95 

10 

50.8” 

54.1” 

PE9863-15 

WR284 

2.60  -  3.95 

15 

31.0” 

30.6” 

PE9863-20 

WR284 

2.60  -  3.95 

20 

17.2” 

16.5” 

PE9864-10 

WR430 

1.70-2.60 

10 

64.8” 

45.4” 

PE9864-15 

WR430 

1.70-2.60 

33  1® 

32.0” 

PE9864-20 

WR430 

1.70-2.60 

20 

17^3^ 

17.4” 

- lypicsl  E-Plane  Pattern 

- Typical  H-Plane  Pattern 

Typical  10dB  Gain 


typical  E-Plane  Pattern 
- typical  H-Plane  Pattern 


Typical  15dB  Gain 


- Typical  E-Plane  Pattern 

■"  —  Typical  H-Plane  Pattern 


Typical  20dB  Gain 


Gain  and  beamwidth  data  are  typicals. 

Additional  data,  such  as  calibration  data,  is  provided  for  an  additional  fee. 


Figure  48.  Pasternack  Sectoral  Horn  Antenna  Data  Sheet. 
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Figure  49.  Pasternack  Sectoral  Horn  Antenna  Outline  Drawing. 


89 


90 


NOTES: 

1.  UNLESS  OTHERWISE  SPECIFIED  ALL  DIMENSIONS  ARE  NOMINAL. 

2.  ALL  SPECIFICATIONS  ARE  SUBJECT  TO  CHANGE  WITHOUT  NOTICE  AT  ANY  TIME. 


2.60-40  GHz 

Standard  Gain  Horns 

•  Primary  Standard  of  Antenna  Gain 
^  j  lyiQjjgIg  2.60  GHz  to  40  GHz 


FREQUENCY 

RANGE 

(GHz) 

MODEL 

SAND 

WAVEGUIDE 

SIZE 

INPUT 

COVER  FLANGE 
EQUIVALENT 

VSWR 

(Max) 

WEIGHT 

Iba,  Kg 

2.60-3.95 

644 

S 

WR-284 

UG-584/U 

1.15 

6 

2,8 

:  .-h.  ..  . 

643 

c 

WR-187 

UG-407/U 

1.15 

2.3 

1.1 

5.4-8.20 

642 

XN 

WR-137 

UG-441/U 

1.15 

1.0 

0.5 

8.20-12.4 

640 

X 

WR-90 

UG-135/U 

1.15 

0.5 

0,23 

12.4-18 

639 

KU 

WR-62 

UG-419/U 

1.15 

0.2 

0.1 

18-26.S 

638 

WR-42 

UG-595/U 

1.15 

0.2 

0,1  : 

26.5-40 

V637 

V 

WR-28 

UG-599/U 

1.15 

0.1 

0.05 

Patterns  for  al  models  in  this  series  conform  to  the  following 
description;  Beam  width  in  E  and  H  plane  varies  from  23''  at  the 
highest  frequency  to  SA-'-'  at  the  lowest  frequency.  Side  lobes  in 
the  H  plane  are  all  more  than  20  dB  down.  First  side  oPes  in  the 
E  p  ane  are  13  dB  down,  second  side  lobes  are  18  dB  down  and 
all  other  E  plane  iooes  are  more  than  20  dB  down. 

Gain  at  Mid  Frequency;  1 6.5  dB  (with  reference  to  isotropic  radiation) 
variation  is  1 .5  dB  over  tota'  band  about  the  mid  band  va  ue. 


narda 


435  MORELAND  ROAD  -  HAUPPAUQE.  NY  11788 
USA  "EL:  (1)  631  231-1700  •  INT'LTEL:  (1)  63  1  231-1390  •  FAXri)  631  231-1711 
E-MAIL:  naraaeasT©L-3COM,Gom  *  wvvw.narGamicrowave.com 


Figure  5 1 .  Narda  Standard  Gain  Horn  Antenna  Datasheet  (page  1). 
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Waveguide  Horns 


Standard  Gain  Horn 


. i-A . i-i . 

. [-■r-B-- . 

' :  'MAX:' '  7  ' ':' 

. . 

1  'MAX  : . 

V637 

1.76 

1.06 

.82 

(44,7) 

(26.9) 

(20.8) 

636 

2.57 

1.51 

1.16 

(36.3) 

(29.4) 

639 

3.48 

2.20 

1.73 

/3) 

(55.8) 

(43.9) 

640 

5.06 

3.09 

3.34 

(128.5) 

(78.4) 

(53,4) 

642 

7,76 

4,67 

3.53 

(197.1) 

(1 1 8.6) 

(83,6) 

B43 

10.47 

6.34 

4.80 

(265.9) 

(161.0) 

(121.9) 

644 

.62 

9.52 

7.16 

(401.8) 

(241.8) 

(181,8) 

Dimensions  in  oarentheses  are  in  miilimeters  and  tor  reference  ony. 
■  ■ :  ■  A'l  dimensions  max,  un;ess  othem/ise  specified/  ■ : 


435  MORELAND  ROAD  •  HAUPPAUQE.  NY  1 1 78S 

USA  TEL:  (I)  631  231-1700  ®  !N~LTEL:  (1)  631  231-1330  •  FAX(I)  631  231-1711 
E-MAIL:  nRrdReafit©L-3COM,Qorr  *  wvvv.'.narnaiTiicroy.'ave.com 


narda 


Figure  52.  Narda  Standard  Gain  Horn  Antenna  Datasheet  (page  2) 


WAVEGUIDE  TO  COAXIAL  ADAPTERS  FOR  STANDARD  GAIN  HORN 


STANDARD  GAIN 

HORN  MODEL 

644 

;  643  : 

WAVEGUIDE  TO  COAXIAL  TYPE 

N  -  FEMALE  MODEL 

61 4A 

61 3A  ;  : 

ADAPTER  SIWA/3.5/2.9 
FEMALE  MODEL 

642 

61 2A 

4602 

. . 

■  :  ■  :  60TA :  ■ 

:  :  :  4601]  :  1  :  1  : 

639 

609 

4609 

.  .  . .  .  .  . 

.  .  . .  . 

V637 

V4607 

narda 


435  MORELAND  ROAD  -  HAUPPAUQE.  NY  11788 
USA  "EL:  (1)  631  231-1700  •  INT'LTEL:  (1)  63  1  231-1390  •  FAXri)  631  231-1711 
E-MAIL:  naraaeasT©L-3COM,Gom  *  www.naraamicrowave.coifi 


Figure  53.  Narda  Standard  Gain  Horn  Antenna  Datasheet  (page  3). 


93 


ABSOLUTE  OAiN  OALieBATIQN 

H0OiL  S4S  STANDARD  6AIN  H0RN 


t.f  i.p  i».e  11.9  11.9  i:.e  is.? 

fHEGUEHCT  fOnZI 


AESOLUTS  OAEN  OALiSRAtlON 
NARDA  MODEL  €43  STANDARD  QAtN  HORN 


i^nEsusHCT  laos 


ABSOLUTS  SA3M  CAUBBATE9R 
NAN  DA  MODEL  @4S  STANDARD  OAIN  HORN 


pRsanEfseY  faHsi 
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Figure  54.  Narda  Standard  Gain  Horn  Antenna  Datasheet  (page  5). 
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APPENDIX  F.  ANTENNA  PATTERN  DATA 


A  total  of  five  antennas  were  used  in  testing.  One  Narda  standard  gain  horn 
antenna  was  used  during  antenna  testing  (patterns)  only.  Two  X-band  horn  antennas 
were  used  during  antenna  testing  (patterns),  then  used  during  ACS  testing  (as  the  transmit 
and  receive  antennas  of  the  radarAfiSlA).  Two  Pasternack  horn  antennas  were  used 
during  antenna  testing  (patterns),  then  used  during  ACS  testing  (as  the  transmit  and 
receive  antennas  of  the  ACS).  Measurements  were  made  at  three  frequencies:  9370, 
9680,  and  9990  MHz  (low,  mid,  and  high  frequency  points  on  the  MWR-05XP  weather 
radar’s  RF  band).  E-plane  and  H-plane  patterns  were  measured.  To  measure  both 
planes,  the  antenna  under  test  was  first  rotated  360°  in  azimuth  with  the  radiator  oriented 
for  horizontal  polarization  (E-plane),  then  rolled  90°  and  the  process  repeated  with  the 
radiator  now  oriented  for  vertical  polarization  (H-plane).  In  total,  two  plane  patterns 
were  measured  for  five  antennas  at  three  frequencies  to  create  thirty  data  sets. 

A,  VNA  ANTENNA  PATTERN  DATA 

In  all  plots,  solid  line  is  for  X-band  Horn  Antenna  #1  (“xbh”)  and  dotted  line  is 
for  X-band  Horn  Antenna  #2  (“ybh”).  The  last  three  plots  are  solid  line  for  the  Standard 
Gain  Horn  (“sgh”)  which  are  included  to  provide  a  gain  reference. 


95 


Figure  55.  X-band  Florn  Antenna  #1  Pattern  at  9370MFIz  (E  plane  in  red  and  FI 

plane  in  blue). 


Figure  56.  X-band  Flom  Antenna  #2  Pattern  at  9370MFIz  (E  plane  in  red).i2 


12  Note  because  the  measured  values  slightly  exceeded  -15  dB  the  gain  values  have  been  normalized  to 
0  rather  than  -15. 
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Figure  57.  X-band  Horn  Antenna  #2  Pattern  at  9370MHz  (H  plane  in  blue). 


Figure  58.  X-band  Horn  Antenna  #1  Pattern  at  9680MHz  (E  plane  in  red  and  H 

plane  in  blue). 
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Figure  59.  X-band  Horn  Antenna  #2  Pattern  at  9370MHz  (E  plane  in  red  and  H 

plane  in  blue). 


Figure  60.  X-band  Horn  Antenna  #1  Pattern  at  9990MHz  (E  plane  in  red  and  H 

plane  in  blue). 
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Figure  61 .  X-band  Horn  Antenna  #2  Pattern  at  9990MHz  (E  plane  in  red  and  H 

plane  in  blue). 


Figure  62.  Standard  Gain  Horn  Antenna  Pattern  at  9370MHz  (E  plane  in  red  and  H 

plane  in  blue). 
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Figure  63.  Standard  Gain  Horn  Antenna  Pattern  at  9680MHz  (E  plane  in  red  and  H 

plane  in  blue). 


Figure  64.  Standard  Gain  Horn  Antenna  Pattern  at  9370MHz  (E  plane  in  red  and  H 

plane  in  blue). 
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B. 


ACS  ANTENNA  PATTERN  DATA 


In  all  plots,  solid  line  is  for  Pasternack  Horn  Antenna  #1  (a.k.a.  “xpe”)  and  dotted 
line  is  for  Pasternack  Horn  Antenna  #2  (a.k.a.  “ype”).  The  last  three  plots  are  solid  line 
for  the  Standard  Gain  Horn  (a.k.a.  “sgh”).  The  pattern  plots  for  the  Standard  Gain  Horn 
have  been  included  to  provide  a  gain  reference  and  in  all  three  plots  the  gain  has  been 
adjusted  to  account  for  the  use  of  a  6.5  dB  attenuator.  Also,  note  all  data  have  been 
shifted  -90°  to  boresight  the  patterns  at  0°.  This  was  done  to  make  these  patterns  look 
like  the  other  patterns  previously  presented.  Additionally,  in  all  plots  the  E-plane  pattern 
is  in  red  and  H-plane  pattern  is  in  blue. 


Figure  65.  Pasternack  Horn  Antenna  #1  Pattern  at  9370MHz  (E  plane  in  red  and  H 

plane  in  blue). 
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Figure  66.  Pasternack  Flom  Antenna  #2  Pattern  at  9370MFIz  (E  plane  in  red  and  FI 

plane  in  blue). 


Figure  67.  Pasternack  Horn  Antenna  #1  Pattern  at  9680MHz  (E  plane  in  red  and  H 

plane  in  blue). 
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Figure  68.  Pasternack  Flom  Antenna  #2  Pattern  at  9680MFIz  (E  plane  in  red  and  FI 

plane  in  blue). 


Figure  69.  Pasternack  Horn  Antenna  #1  Pattern  at  9990MHz  (E  plane  in  red  and  H 

plane  in  blue). 
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Figure  70.  Pasternack  Horn  Antenna  #2  Pattern  at  9990MHz  (E  plane  in  red  and  H 

plane  in  blue). 


Figure  7 1 .  Standard  Gain  Horn  Antenna  Pattern  at  9370MHz  (E  plane  in  red  and  H 

plane  in  blue). 
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Figure  72.  Standard  Gain  Florn  Antenna  Pattern  at  9680MHz  (E  plane  in  red  and  H 

plane  in  blue). 


Figure  73.  Standard  Gain  Horn  Antenna  Pattern  at  9990MHz  (E  plane  in  red  and  H 

plane  in  blue). 
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C.  MATLAB  PLOTTING  FILES 

plotpatterns.m _ 


polarDB(xbh9370H(; ,  1  ),xbh9370H(:  ,2),’r') 
polarDB(xbh9370 V(; ,  1  ),xbh9370 V(:  ,2),’b') 

figure 

%%  normalized  ro  to  0  not  -15 

polarDB(ybh9370H(:,l),ybh9370H(:,2),’r;') 

figure 

polarDB(ybh9370V(;,l),ybh9370V(;,2),'b;') 

figure 

polarDB(xbh9680H(;,l),xbh9680H(;,2),V) 

polarDB(xbh9680V(;,l),xbh9680V(;,2),’b') 

figure 

%%  both  normalized  to  0 

polarDB(ybh9680H(:,l),ybh9680H(:,2),'r;') 

polarDB(ybh9680V(;,l),ybh9680V(:,2),’b;’) 

figure 

polarDB(xbh9990H(; ,  1  ),xbh9990H(:  ,2),’r') 
polarDB(xbh9990 V(: ,  1  ),xbh9990 V(:  ,2),’b’) 

figure 

polarDB(ybh9990H(;,l),ybh9990H(:,2),'r;') 
polarDB(ybh9990 V(; ,  1  ),ybh9990 V(;  ,2),'b;') 

figure 

polarDB(sgh9370H(;,l),sgh9370H(;,2),'r-') 
polarDB(sgh9370 V(; ,  1  Xsgh9370 V(;  ,2),'b-') 

figure 

polarDB(sgh9680H(:,l),sgh9680H(;,2),'r-') 

polarDB(sgh9680V(;,l),sgh9680V(;,2),'b-') 

figure 

polarDB(sgh9990H(; ,  1  ),sgh9990H(;  ,2),'r-') 
polarDB(sgh9990 V(: ,  1  ),sgh9990 V(:  ,2),’b-') 


106 


plot  pattems2.m 


%%  -90  correction  for  90  degree  offset  when  patterns  taken  (-90  moves 
%%  boresight  to  0  degrees) 
polarDB(xpe9370H(;,l)-90,xpe9370H(:,2),'r') 
polarDB(xpe9370V(;,l)-90,xpe9370V(;,2),'b') 

figure 

polarDB(ype9370H(:,l)-90,ype9370H(;,2),'r;') 

polarDB(ype9370V(;,l)-90,ype9370V(;,2),'b;') 

figure 

polarDB(xpe9680H(;,l)-90,xpe9680H(;,2),'r') 

polarDB(xpe9680V(;,l)-90,xpe9680V(;,2),'b') 

figure 

polarDB(ype9680H(;,l)-90,ype9680H(;,2),'r;') 

polarDB(ype9680V(:,l)-90,ype9680V(:,2),'b;') 

figure 

polarDB(xpe9990H(:,l)-90,xpe9990H(;,2),'r') 

polarDB(xpe9990V(;,l)-90,xpe9990V(;,2),'b') 

figure 

polarDB(ype9990H(;,l)-90,ype9990H(;,2),'r;') 

polarDB(ype9990V(:,l)-90,ype9990V(:,2),'b;') 

%%+6.5  correction  for  6.5dB  attenuator  used  in  sgh  measurement,  not  used 

%%in  pe  measurements 

figure 

polarDB(sgh9370H(:,l),sgh9370H(:,2)-l-6.5,'r-') 

polarDB(sgh9370V(:,lXsgh9370V(;,2)+6.5,’b-') 

figure 

polarDB(sgh9680H(:,l),sgh9680H(;,2)-l-6.5,'r-') 

polarDB(sgh9680V(;,l),sgh9680V(:,2)+6.5,'b-') 

figure 

polarDB(sgh9990H(;,l),sgh9990H(;,2)-l-6.5,'r-') 

polarDB(sgh9990V(:,l),sgh9990V(:,2)+6.5,’b-') _ 
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polarDB.m 


function  polardb(ang,pat,line_style) 

%  NON-NORMALIZING  VERSION  THAT  DOES  NOT  REQUIRE  POEAR2 
%  take  angle  and  dB  data  and  generate  polar  plot 
%  angle  data  in  ang(;)  is  degrees;  dB  data  in  pat(;) 
if  any(size(ang)  ~=  size(pat)) 

error('THETA  and  RHO  must  be  the  same  size.'); 
end 

pmax=max(pat); 

dynr=60;  %  dynamic  range  of  dB  scale 
rticks=4;  %  number  of  rings  in  polar  plot 

dbstep=dynr/rticks;  %  dB  spacing  between  rings  —  must  be  integer 
top=floor(pmax/dbstep)*dbstep+dbstep;  %  convenient  outer  ring  dB  value 
theta=ang*pi/180; 

%  set  a  floor 
for  i=l  ;length(pat) 

if  pat(i)  <  top-dynr,  pat(i)=top-dynr;  end 
end 

rho=p  at+dynr-top ; 
hold  on 

%  define  a  circle 

th  =  0;pi/50;2*pi; 
xunit  =  cos(th); 
yunit  =  sin(th); 

%  now  really  force  points  on  x/y  axes  to  lie  on  them  exactly 
inds  =  [l;(length(th)-l)/4;length(th)]; 
xunits(inds(2;2;4))  =  zeros(2,l); 
yunits(inds(l;2;5))  =  zeros(3,l); 
rmax=dynr;  rmin=0; 
rinc  =  (rmax-rmin)/rticks; 
for  i=(rmin+rinc);rinc;rmax 

plot(xunit*i, yunit*i,'-', 'color', 'k','linewidth',l); 

text(0,i+rinc/10,['  '  num2str(i-dynr-l-top)],'verticalalignment', 'bottom'  ); 
end 

%  plot  spokes 

th  =  (l;6)*2*pi/12; 

cst  =  cos(th);  snt  =  sin(th); 

cs  =  [-cst;  cst]; 

sn  =  [-snt;  snt]; 

plot(rmax*cs, rmax*sn,'-', 'color', 'k','linewidth',l); 

%  annotate  spokes  in  degrees 


108 


rt  =  1 . 1  *rmax; 

for  i  =  1  :niax(size(th)) 

text(rt*cst(i),rt*snt(i),mt2str(i*30),'horizontalalignment', 'center' ); 
if  i  ==  max(size(th)) 

loc  =  int2str(0); 

else 

loc  =  int2str(l  80+1*30); 
end 

text(-rt*cst(i),-rt*snt(i), loc, 'horizontalalignmenf, 'center'  ); 
end 

%  set  viewto  2-D 
view(0,90); 

%  set  axis  limits 

axis(rmax*[-l  1-1.1  1.1]); 

%  transform  data  to  Cartesian  coordinates 
XX  =  rho.*cos(theta); 
yy  =  rho.*sin(theta); 

%  plot  data  on  top  of  grid 

q  =  plot(xx,yy,line_style); 
hpol  =  q; 

axis('equar);axis('off); 
hold  off 
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D.  ABSOLUTE  GAIN  AND  3  DB  BEAMWIDTH  FOR  X-BAND  AND 
PASTERNACK  HORN  ANTENNAS  (CALCULATIONS  AND 
MEASUREMENTS) 


Standard  Gain  Horn 

Data  from  Manufacturer's  Chart 

Frequency  (GHz) 

Gain  (dB) 

9.37 

15.66 

9.68 

15.90 

9.99 

16.15 

Standard  Gain  Horn  (sgh) 

Data  Measured  in  Chamber 

Frequency  (GHz) 

Plane 

Location  (Degrees 

Gain  (dB) 

(Eor  H) 

from  Boresight) 

9.37 

E 

-1 

-12.43 

9.37 

H 

2 

-12.61 

9.68 

E 

0 

-12.07 

9.68 

H 

2 

-12.65 

9.99 

E 

0 

-12.05 

9.99 

‘slight  variation  in  setup 

H 

0* 

-11.75 

X-band  Horn  1  (xbh) 

M=Measured,  C=Calculated 

M 

C 

C 

M 

M 

M 

M 

C 

Frequency  (GHz) 

Plane 

Location  (Degrees 

Relative  Gain 

A  sgh  (dB) 

Absolute 

First 

-3dB  Point 

Second 

-3dB  Point 

3dB  Beamwidth 

(Eor  H) 

from  Boresight) 

(MAX  Value,  dB) 

Gain  (dB)  (Loc,  Deg) 

(Rel  Gain,  dB) 

(Loc,  Deq) 

(Rel  Gain,  dB) 

(Degrees) 

9.37 

E 

2 

-15.18 

2.75 

12.91 

-15 

-18.04 

17 

-18.08 

32 

9.37 

H 

3 

-15.66 

3.04 

12.61 

-21 

-18.64 

24 

-18.62 

45 

9.68 

E 

3 

-15.10 

3.02 

12.88 

-14 

-18.09 

16 

-18.03 

30 

9.68 

H 

2 

-15.11 

2.46 

13.44 

-19 

-18.14 

22 

-18.10 

41 

9.99 

E 

0&1 

-15.60 

3.55 

12.60 

-14 

-18.42 

16 

-18.56 

30 

9.99 

H 

1 

-15.72 

3.97 

12.18 

-18 

-18.62 

22 

-18.73 

40 

X-Band  Horn  2  (ybh) 

9.37 

E 

2 

-14.71 

2.28 

13.38 

-14 

-17.43 

17 

-17.58 

31 

9.37 

H 

3 

-15.02 

2.41 

13.25 

-21 

-18.02 

24 

-18.06 

45 

9.68 

E 

3 

-14.89 

2.82 

13.08 

-13 

-17.81 

17 

-18.04 

30 

9.68 

H 

2 

-14.80 

2.15 

13.75 

-19 

-17.93 

22 

-17.74 

41 

9.99 

E 

1 

-15.12 

3.06 

13.09 

-14 

-18.09 

16 

-18.11 

30 

9.99 

H 

2 

-15.34 

3.59 

12.56 

-18 

-18.38 

22 

-18.29 

40 

Table  8.  Datafile  for  X-band  Horn  Antennas. 
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Figure  74.  Plot  of  Absolute  Gain  for  X-band  Flom  Antennas. 
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Figure  75.  Plot  of  3dB  Beamwidth  for  X-band  Florn  Antennas. 


Standard  Gain  Horn 

Data  from  Manufacturer's  Chart 

Frequency  (GHz) 

Gain  (dB) 

9.37 

15.66 

9.68 

15.90 

9.99 

16.15 

Standard  Gain  Horn  (sgh) 

Data  Measured  in  Chamber  [minus  -6.5dB  attenuator  in  line] 

Frequency  (GHz) 

Plane 

Location  (Degrees 

Gain  (dB) 

(Eor  H) 

from  Boresight) 

9.37 

E 

-1 

-5.93 

-12.43 

9.37 

H 

2 

-6.11 

-12.61 

9.68 

E 

0 

-5.57 

-12.07 

9.68 

H 

2 

-6.15 

-12.65 

9.99 

E 

0 

-5.55 

-12.05 

9.99 

‘slight  variation  in  setup 

H 

0* 

-5.25 

-11.75 

Pasternack  Horn  1  (1PE) 

M=Measured,  C=Calculated 

M 

C 

C 

M 

M 

M 

M 

C 

Frequency  (GHz) 

Plane 

Location  (Degrees 

Relative  Gain 

A  sah  fdBl 

Absolute 

First 

3dB  Point 

Seconc 

-3dB  Point 

3dB  Beamwidth 

(Eor  H) 

from  Boresight) 

(MAX  Value,  dB) 

Gain  (dB)  (Loc,  Deg) 

(Rel  Gain,  dB) 

(Loc,  Deq) 

(Rel  Gain,  dB) 

(Degrees) 

9.37 

E 

92 

-7.24 

1.31 

14.35 

76 

-10.30 

109 

-10.32 

33 

9.37 

H 

93 

-8.85 

2.74 

12.91 

77 

-12.04 

108 

-12.10 

31 

9.68 

E 

92 

-6.87 

1.30 

14.60 

77 

-9.93 

108 

-9.91 

31 

9.68 

H 

92 

-8.16 

2.01 

13.89 

79 

-10.96 

106 

-10.96 

27 

9.99 

E 

91 

-7.02 

1.47 

14.68 

76 

-10.02 

105 

-10.19 

29 

9.99 

H 

93 

-8.43 

3.18 

12.97 

78 

-11.61 

107 

-11.57 

29 

Pasternack  Horn  2  (2PE) 

9.37 

E 

93 

-7.38 

1.45 

14.20 

77 

-10.21 

109 

-10.30 

32 

9.37 

H 

93 

-8.86 

2.74 

12.91 

78 

-11.67 

107 

-11.64 

29 

9.68 

E 

93 

-6.71 

1.14 

14.76 

77 

-9.87 

109 

-9.92 

32 

9.68 

H 

92&93 

-8.07 

1.92 

13.98 

79 

-10.91 

107 

-11.24 

28 

9.99 

E 

93 

-6.85 

1.29 

14.86 

78 

-9.98 

108 

-10.19 

30 

9.99 

H 

93 

-8.17 

2.92 

13.23 

79 

-10.96 

107 

-11.27 

28 

Table  9.  Datafile  for  Pasternack  Flom  Antennas. 
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Figure  76.  Plot  of  Absolute  Gain  for  Pasternaek  Florn  Antennas. 
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3dB  Beamwidth 


-♦-Pasternack  Horn  1  (1  PE)  -  H-Pol/E-Field 
-■-Pasternack  Horn  1  {1PE) -  V-Pol/H-Field 
-♦-  Pasternack  Horn  2  (2PE)  -  H-Pol/E-Field 
Pasternack  Horn  2  (2PE)  -  V-Pol/H-Field 


Frequency  (GHz) 


Figure  77.  Plot  of  3dB  Beamwidth  for  Pasternaek  Florn  Antennas. 


E.  EXPLANATION  OF  DATAFILE  MEASUREMENTS  AND 
CALCULATIONS 

The  first  entry  in  both  datafiles  is  the  absolute  gain  of  the  referenee  horn  antenna, 
the  Narda  standard  gain  horn  antenna.  These  values  were  obtained  from  the  Narda 
datasheet,  page  5.  See  below  for  a  repliea  of  this  graph  with  the  look-up  proeess 
annotated. 
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Figure  78.  Narda  Standard  Gain  Horn  Absolute  Gain  Plot  (zoomed  to  area  of 
interest  and  annotated  with  values  of  interest). 


The  second  entry  in  both  datafiles  is  the  measured  power  received  when  using  the 
standard  gain  horn  antenna  during  antenna  testing  as  the  antenna  under  test.  In  the 
Pasternack  datafile  these  values  have  been  adjusted  by  6.5  dB  because  the  measured 
values  included  an  additional  6.5  dB  of  loss  from  an  in-line  attenuator  that  was  required 
to  prevent  the  power  level  from  saturating  the  measurement  device  (Vector  Network 
Analyzer  or  VNA). 

The  third  entry  in  both  datafiles  is  the  measured  power  received  when  using  either 
the  X-band  or  Pasternack  horn  antenna  during  antenna  testing  as  the  antenna  under  test. 
The  maximum  relative  amplitude  as  measured  is  listed.  Then  this  value  is  subtracted 
from  the  corresponding  standard  gain  horn  measurement  value  to  get  a  delta.  Note  the 
value  is  subtracted  because  the  power  values  are  in  negative  dBs  and  the  X-band  and 
Pasternack  horn  values  were  less  than  the  standard  gain  horn  values.  This  delta  is  then 
subtracted  from  the  value  of  the  absolute  gain  of  the  reference  horn  antenna  obtained 
from  the  graph.  Note  the  values  are  subtracted  because  the  delta  though  positive  is  the 
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amount  below  in  power.  The  result  is  the  absolute  gain  for  that  antenna  at  the  given 
frequency. 

Then  the  value  for  the  first  measurement  to  exceed  3  dB  below  the  maximum 
relative  amplitude  is  listed  along  with  its  angle  in  degrees.  This  was  done  for  both  the 
first  value  above  and  below  in  angle  of  the  maximum  relative  amplitude  angle.  These 
angle  values  were  subtracted  from  one  another  to  find  the  3dB  Beamwidth. 

Example  calculation  from  X-band  Horn  #1  at  9.37  GHz  in  the  E  plane  (as 
described  above); 

Absolute  Gain; 

(-12.43115)-(-15.18115)  =  2.75  dB 
15.655-2.75  =  12.905  dB 

3dB  Beamwidth; 

17 -(-15)  =  32  dB 

F.  ANTENNA  ISOLATION  TEST  DATA 

At  80  cm  separation  from  each  other  oriented  main  beam  parallel  to  each  other 
both  at  horizontal  polarization  the  isolation  between  the  antennas  was  found  to  be  -42  dB. 
At  cross  polarization  (one  horizontal  and  one  vertical)  the  isolation  was  found  to  be  -57 
dB.  At  slant  45°  degree  polarization  (parallel,  then  perpendicular)  the  isolation  was 
found  to  be  -47  dB.  At  twice  the  separation  distance  (160  cm),  the  isolation  increased  by 
6  dB  (to  -48  dB  from  -42  dB  for  co-polarization  horizontal).  This  demonstrated  the  drop¬ 
off  from  spreading  loss  alone  (where  twice  the  distance  reduces  power  by  a  quarter).  A 
distance  of  6  ft  (same  as  lumber  length)  was  chosen  for  the  antenna  separation  and 
mounting  on  the  ACS.  This  was  sufficient  isolation  to  prevent  regenerative  feedback.  In 
addition  to  this  minimum  antenna  separation  distance  an  RE  isolator  was  used  to  make 
the  ACS  transmission  path  regenerative.  There  is  a  need  for  an  RE  isolator  because  all  of 
the  RE  components  are  bidirectional.  When  the  RE  amplifier  was  used  (as  in  the 
anechoic  chamber  testing)  the  isolator  was  not  required,  as  an  amplifier  is  a  uni¬ 
directional  device.  If  the  signal  path  through  the  ACS  is  allowed  to  be  bi-directional, 
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then  the  horns  beeome  radiators  in  a  two  element  array  with  large  eleetrieal  spaeing.  This 
results  in  a  radiation  pattern  with  numerous  narrow  lobes.  This  eondition  must  be 
eliminated  if  the  ACS  radiation  pattern  is  to  be  that  of  a  single  horn. 
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APPENDIX  G.  ANECHOIC  CHAMBER 


This  appendix  provides  information  on  the  setup  of  the  Naval  Postgraduate 
SehooTs  aneehoie  ehamber  on  the  roof  of  Spanagel  Hall.  This  is  the  aneehoie  ehamber 
owned  by  the  Eleetrieal  and  Computer  Engineering  Department.  This  was  the  only 
aneehoie  ehamber  used  in  this  thesis  (to  inelude  Antenna  Testing  and 
Amplitude/Refleetivity  Measurements). 
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Anechoic  Chamber  Calculations 


Distance  between  homs  -  19  ft.  about  18.5  ft  for  the  big  horn 


Path  loss  @  5.8GHz  =  10 log 


10 


4^' 

^  c  J 


=  63dB 


General  Friis  Equation 

=Pi+Gi+G^-Lp+  gains  -  losses 
In  the  chamber 

“  ^palh  ~  Variable  all  ~  pad 

P^^  =  P.-l0dB^,^^,^^-20dB^,  .s.SdB,„^-ldB^, 

The  chamber  measures  P^  -  P^^  or  the  system  loss. 
Solving  for  P^  and  substituting  into  the  equation  for  P^, 
after  a  little  algebra,  we  find  that: 


[p  -P  . 
y  r  ref 

measured  value 


-G,+L^+L 


variable  an 


-30dB 


MPS7300C 


Planned  configuration  after  repairs 
March  2006 


Figure  79.  Anechoic  Chamber  Diagram. 
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APPENDIX  H.  COMPACTRIO 


Figure  80.  Figure  of  CompactRIO  Software  Architecture  from  ni.com. 


A,  COMPACTRIO  REAL-TIME  EMBEDDED  CONTROLLER  (NI  CRIO- 

9004,  64/512  MEMORY) 

This  model  embedded  controller  runs  Lab  VIEW  Real-Time  with  64  MB  of 
nonvolatile  CompactFlash  storage  and  512  MB  of  DRAM  memory,  has  a  195  MFlz 
processor,  and  is  equipped  with  a  lO/lOOBaseT  Ethernet  port. 

Erom  the  CompactRIO  ADVISOR  on  the  webpage  for  this  device, 
http://sine.ni.eom/nips/cds/view/p/lang/en/nid/I4I55  “National  Instruments  CompactRIO 
real-time  embedded  controllers  offer  powerful  stand-alone  embedded  execution  for 
deterministic  EabVIEW  Real-Time  applications.  The  NI  CompactRIO  controllers  are 
designed  for  extreme  ruggedness,  reliability,  and  low  power  consumption,  with  dual  9  to 
35  VDC  supply  inputs  that  deliver  isolated  power  to  the  CompactRIO  chassis/modules 
and  a  -40  to  70  °C  operating  temperature  range.  The  NI  cRIO-90Ix  series  controllers 
feature  a  400  MHz  industrial  processor  balancing  low  power  consumption  and  powerful 
real-time  floating-point  signal  analysis,  control  and  logging.  The  cRIO-900x  controllers 
feature  a  200  MHz  industrial  processor  for  real-time  analysis,  control  and  logging.” 
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Figure  81.  Picture  of  cRIO-9004  from  ni.com. 
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CompactRIO  Real-Time  Embedded  Controllers 
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Owerwiew  and  Applications 

National  Instruments  cRIO-900x  real-time  embedded  controllers  offer 
powerful  stand-alone  embedded  execution  for  deterministic  LabVIEVV 
Real-Time  applications.  The  Nl  cRiO-9002  includes  32  MB  of  DRAM 
memory  and  C4  MB  of  nonvolatile  flash  storage  for  file  storage,  [he 
cRIO-9004  it^cluiles  64  MB  of  DRAM  memory  arid  51  ?  MB  of  ncnvulalile 
Hash  storage  lor  data-logging  applications.  Rolh  contreilers  arc  designed 
for  extreme  ruggedness,  reliability,  and  low  power  consumption  with 
dual  Oto  dBVDC  supply  inputs  that  deliver  isolated  power  to  the 
CompactRIO  chassis/mcdules  and  a  -dOtc  70  '"C  temperature  range, 

A 195  MHz  industrial  processor  balances  low  pcwiTr  consumption  with 
powertul  reaMimo  lloaling- point  signal  processing  and  analysis 
capabilihes  for  deterministic  cunlroi  loops  exceeding  1  kHz, 


figure  1.  CGnitj<ict.RiO  Harriwsre  Architecture 


System  Configuration 

Iho  CompaclRIO  real  time  controller  conncclsto  any  4  or  8  slot 
CornpaclRiO  reconfiyurabie  chassis.  Tlie  user-defined  FPGA  circuitry  in 
the  chassis  ccnlrols  each  I/O  module  and  passes  dala  lo  the  ccntrollor 
through  a  local  PCI  bus,  using  built-in  communication  functions. 


Embedded  Software 

You  can  synchronize  embedded  code  execution  to  an  FPGA-generated 
interrupt  request  (IRQ)  or  an  internal  millisecond  real-time  clock  source. 
The  LabVIFW  Real-Time  F.TS  OS  provides  reliatiilily  and  simplifies  tlie 
development  ol  complete  embedded  applical.Icns  [hat  include  l.ime-icritical 


ccnlrul  and  acquisition  loofis  in  additiun  to  lower-priority  loops  for 
postprocessing,  data  logging,  and  Ethernet/serial  ccmmunication. 

Built-in  elemental  I/O  functions  such  as  the  FPGy\  Read, /Write  function 
provide  a  ccmmunication  interface  to  the  highly  optimized  reconfigurabSe 
PHOA  circuitry.  Data  values  are  read  trem  the  PPb'A  in  integer  icrmal,  and 
llinii  corsverled  lo  scaled  Kigineeriny  units  in  the  ccntruiier. 
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Figure  82.  Datasheet  for  CompaetRIO  Real-Time  Embedded  Controllers  from 

ni.eom. 
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B. 


COMPACTRIO  RECONFIGURABLE  EMBEDDED  CHASSIS  (NI  CRIO- 
9101,  IM  GATE  4-SLOT  CHASSIS) 


From  the  CompactRIO  ADVISOR  on  the  webpage  for  this  device, 
http://sine.ni.eom/nips/cds/view/p/lang/en/nid/14154  “The  National  Instruments 
CompactRIO  embedded  chassis  are  the  heart  of  the  CompactRIO  system  because  they 
contain  the  reconfigurable  I/O  (RIO)  FPGA  core.  The  RIO  core  has  an  individual 
connection  to  each  I/O  module  and  is  programmed  using  easy  to  use  elemental  I/O 
functions  to  read  or  write  signal  information  from  each  module.  Because  there  is  no 
shared  communication  bus  between  the  RIO  FPGA  core  and  the  FO  modules, 
input/output  operations  on  each  module  can  be  precisely  synchronized  with  25  ns 
resolution.  The  RIO  core  can  perform  local  integer-based  signal  processing  and  decision 
making  and  directly  pass  signals  from  one  module  to  another.” 


Figure  83.  Picture  of  cRIO-9103  from  ni.com. 
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CompactRIO 
Reconfigurable  Chassis 


•  Design  hardware  using  Lab  VIEW 
programming  skills 

•  4  or  8-slot  chassis  for  any 
CompactRIO  I/O  modules 

•  1  M  or  3  M  gate  RIO  FPGA  core 
for  normal  or  extended  RIO 
processing  power 

•  DIN-raiL  mounting,  19  in.  rack 
mount,  and  panel  mounting  options 

•  NI  CompactRIO  Extreme  Industrial 
Certifications  and  Ratings'- 


•  Program  in  easy-to-use 
LabVIEW  FPGA  graphical 
development  environment 
to  automatically  synthesize 
an  optimized  high- 
performance  electrical 
circuit  implementation 

of  your  application 

•  RIO  FPGA  core  executes 
LabVIEW  control  logic  at 
rates  up  to  40  MS/s  using 
single-cycle  timed  loops 


Product 

■  'cRiOiaioi' 
cRIO-9102 


'  cRIO/9103' 
cRIO-9104 


FPfiA  System  Gates 

'liM 


RAM  (KB) 

;  ■  82'  ■  ■ 


DefaultTimebase(MHz} 


Maximum  Power  Consumption  (W) 


Built-In 

Panei  Mounting  Hoies 


Overview  and  Applications 

The  National  Instruments  CompactRIO  reconfigurable  chassis  are 
the  heart  of  the  CompactRIO  system  because  they  contain  the 
reconfigurable  I/O  (RIO)  core.  The  RIO  FPGA  core,  which  has  an 
indi-vidual  connection  to  each  I/O  module,  is  programmed  with  easy- 
to-use  elemental  I/O  functions  to  read  or  write  signal  information 
from  each  module.  Because  there  is  no  shared  communication  bus 
between  the  RIO  FPGA  core  and  the  I/O  modules,  I/O  operations  on 
each  module  can  be  precisely  synchronized  with  25  ns  resolution.  The 
RIO  core  can  perform  local  integer-based  signal  processing  and 
decision-making  and  directly  pass  signals  from  one  module  to  another. 

The  RIO  core  is  also  connected  to  the  CompactRIO  real-time 
controller  through  a  local  PCI  bus  interface.  The  real-time  controller 
can  retrieve  data  from  any  control  or  indicator  on  the  front-panel  of 
the  RIO  FPGA  application  throt^  an  easy-to-nse  FPGA  Read/Write 
function.  The  RIO  FPGA  can  also  generate  interrupt  requests  (IRQs) 
to  s-yndironize  the  real-time  software  execution  -with  the  RIO  FPGA. 
Typically,  the  real-time  controller  is  used  to  convert  the  integer  based 
I/O  data  to  scaled  floating-point  numbers.  In  addition,  the  real-time 
controller  typically  performs  single-point  control,  waveform  analysis, 
data  losing,  and  Ethemet/serial  communication. 

The  reconfigurable  chassis,  real-time  controller,  and  I/O  modules 
combine  to  create  a  complete  stand-alone  embedded  system. 
Application  development  consists  of  three  steps: 

1.  Target  the  reconfigurable  chassis  to  automatically  detect  the 
I/O  modules  and  develop  the  RIO  FPGA  application, 

2.  Compile  the  RIO  application  to  automatically  synthesize 

an  optimized  high-performance  electrical  circuit  implementation 
of  your  application, 

*See  CompactRIO  Overview  on  page  366  for  details. 


3.  Develop  the  real-time  controller  application  to  add 
floating-point  control,  signal  processing,  data  logging, 
and  communication 

Key  Features 

•  Create  any  local  or  multichassis  timing,  tri^ering,  and 
synchronization  scheme  with  25  ns  resolution 

•  Use  multiple  while  loops  to  create  a  parallel  processing  application 
for  high-performance  signal  processing  or  multirate  control  systems 

•  Bufit-in  PID  control  fimctions  for  control  system  loop  rates 
greater  than  100  kHz 

•  Generate  waveforms  or  implement  nonlinear  look-up  tables 
(LUTs)  using  LabVIEW  FPGA  express  Vis 

» Integrate  -widely  available  third-party  HDL  cores  using  the 
LabVIEW  FPGA  Module  HDL  Node 

» Enforce  critical  logic  and  interlocks  in  siiicon  hardware  drcuitry,  or 
use  the  parallel  RIO  architecture  to  create  duai,  triple,  or  quadruple 
redundant  systems 

Visit  m.com/cGmp&cJjiG  for  example  programs,  application  notes, 

and  other  developer  tools. 

For  ordering  information,  see  page  381. 


.®«  ®  d©-  <g>  g  C€  © 

h=?sRi;ii(  I  xsKjr?  (pending) 

For  more  information,  see  page  379. 

BUY  ONLINE  ot  CALL  '866.'  266-9891!  f :  ni.coni 
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Figure  84.  Datasheet  for  CompaetRIO  Reeonfigurable  Embedded  Chassis  from 

ni.eom. 
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C.  COMPACTRIO  I/O  MODULES  FOR  DIGITAL  INPUT  AND 
COUNTER/TIMER  (NI  9401,  8  CH  5  V/TTL  HIGH  SPEED 

BIDIRECTIONAL  DIGITAL  I/O  MODULE) 


From  the  CompactRIO  ADVISOR  on  the  webpage  for  this  device, 
http://sine.ni.eom/nips/cds/view/p/lang/en/nid/14168  “National  Instruments  CompactRIO 
digital  input  modules  are  compatible  with  any  CompactRIO  chassis  and  work  with  5 
V/TTL,  industrial  and  universal  logic  levels  and  signals  for  direct  connection  to  a  wide 
array  of  devices  such  as  industrial  switches,  encoders  and  transducers.  Each  module 
features  an  integrated  connector  junction  box  with  a  screw  terminal  or  industrial  D-Sub 
connector  for  flexible  and  low  cost  signal  wiring.  All  modules  feature  the  NI 
CompactRIO  Extreme  Industrial  Certifications  and  Ratings.  Each  CompactRIO  EO 
module  has  a  direct  connection  to  the  reconfigurable  EPGA  in  CompactRIO  systems. 
Using  the  EabVIEW  EPGA  module,  users  can  configure  and  reconfigure  all  CompactRIO 
digital  input  modules  to  implement  counter/timer,  encoder,  digital  protocol,  and/or  static 
digital  inputs.” 


Eigure  85.  Picture  of  cRIO-9401  from  ni.com. 
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Owerwiew 

High-pert orniancu  digital  output  and  switching  modules  tor  National 
Instruments  CcmpactRIQ  embedded  systems,  R  Series  expansion 
chassis,  and  Nl  CompactDAQ  systems  provide  e.xteiided  voltage  ranges 
and  high  current-switching  capacityfor  direct  control  of  a  vvide  array  of 
induslrial  and  aulomutive  aclualcrs.  Faoh  rnodulB  leaUires  an  inlegraled 
connector  junction  box  with  screvv-lerminal  or  cable  options  for  flexible, 
luw-cusl  signal  wiring.  All  modules  feature  [he  Nl  CompaclRSQ  Fxdreme 
Industrial  Certifications  and  Ratings. 

System  Compatibility 

You  can  use  N!  C  Series  modules  interchangeably  in  CcmpactRIO  and 
Nl  CompactD.AQ  chassis.  Many  of  the  advanced  timing  features 
described  apply  only  I o  CompaciRIO  reconfigurable  l/'O  systems  and  not 
to  Nl  CampaclDAQ, 

Advanced  Features 

When  used  in  CcmpactRIO,  C  Series  digital  output  modules  connect 
directly  lu  rBcunfigyrablG  I/O  (RIO)  FPGA  hardware  to  create  high- 
performance  embedded  systems.  The  reconfigurable  FPG.A  hardware 
within  CcmpactRIO  provides  a  variety  of  options  for  timing,  triggering, 
synchronization,  digital  waveform  generation,  or  digital  communication. 
For  instance,  with  CompaciRIO  you  can  implement  a  circuit  lo  generate 


pulse-widtii  mudulalion  (PWM)  outputs  for  controlling  mclors,  liealers, 
or  tans,  as  well  as  perform  pulse  code  modulation  encoding  (PCMF.)  for 
wireless  telemetry  applications. 

Key  Features 

•  High-performance  digital  output  switching  for  any  CompactRiO  embedded 
system,  R  Series  e.xpansion  chassis,  or  Nl  CompactD.AQ  system 
Screw'-terminal,  strain  relief,  high-voltage,  cable,  solder-cup 
backsholl,  and  other  connoctivity  options 

•  Channrd-to-eartli  ground  double-isolation  barrier  for  safely 
and  noise  immunity 

•  Nl  CompactRiO  Extreme  Industrial  Certifications  and  Ratings 

Visit  iii.cofii/compactrw  or  for  up-to-date 

information  on  module  availability,  example  programs,  application  notes, 
and  other  developer  tools. 

Typical  certificatiuns  -  Actual  specincations  vary  from  product  to 
product.  Visit  niconi/certificatbii  for  details. 
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Figure  86.  Datasheet  for  C  Series  Digital  Output  and  Relay  Modules  from  ni.com. 
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D. 


NI  DEVELOPER  SUITE  WITH  REAL-TIME  AND  FPGA  SOFTWARE 
FOR  COMPACTRIO 


From  the  CompactRIO  ADVISOR  on  the  webpage  for  this  software  paekage, 
http ://sine .ni. com/nips/eds/view/p/lang/ en/nid/ 1 4294  “The  NI  Developer  Suite  Core 
Package  with  the  Real-Time  Deployment  and  FPGA  Deployment  options  include  all  of 
the  tools  you  need  for  building  a  custom  embedded  system  -  Lab  VIEW  FPGA  for 
customizing  hardware  I/O  in  the  user-configurable  FPGA  core  and  Lab  VIEW  Real-Time 
for  building  deterministic  and  reliable  real-time  applications.  The  NI  Developer  Suite 
Core  Package  with  the  Real-Time  Deployment  and  EPGA  Deployment  options  also 
include  a  comprehensive  set  of  Lab  VIEW  add-on  tools  for  algorithm  development,  real¬ 
time  development  analysis,  Microsoft  Office  report  generation,  and  database 
communication.” 

Of  particular  note  in  this  software  package  is  the  EPGA  Deployment  Option, 
which  includes  the  LabVIEW  EPGA  Module  and  NI  RIO.  These  are  essential 
programming  components  to  the  LabVIEW  software  library  used  with  the  ACS.  NI  RIO 
is  needed  to  communicate  with  the  CompactRIO  Chassis(EPGAs)  and  the  LabVIEW 
EPGA  module  for  programming  commands  that  can  be  auto-compiled  into  VHDL  for  the 
CompactRIO  Chassis(PPGAs). 
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APPENDIX  I.  LAB  VIEW  FOR  THE  DIGITAL  PHASE  SHIFTER 


How  the  digital  phase  shifter  is  controlled  by  code  written  in  Lab  VIEW  and  run  on  the 
CRIO  is  explained  here. 

The  Lab  VIEW  Virtual  Instruments  (Vis)  described  here  include:  ManualControl, 
PhaseShifterControl,  MasterControl,  and  SlaveControl.  Together  they  provide  three 
distinction  functionalities  in  the  development  of  the  ACS’s  software.  These  are  the  main 
Vis.  All  Vis,  including  Block  Diagram  and  Eront  Panel  files  for  each,  are  included  after 
the  description  of  these  main  Vis.  The  Vis  not  described  are  test  components  only  and 
not  used  with  the  ACS’s  current  functionalities,  although  they  did  play  a  part  in  its  past 
development  and  possibly  future  development  and  are  included  here  for  that  purpose. 
Primarily  they  were  used  to  test  certain  functions  that  were  later  incorporated  into  the 
main  Vis  or  demonstrate  alternative  ways  of  implementing  the  same  functionality  as  was 
incorporated  into  the  main  Vis  as  part  of  testing  method  efficiency.  These  Vis  include: 
datalO,  fpgaRX_test,  fpgaTX_test,  NetworkRX_test,  NetworkTX_test,  and  tDcalc. 
KNORR-CAL_project  is  an  overall  view  of  the  Lab  VIEW  project  with  all  Vis  and 
hardware  indicated.  The  VI  PSC6-bit  is  merely  an  example  of  how  to  modify  the  files 
for  a  different  number  of  control  bits.  When  the  ACS  was  first  developed  a  6-bit  digital 
phase  shifter  rather  than  an  8-bit  digital  phase  shifter  was  used.  The  basic  principle  of 
this  modification  is  to  globally  find  and  replace  the  number  representing  2^(current 
number  of  bits)  with  2^(future  number  of  bits).  Eor  example,  replace  255  (for  8-bit)  with 
63  (for  6-bit). 

Observe  the  Block  Diagram  for  each  VI  and  follow  along  in  the  diagram  left  to  right  as 
the  description  explains  the  graphical  code. 

A,  PHASESHIFTERCONTROL 

First,  the  count  is  reset  (initialized)  to  zero.  Then,  the  time  delay  (tD)  is  read  in  to 
become  the  time  between  phase  state  changes.  This  is  also  referred  to  as  ‘ticks’.  1  tick  = 
1  bit  change  =  25  ns  (the  time  of  one  40MHz  clock  cycle).  The  count  is  checked  so  as  to 
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not  exceed  the  number  of  phase  shifter  bits  (phase  states).  Here  there  are  8-bits  or 
2^8=256  phase  states.  Since  0  is  considered  the  first  phase  state,  the  count  will  increment 
through  all  256  phase  states  and  stop  at  255. 

If  the  Doppler  shift  is  going  to  be  an  up  Doppler  shift  (increase  in  frequency)  the 
count  proceeds  unaltered,  counting  up.  If  the  Doppler  shift  is  to  be  a  down  Doppler  shift 
(decrease  in  frequency)  the  count  proceeds  inverted,  counting  down.  The  integer  value  of 
the  count  is  then  converted  into  a  Boolean  value.  For  example,  the  integer  value  of  157  = 
10011101.  The  Boolean  value  is  then  broken  into  each  of  its  individual  bits  and  routed  to 
the  corresponding  Digital  I/O  (DIO)  pin  to  apply  or  remove  voltage  from  the  appropriate 
phase  shifter  control  pin.  In  this  Vi’s  current  state,  both  DIO  modules  receive  the  same 
signal. 

See  Chapter  on  Future  Work.  Once  the  programmable  attenuator  is  received,  the 
second  group  of  DIO  pins  will  be  assigned  to  control  this  device  with  dependent 
statistics.  For  independent  statistics,  a  second  PhaseShifterControl  VI  will  need  to  be 
created( duplicated  and  renamed  AttenuatorControl).  The  only  difference  between  the 
original  and  its  duplicate  is  that  the  duplicate  DIOs  will  be  deleted  one  from  each  VI  so 
that  each  VI  contains  only  one  unique  set  of  DIO  assignments.  For  dependent 
(correlated)  statistics  the  DIOs  will  be  controlled  from  within  the  same  VI  in  a 
synchronized  time  delay  (tD). 

Once  the  commands  for  the  appropriate  voltages  have  been  sent,  the  count 
increments,  and  the  process  begins  all  over  again.  This  process  continues  until  one 
complete  cycle  through  all  the  phase  states  is  completed  (one  period  of  2n).  Then,  the 
loop  encounters  an  exit  condition  and  the  count  is  reset  and  the  whole  sequence  begins 
again  (the  next  period  of  2n). 

B,  SLAVECONTROL 

Paired  with  MasterControl  provides  automatic  PhaseShifterControl  to  work  with 

data  read  in  from  a  file.  The  difference  is  that  an  array  of  integer  values  corresponding  to 

the  tD  is  first  read  and  stored.  Then,  the  individual  time  delay  integer  values  are  used  to 

provide  the  appropriate  phase  shifter  Doppler  shift  values  as  opposed  to  the  manual  entry 

130 


input  by  a  user  in  the  GUI  of  PhaseShifterControl.  In  addition,  instead  of  a  switeh  being 
flipped  by  a  user  in  the  GUI  to  deeide  up  or  down  Doppler  shift,  the  sign  of  the  integer 
value  is  used.  A  negative  sign  is  used  to  indieate  a  down  shift  and  the  absenee  of  a 
negative  sign  (i.e.  positive)  to  indieate  an  up  shift.  Zero  indieates  no  Doppler  shift.  The 
proeess  is  then  exaetly  the  same  as  that  implemented/programmed  in  PhaseShifterControl 
(See  seetion  above  by  the  same  name). 

*NOTE;  the  number  of  Doppler  frequeneies  eannot  be  automatieally  deteeted  from  the 
data  file  stored  on  the  PC.  The  number  of  frequeneies  must  be  defined  from  within  the 
VI.  This  allows  the  appropriate  sized  array  to  be  ereated  in  the  software  even  before  it 
starts  running.  This  is  noted  in  the  VI  Front  Panel  (GUI)  as  well. 

C.  MASTERCONTROL 

This  VI  reads  in  the  value  of  the  Doppler  frequeneies  from  a  tab  delaminated  file 
(sueh  as  how  the  array  of  variables  would  be  exported/saved  from  the  Weather  Signal 
Simulator  in  MATLAB).  The  numbers(frequeneies)  are  read  in  as  double  floating  point. 
numtDs  is  the  number  of  Doppler  frequeneies  read  in.  The  Doppler  frequeneies  (DF)  are 
then  used  to  ealeulate  the  Radial  Veloeity  (Vr),  Period  (T),  and  Bit  Delay  (t)  aeeording  to 
the  following  relations: 

Q 

Equation  1  Vr  =  DF  x -  in  meters  per  seeond 

RF  x2 

Equation  2  T  =  — x  1 000  in  milliseeonds 
DF 

T  lO*” 

Equation  3  t  = - in  ‘tieks’  (of  25ns) 

n  25 

Where  the  RF  Frequeney  (RF)  is  defined  to  be  9990000000Hz  or  9.99GHz,  the  Speed  of 
Fight  (e)  is  299792458  m/s,  and  the  number  of  Phase  Steps  (n)  is  63  (for  6  bits). 

However,  sinee  only  integers  (INT)  ean  be  passed  to  the  FPGA,  the  doubles 
(DBF)  are  eonverted  to  integers  and  sent  to  the  FPGA.  The  Period,  DF,  and  Vr  are  then 
reealeulated  based  on  these  16-bit  integer  values.  This  is  to  show  the  user  what  aetual 
values  will  be  produeed  by  the  ACS.  This  only  affeets  values  stored  in  the  data  file 
whieh  may  have  been  truneated  by  the  eonversion.  This  is  a  potential  software  limitation. 
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However,  this  will  not  affeet  the  ealibration  process,  only  “real  world”  weather 
simulation,  because  any  variation  from  the  input  values  that  have  been  truncated  are  now 
recomputed  and  displayed  to  the  user.  The  radar  should  be  tested/calibrated  to  agree  with 
these  values,  not  the  values  output  by  the  Weather  Signal  Simulator.  This  is  a  limitation 
of  the  calibrator’s  software  when  interacting  with  the  FPGA.  Although  values  may  be 
calculated  accurately  by  the  Weather  Signal  Simulator  to  agree  with  “real  world”  weather 
patterns,  the  calibrator  cannot  always  reproduce  these  exact  values.  This  is  a  hardware 
limitation  due  to  the  FPGA  and  seen  in  the  software.  In  practice,  this  limitation  is  similar 
to  other  hardware  limitations  already  discussed.  See  Hardware  Limitations  Section. 

*NOTE:  before  executing  the  VI,  the  array  size  in  SlaveControl  must  be  defined  to  the 
number  of  Doppler  frequencies  intended  to  be  read  from  the  data  file.  This  is  not 
automatic.  It  cannot  predetermine  the  number  of  entries.  There  is  no  “read  until  end  of 
file.”  See  also  SlaveControl  NOTE. 

*NOTE;  the  path  on  the  PC  to  the  data  file  must  be  entered  into  the  VI.  The  data  file 
must  be  a  single  row  of  tab  delaminated  Doppler  frequency  entries  stored  in  a  text  file 
(.txt)  where  a  negative  sign  indicates  a  down  Doppler  shift  and  no  sign  indicates  an  up 
Doppler  shift.  Do  not  use  +  signs. 
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LabView  Code  (#  indicates  state  change) 
User  Interface 


CODE 

Ul 


datalO 


reads  in  data  file  containing  Doppler  Shift  values  input  from 
Ulf  Shroder's  MATLAB  weather  simulator 


fpgaRX_test  test  file  to  send  data  read  on  laptop  to  fpga 

fpgaTX_test  test  file  to  receive  data  sent  from  laptop  on  fpga 

KNORR-CAL  project  overall  software  view.  LabView  Project. 

vi  to  manually  control  individual  phase  shifts.  High  or  low 
any  one  or  combination  of  bits 

vi  to  automatically  control  phase  shifts.  Runs  on  laptop, 
test  file  to  send  data  read  on  laptop  to  CRIO 
test  file  to  receive  data  sent  from  laptop  on  CRIO 
only  difference  from  PhaseShifterControl  is  6-bit  (63)  test 
PhaseShifterControl6-bit  vs.  8-bit  (255) 

fpga  resident  to  control  phaseshifter  through  a  series  of 
PhaseShifterControl  phase  shifts 

SlaveControl  version  of  PhaseShifterControl  to  work  with  MasterControl. 

LabVIEW  computer  for  the  mathematical  relations  between 
time  delay  tD,  radial  velocity,  and  Doppler  frequency. 
MATLAB  Script  Node  is  used  in  MasterControl  instead  of 

tDcalc  this. 


ManualControl 
MasterControl 
NetworkRXtest 
NetworkTX  test 


Table  10.  Table  Summary  of  Lab  VIEW  Vis. 
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Figure  87.  datalO  CODE 
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Figure  88.  datalO  UI 


Figure  89.  fpgaRX  test  CODE 
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fpgaRX test.vi  Front  Panel  on  FPGA  Target 


Figure  90.  fpgaRX  test  UI 


Figure  9 1 .  fpgaXX  test  CODE 


Figure  92.  fpgaTX_test_UI 
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1  iK  Project  Explorer  -  KNORR-CAL.Ivproj 

i 

1 

Fiie  Edit  View  Project  Operate  look  Window  Help 

jj  “b 

X 

■ 

B  Prolect:  KNORR-CAL.Ivpro1 

$  B  cRBn 

I  MasterControl.vi 

I  I  datalO.vi 

I  fpgaTX_test.vi 

I  f  NetvvorkTXJ:est.vi 

I  host.vi 

I  f- Dependencies 
^  ’iti  Build  Specifications 
6  ^  CRIO-KNORR  (169.254.0.8) 

$  H  FPGA  Target  (RIOO::INSTR,  cRIO-9101) 
I  B  Digital  Line  Input  and  Output 
M  B  ^  NI9401 

1  hi*  DIO0(NI9401) 
i  i*  DIO  1(NI  9401) 

I  hi*  DIO2(NI9401) 

I  i*  DIO3(NI9401) 

I  hi*  DIO4(NI9401) 
h  i*  DIOS  (NI 9401) 

I  hi*  DIO6(NI9401) 

'  i*  DI07  (NI  9401) 
i  B  Q  NI 9401 2 

1  i*  DIOO  2  (NI  9401  2/DIOO) 
h  i*  DIOl  2  (NI  9401  2/DIOl) 

I  i*  DI02  2  (NI  9401  2/DI02) 
h  i*  DI03  2  (NI  9401  2/DI03) 
h  i*  DI04  2  (NI  9401  2/DI04) 
h  i*  DI05  2  (NI  9401  2/DI05) 

I  i*  DI06  2  (NI  9401  2/DI06) 

-  i*  DI07  2  (NI  9401  2/DI07) 

I  j!fe  40  MHz  Onboard  Clock 
I  h  si’  PhaseShifterControLvi 
h  l?ii  PhaseShifterControl6-bit.vi 
i  r  SlaveControl.vi 
h  isi  ManualControl.vi 
i  1  si  fpgaRX_test.vi 
!  NI9401(Stot2) 

h  NI  9401  2  (Sbt  3) 
h  ^  Dependencies 
'  fc  Build  Specifications 
I  ■  si  tDcat.vi 
hsi  Net\vorkRX_test.vi 
i  target  -  single  rate  -  variables.vi 
f  1^  support  -  acquire  data.vi 
variables  -  network  -  RT.Ivlib 
i  hJ^  Array 
!■  SfiJ  data 
i  ’S;:  stop 
i  ^  Dependencies 
'  iti  Build  Specifications 


Figure  93.  KNORR-CAL_project 
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Figure  94.  ManualControl  CODE 
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^  ManualControLvi  Front  Panel  on  FPGA  Target  ^ 


Figure  95.  ManualControl  UI 


Figure  96.  MasterControl  CODE 
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Figure  97.  MasterControl  UI 


Figure  98.  NetworkRX_test_CODE 
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Figure  99.  NetworkRX_test_UI 


Figure  100.  NetworkTX_test_CODE 
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Figure  101.  PhaseShifterControl6-bit_CODEl 


Figure  102.  PhaseShifterControl6-bit_CODE2 
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Figure  103.  PhaseShifterControl6-bit_UI 


Figure  104.  PhaseShifterControl  CODEl 
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Figure  105.  PhaseShifterControl_CODE2 


Figure  106.  PhaseShifterControl  UI 
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Figure  107.  SlaveControl  CODEO 


Figure  108.  SlaveControl  CODEl 
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Figure  109.  SlaveControl_CODE2 


Figure  110.  SlaveControl  Ul 
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Figure  111.  tDcalc  CODE 
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Figure  112.  tDcalc  UI 
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